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VORTRAGSBERICHTE 


A. Strömungslehre 


Die Reflexion schwacher Störungen an Unstetigkeitsflächen 
einer ebenen Unterschallströmung 


Von Fr. Vandrey in Göttingen 


Der Streifen 0 < y.< y, einer Strömungsebene sei nach oben abgeschlossen durch eine 
gerade feste Wand y = y,, nach unten durch eine schwach kurvenförmige, gleichfalls feste Wand 
y= y(z) <y,. In diesem Streifen ströme in x-Richtung ein Gas mit der Gesehwindigkeit U, 
im unteren Teil 0 <y <y, und mit der Geschwindigkeit U, im oberen Teil y, <y < y,, so 
daß die Grundströmung bei y= y, eine Unstetigkeitsfläche enthält. 

Betrachtet man die beiden Grenzfälle U, = U,und U, = 0, so hat man im ersten einen 
oben durch eine feste Wand Yy= y, begrenzten Streifen, im zweiten einen oben durch eine freie 
Grenze y= y, begrenzten Streifen. Da nun die Reflexion schwacher Störungen an festen und 
freien Grenzen mit verschiedenem Vorzeichen erfolgt, kann man erwarten, daß die Einflüsse 
von der oberen Wand und von der Unstetigkeitsfläche sich im unteren Teilstreifen entgegen- 
wirken und sich teilweise, in besonders günstigen Fällen sogar völlig aufheben können, wenn U, 
zwischen 0 und D, liegt. 

Zur Berechnung der Strömung in den beiden Teilstreifen dienen die folgenden Randbedin- 
gungen: An der oberen und an der unteren Grenze muß die Normalgeschwindigkeit verschwinden, 
ferner muß auf der Unstetigkeitsfläche gleichfalls die Normalgeschwindigkeit auf beiden Seiten 
Null und der Druck in jedem ihrer Punkte auf beiden Seiten gleich sein. Die ersten drei Bedin- 
gungen lassen sich durch die Gleichungen ausdrücken 


en dt, u, = dydelüry=0O: v,/0,=%,10, fürys=y: 
Die noch übrige Druckbedingung erhält man aus der Druckgleichung wdw +dp/oe=0 zu: 


U, d,,0 = U1d,, 0, für y= Y,. Diese Gleichung läßt sich mit Hilfe der für alle -Gase gültigen 
Beziehung c? = x : p/o umformen in 


Dyt„le = Vıv,li für y=Y- 


Außerdem muß das Potential der Strömung im Inneren der beiden Teilstreifen der linearisierten 
kompressiblen Potentialgleichung genügen. 

= Eine Grundlösung der Aufgabe erhält man, wenn man y(x) als einfache Sinuswelle voraus- 
setzt, aus dieser Grundlösung kann man dann die Lösung für beliebiges y(&) durch Überlagerung 
in Form Fourierscher Reihen oder Integrale gewinnen. Für y(2) = a, sin q% kann man für 
das Potential ® den Ansatz machen 
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wo M, — U,jc, die Machsche Zahl ist. Der Ansatz erfüllt die Potentialgleichung und für 


Ys <a, was bei genügender Breite des Streifens der Fall ist, auch die Randbedingung bei 
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a=0. Aus den abiigen drei e Randbedingungen eat man da 
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folgt. a* ist ein Maß für den Einfluß der a Wand und der Unstetigkeitsläche an der 
unteren Grenze y= yie): Dieser Einfluß verschwindet, wenn 


era) — At 1-2 
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ist. Ein Diagramm (Kurven q(Y,— %,) = const in der M,, M,-Ebene) zeigt, daß sich für diesen | 2 
Fall der Einfluß der Kompressibilität-dahin auswirkt, daß die Geschwindigkeit U, im oberen 
Teilstreifen größer sein muß als bei inkompressibler Strömung und daß sie sich mit ENTE we 


von U, an die Schallgeschwindigkeit gleichfalls dieser nähert. 


Näherungsyerfahren zur Berechnung kompressibler, 
Unterschallströmung 


Von E. Krahn in Göttingen 


‘ Für die Berechnung der kompressiblen Strömung um einen gegebenen Körper sind bisher 


keine exakten Lösungen bekannt. Über die Prandtlsche Regel hinausgehende Näherungs- 
methoden sind meist mit sehr umfangreichen Rechnungen verbunden und sind nicht immer 


. bei der Umströmung beliebig vorgeg&bener Körper durchführbar. Hier wird ein Verfahren be- 


handelt, nach dem zwei Näherungen einer zweidimensionalen Umströmung eines zylindrischen 
Körpers bis zu numerischen Resultaten durchführbar sind, und auf Grund vieler gerechneter 
Beispiele dann eine Näherungsformel für die Geschwindigkeit an der Körperoberfläche gegeben. 


. . In der «-, y-Ebene sei eine Kontur vorgegeben, die parallel der #-Achse angeströmt wird. 
Es wird angenommen, daß die Geschwindigkeit ® der inkompressiblen Strömung bekannt ist. 
Der Geschwindigkeitsvektor der kompressiblen Strömung soll mit p bezeichnet werden. Als 
erste Näherung wird 
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gesetzt, wo v, und V„ die Anströmgeschwindigkeiten der kompressiblen und inkompressiblen 
Strömung sind, o und og, die örtliche und die der Anströmung entsprechende Dichte. 
Um eine zweite Näherung!) zu berechnen, wird ein Korrekturglied d’ durch folgende 
Gleichung eingeführt 
p’ Se R -H vd’ 
Yan 5 Va 


Setzt man diesen Ausdruck in die Gleichungen der Kontinuität und der. Wirbelfreiheit, so 
findet man 


1 4 
divb’ = na (®-+v’) . gradlog 0; rtv — 5 


Die #- und y-Komponente des Vektors p’ seien v, und v,und die von®+b’ V, und P,. 
Wir führen nun neue Variable durch die folgenden Gleichungen ein 


a = Fr hPa 
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(®+v’) x gradlogo. 


sHtiy=2 a —iy=2 ———v (2%) 
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!) Diese ist von A. Betz und dem Vortragenden gemeinsam durchgeführt worden. 
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läßt sich. 
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Die Gleichung wird nun zur Berechnung einer zweiten Näherung für die Geschwindigkeit - 
‚ benutzt, indem auf der rechten Seite v’ gegenüber ® vernachlässigt und für V die Geschwindig- Fi 
keit der inkompressiblen Strömung gesetzt wird, von der wir wissen, daß sie eine Funktion 


® von 2 allein ist., Man hat dann GR 
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wo Fiz ) so zu bestimmen ist, daß den Randbedingungen genügt wird. Berücksichtigt man die $ 


- Entwicklung von A, so hat man Ausdrücke dr Form 0.0.2.0... 
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Ri zu berechnen. Bedeutet d den Winkel der Profiltangente mit der x-Achse und s.die Bogenlänge 
längs der Kontur, so wird $r BR 
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wo R der Krümmungsradius des Profiles an der betreffenden Stelle ist. In Real- und Imaginär- 
teil des gefundenen Ausdruckes haben wir nun die Tangential- und Normalkomponente des zu 
2... , diesem ’n gehörigen Korzekturanteils der Geschwindigkeit. Die Normalkomponente bringt man, 
durch Hinzufügen einer Funktion von 2 zum Verschwinden. Unter Zuhilfenahme der Abbil- 
- dungsfunktion des Profiles auf den Kreis führt man die Rechnung in der Kreisebene durch, 
indem’ man die Normalkomponente nach dem Kreiswinkel in eine trigonometrische Reihe ent- 
wickelt. Dadurch, daß maıt die Normalkomponente zum-Verschwinden bringt, erhält man dann 
zusätzlich eine Tangentialkomponente, die dann wieder auf das Profil umgerechnet wird. TEE 
Eine Zirkulation um das Profil kann mit berücksichtigt werden und ergibt keine weiteren . er 
4 Rechenschwierigkeiten. | Nr EHER DE 
Nach der dargestellten Methode kann auch die Geschwindigkeit im Felde außerhalb der 
umströmten Kontur berechnet werden. 3“ 
Nach dem hier besähriebenen Rechenverfahren ist die Geschwindigkeit an der Oberfläche 
“ einer ganzen Reiche von Profilen berechnet worden, unter anderen für elliptische Zylinder vom 
e Dickenverhältnis 0,2 bis 1 für verschiedene Machsche Zahlen der Anströmung bis zum Er- 
reichen der Schallgeschwindigkeit im Strömungsfeld. Die Ergebnisse dieser numerischen Rech- 
‘nungen können in folgender empirischen Näherungsformel zusammengefaßt werden: 
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w=W m 0,44 Ma&° W2 (W2-—-1)?—0,4Ma, = -,) (1 — 7) w°»5 aM ; 
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wo w und W den Betrag der auf die Anströmgeschwindigkeit bezogenen kompressiblen und 
inkompressiblen Geschwindigkeit an der Oberfläche des Profils bedeuten, Ma, die Machsche 
3 - Zahl der Anströmung und ! die Tiefe des Profiles ist. x ist die von der Vorderkante nach der 
| Hinterkante laufende Koordinate und ö die Entfernung in «-Richtung des Punktes der maxi- 


malen Geschwindigkeit von der Vorderkante. 
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gende neue Veränderliche ein: 


Cesekmanlgkeiteh der s cha. 

Von Klaus Oswatitsch in x 
Im Bienen wird die ebene, stationäre Umströmung flacher Droieb 

schwindigkeiten behandelt. Diese kann nach Guderley [1] stets als d 


werden. Sie wird durch- die bekannten Gleichungen für die Drehur 
nuität beschrieben. Führt man mit ß als FRARCNAEN, .c* als kritische 
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(u, = Kucbröihgeschwändigken)‘ so lauten für ausreichend schlanke Bine die euer 38 
und Randbedingungen: : 
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_ Für nn werden die Definitionen zweckentsprechend etwas geändert: 
ß®=M2—1, U= —=, usw., woraus ein a en Gleichungssystem folgt. 

Nach Gl. (1) ae ee Randbedingungen in Gl. (2) unverändert, wenn das Profil affin 
verdickt und die Anströmgeschwindigkeit gleichzeitig entsprechend verringert wird. Aus Gl. (2), 
welche auch für lokale Überschallgebiete gilt, folgt also ein Ähnlichkeitsgesetz, welches die 
Strömung um ein Profil auf die Strömung um ein flacheres Profil bei höherer Machzahl der 
Anströmung abbildet. Dabei ändern sich die Störgeschwindigkeiten wie 7”?’? (7 = Dickenver- 
hältnis). Das Ähnlichkeitsgesetz gilt entsprechend bei Überschallströmung. Verdichtungsstöße 
werden nur richtig innerhalb des Bereiches, in dem die Verluste mit der dritten Potenz der 
Druckanstiege wachsen, wiedergegeben. Das Gesetz wurde zuerst von G. Guderley [1] er, 
wähnt. Unabhängig davon wurde es in [2] ausführlich abgeleitet und angewendet (kritische 
Machzahl). Inzwischen wurde es auch von v. Kärmän im Journal of Aeronautical Sciences. 
kurz wiedergegeben. Für Schallanströmung (u, = c*) wurde ein entsprechendes Gesetz von 
v. Kärmän in Paris (1946) vorgetragen. 

Dem Gesetz entspricht eine Näherung der Charakteristik im Hodographen am Schallkreis 
durch eine Neilsche Parabel [3]. Es ergibt sich daraus eine stärkere Winkelabhängigkeit des 
Druckes als bei den Näherungen, welche den Auftriebs- und Widerstandsformeln von Ackeret 
und Busemann zugrunde liegen. Diese ergeben also zu niedrige Werte in Schallnähe, wie die 
Prandtlsche Regel zu niedrige Werte bei schallnahem Unterschall ergibt. Die verbesserten 
Formeln (Krahn oder v. Kärmän-Tsien) ergeben richtigen Auftrieb, aber bei lokalem Über- 
schallgebiet noch keinen Widerstand, 

Ein Vergleich dieser Formeln mit den gemessenen Werten zeigt richtige, aber zu weit 
stromauf Tiegende Maximalwerte der Geschwindigkeit. Das Stromabwirken der Überschall- 
strömung kann grob angegeben werden, indem Höhe des Überschallgebietes, mittlerer Mach- 
scher Winkel und Profilkrümmung abgeschätzt und durch den mittleren Druckbeiwert © Op; den 
kritischen Druckbeiwert en und den maximalen Druckbeiwert bei inkompressibler Strömung 


(Ma = 0) 6,,max Ausgedriickt wird. Für den Wert c, ergibt sich dann folgende Abwärts- 
verschiebung As (t Flügeltiefe) [4]: 


da Ras 4 k+1— 6) Ya —% 
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Der mit diesen Verschiebungen gerechnete Widerstandsanstieg stimmt mit den Versuchen 
überraschend gut überein. Ein Anstieg ergibt sich erst dann, wenn sich das Überschallgebiet 
über stärkere Neigungen der Profiloberfläche erstreckt. Für die Beschreibung des dann ein- 
tretenden starken Widerstandsanstieges ist aber offenbar nur eine rohe Formel erforderlich. 


Das Ähnlichkeitsgesetz gestattet die Übertragung‘ der Widerstandskurve auf verschiedene 
Profildicken. 
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; Über die Berechnung. ar ee) von Stö ränben‘ in ee x a 
_ Überschallströ ömung nach der direkten Hodographenmethode 
Von Manfred Schäfer (Max-Planck- Institut für Strömungsforschung, Göttingen) x 


m einer reibungsfreien Überschallströmung pflanzen sich bekanntlich Singularitäten a; 
Ba: le, die der Strömung etwa durch gewisse Bedingungen an einer festen Wand aufgezwungen ES. 
werden, längs der „‚charakteristischen Grundkurven“ (Machschen Wellen) fort. ‚Obwohlesdaher 
' naheliegend erscheint, für die Berechnung der Ausbreitung von Störungen diesen anschaulichen. E 
 Physikalischen Sachverhalt, daß die Machschen Wellen Träger der Singularitäten sind, un- 5 
mittelbar zum Ausgangspunkt zu machen und sich dabei des direkten Zusammenhanges zwischen 
ei Geschwindigkeitskomponenten und Ortskoordinaten zu bedienen (‚direkte Hodographen- Sr 
methode“), sind Ansätze in dieser Richtung! bisher nicht bekannt geworden. Für die ebene Zirku- 
lationsströmung hat G. I. Taylor!) vor nahezu 20 Jahren die genannte Fragestellung nach der 
Zeig „Methode der kleinen Schwingungen‘ behandelt, aus neuerer Zeit liegen von G. Gu derley? = 
_ Ansätze für eine Methode vor, die sich auf eine Entwicklung i in der Umgebung eines ausgezeich- = 

I neten Punktes des Strömungsfeldes gründet. . 

Be .9 Setztuman voraus, daß die eingeführten polaren sa eek ynd 9 
Br; LP Geschwindigkeitsbetrag, © Stromlinienwinkel gegen eine feste Richtung) bei Überquerung 
der Machschen Welle 7 = const. stetig bleiben (man schließt also sprunghafte Richtungsände- 
rungen aus, wie sie z. B. an Stoßfronten auftreten), so bleiben die zugelassenen endlichen Sprünge 
i [9,] i in den Ableitungen von © quer zur Machschen Welle längs dieser konstant, etwa gleich A, 
- » wie kürzlich W. Tollmien?) nachgewiesen hat. Für die Rechnung von Bedeutung ist eine all- 
. „gemeine Beziehung zwischen dem Krümmungssprung einer Stromlinie (speziell der festen Wand) 
- und den zwei Sprunggrößen A, und A, bezüglich der beiden die Stromlinie treffenden Machschen 
- Wellen. Soweit ist alles ganz streng. Für schwache Störungen liegt nun folgende Näherungs- 
methode nahe, bei der die Störungen (gekennzeichnet durch Verwendung kleiner Buchstaben) 
einer ungestörten Grundströmung (große Buchstaben) überlagert werden. Gemäß diesem Nähe- 
| rungsprinzip ist der Stromlinienwinkel der Gesamtströmung gleich © + d, während für den RER 
SEE Geschwindigkeitsbetrag V -+ v gesetzt wird; außerdem bleibt der Störungseinfluß auf dieMach- Br. 
schen Wellen unberücksichtigt, so daß die ‚‚charakteristischen Koordinaten‘ der ungestörten Be 
- Überschallströmung für die Gesamtströmung unverändert beibehalten werden können. Dieses _ 

-  Überlagerungsprinzip gibt dann die Möglichkeit, auch den Einfluß einer kontinuierlichen An- 
ordnung von schwachen Störungen auf eine Strömung zu studieren, indem diese Anordnung | 
aufgefaßt wird als Grenzwert der Summe von diskret angeordneten Störungen. SER 
Die Durchführung dieses Gedankens wurde zunächst an dem ganz einfachen Beispiel der Er 
ebenen Zirkulationsströmung um einen a vom Radius / gezeigt, deren Geschwindig- ec 
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Oak +75 ” Tautet (r, » ebene Polarkoordinaten ; a* kritische Schall- 


geschwindigkeit, ka* Geschwindigkeit el dem Zylindermantel). Die Schallgrenze liegt hier 

" unabhängig von Randbedingungen von vornherein fest. Durch geeignete Wahl des ‚‚Eigenwertes“ 
k läßt sich erreichen, daß eine vom Zylindermantel ausgehende Machsche Welle nach — unter = 
Umständen mehrfacher — Reflexion an Schallgrenze und Zylinderwand schließlich an ihren 
Ausgangspunkt zurückkehrt. Durch diese ‚Periodizitätseigenschaft‘“ zusammen mit dem Toll- 
mienschen Satz über die Konstanz der Sprunggröße A wird der Verlauf einer Singularität ein- 


3) @. I.Taylor, Strömung umeinen Körperin einer kompressiblen Flüssigkeit. Z.angew. Math. Mech. Bd. 10 (1930), 
5. 334— 345. 

2) G, Guderle-y, Störungenin ebenen und rotationssymmetrischen Schall- und Überschallparallelstrahlen. Z. angew. 
Math. Mech. Bd. 25/27 (1947), 8. 190—195. - 

») w.Tollmien, The direct hodograph method in the theory of the flow of compressible fluids. Erscheint in einem 
Sammelband Januar 1949 im Verlage J.'w. Edwards, Ann Arbor, Michigan U.S.A. 


. erwünschte Möglichkeit, den Vergleich mit den Taylorschen Ergebnissen anzustellen. Erzeigt 


R: deutig a de man eat. ale 36 esamte, 
eine periodische kontinuierliche Singularitätenverteilung 
- nächstliegend ein cosinus-Ansatz an. Hieraus 
besonders bemerkenswert, daß durch Anwendung der direkten F 
- Störungsverlauf grundsätzlich durch bloße Quadraturen bestimmen läßt, 
"bisher benutzten Methoden (z. B. die ‚Methode der kleinen 


lassen sich en 


Schwingungen“) 
Differentialgleichung erfordern. Für die radiale en "besteht. 
eine qualitativ durchaus befriedigende, dngegen quantitativ nur USER. OT 


Ein Pohlhausen-Verfahren zur ‚Berechnung laminarer 
kompressibler Grenzschichten EN: 
Von J. Ginzel in Göttingen | u ee = 

Die orehnge in der laminaren Grenzschicht werden beherrscht N nr 
chung, Kontinuitätsgleichung und Energiegleichung. Durch Integration in y-Richtung (senk- 


recht zur Wand) von der Wand bis in die Außenströmung bekommt man den bekannten Impuls- 
integralsatz und einen ra sie Dem Impulsintegralsatz geben wir die Form 


wobei-, 
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wund U sind die Geschwindigkeitskomponenten in x-Richtung (parallel zur Wand) in der Grenz- 

schicht: bzw. in der äußeren Strömung, oe und 0, die Dichte in der Grenzschicht und in der 

Außenströmung; 9, heißt Verdrängungsdicke, 9, Impulsverlustdicke. Der Index w deutet die 

Werte an der Wand an. 
' Für den Energieintegralsatz geben wir eine Zwischenform an: 


z len wars [uw oa- le), [nk SR Re, 


i„ und i, sind die Wärmeinhalte in der Grenzschicht und in der Außenströmung, « der Zähig- 
keitskoeffizient, o die Prandtlische Zahl. Wir definieren eine Enthalpieverlustdicke 


ou D 


und eine entsprechende Länge für die Reibungsverluste 


a 


2 


dy. 
Setzt man 


so erhält man den Energieintegralsatz in der Form einer gewöhnlichen linearen Differential- 
gleichung, in deren Lösungsformel sich eine Integration durchführen läßt. Er heißt dann 
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res der durch | ntegration 
ERRE: ent ‘ ht, so erhält ‚man 23 \ 


in der c das Gleichgewicht zuischen k kisstlieher Brig, re et Wärme- BER* 
ergang v von Se a in die or besonders deutlich ie 2. der Koeltizient Sa ne 
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\ Gräben Fe und q als vPüuktisnen eines noch zu Ferien: Freenet ns der Mach- 


. Zahl M und des. Verhältnisses- BrLE aus Tomperaturgrenzschichtdicke und Geschwindigkeit. 
 grenzschichtdicke kurvenmäßig r- —n == : ev 
fürs ein bestimmtes Verhältnis aus 
_ Wandtemperatur zu Ruhetempe- 
% ratur festgelegt werden. (Nach _ 
‚der Form der Gleichungen treten ’ 
die Grenzschichtdicken selbst zu- 
Bi nächst nicht auf, sondern statt 
dessen die Impulsverlustdicke 9, 
unddie Enthalpieverlustdicke d,, 
- dieden Grenzschichtdicken öund 


 ö, proportional sind.) Darin, 20 SS 

daß jede der Gleichungen von ee : 
beiden Grenzschichtdicken ab- _ Impuintegrauatt  — —  Energieintgraliar 

. hängt, zeigt sich die Koppelung BE {ar Bild ı. Sr 

‚der Gleichungen. Der Berech- ; ae \ 


nung der Kurvenblätter (von denen Bild 1 ein n Beispiel zeigt) wurde eine Klasse von Potenz- 2 
_ profilen zugrunde gelegt, die außer von der Mach-Zahl und dem Verhältnis der Grenz- 


rd 
‚schichtdicken von einem Parameter Ara Ö? Oa/lta ( (Pohlhausen- Parameter) - abhängt. 


Diese Kurvenblätter liegen vor für zwei Wandtemperaturen (T, »100°G 4.1: 13 

und 7, © — 25°C, d.i. w/% : 
—=0,85)im Bereich 1,2<M<2, 
1<6/ö<15 und für A, 
Bereiche bis zur Ablösung. Man 
„ rechnet abwechselnd einen 
Schritt in der Lösung dereinen 
und einen Schritt in der Lö- 
sung der anderen Gleichung, 
| wobei die erste Gleichung mit 
E -z auch ö und die zweite mit 2, 
- auch ö, liefert. Das berechnete 
Verhältnis ö,/ö dient als Para- 
meterwert für den nächsten 
Schritt. Bild 2 zeigt für eine 
ebene Wand der oben im Bild 
angegebenen Form die ört- 


4 


liche Schubspannung Ki 


Be ns eine örtliche Wirmeibergangsah a gemäß der Gliederung 
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e- ? Für die Zar Wand wurde die Rechnung mit 

EN “ durchgeführt (Bild2%a und 2b) und geringer Einfluß dieses Gesetzes auf 

0 festgestellt: Bild 2c gehört zur gekühlten Wand und zeigt die praktische I $ 
Be. der Wandschubspannung von der Wandtemperatur. Die Wärmeübergangszahl hängt stark: r 
E; Verhältnis der Grenzschichtdicken ab und wird für große ö,/ö ‚ auch wenn die Bet 


temperatur noch niedriger ist als die Wandtemperatur. Dieser Fall ist im im rechten Teil des cent: 


Br, Bildes, zu erwarten, genauere Zahlenwerte RER noch nicht vor. | 

B?: Erzeugung von Zirkulationen beim Schütteln von Gefäßen 

> ? Von L. Prandil in Göttingen. - 

e Die Anregung zur Beschäftigung mit dieser Frage leitet sich von einer Anekdote her, die 


mir mein Freund Carl Runge vor reichlich 40 Jahren erzählte. Sie handelt von Laplace. 
r Auf diesen sprach eines Tages bei einem Glas Wein ein Besucher mit emphatischen Worten ein, 
2 wie herrlich es sei, daß man nun die physikalischen Vorgänge in allen Einzelheiten berechnen 
322 könne. Laplace machte jedoch eine ablehnende Handbewegung; er nahm die Flasche mit dem 
Rest des Weines in die Hand und sagte, während er die Flasche kräftig im Kreis herumschüttelte: 
„Rechnen Sie doch das!“ 

Diese sicher wahre Anekdote, für die ich ee die Verantwortung dem nun schon 
seit 20 Jahren toten Runge überlassen muß, hat mich in der Zwischenzeit noch oft verfolgt. 
Die ältere Theorie, die die Flüssigkeiten als reibungslos voraussetzte, hatte für die entstehende 
Umlaufbewegung keine Erklärung, sie gab vielmehr für jede Translatationsbewegung eines Ge- 
fäßes das Resultat, daß hier nur Potentialbewegungen möglich seien, eine Zirkulation also nicht 
entstehen könne. Auch die Einbeziehung einer kleinen Reibung, die nur Abweichungen in den 
Grenzschichten an den Gefäßwänden hervorbrachte, änderte nichts an dem Zustand der Flüssig- 
keit, abgesehen von dem geringen, von den Grenzschichten eingenommenen Raum. Die Be- 
wegung blieb auch hier eine Potentialbewegung ohne Zirkulation, es sei denn, daß ein Zusammen- 
fließen von vorher getrennten Flüssigkeitsteilen stattfand. Ein solches Zusammenfließen tritt 
freilich immer dort ein, wo eine Welle sich überschlägt. Da die Yorwärtsbewegung im Wellen- 
berg stärker ist als die Rückwärtsbewegung im Wellental, bleibt nach dem Zusammenfließen 
ein Impuls in 10 Richtung der Wellenbewegung übrig, siehe Bild 1. 


Zur Aufklärung der theoretischen Frage entschloß Ro mich vor etwa zehn Jahren, einen 
kleinen Apparat zu bauen, mit-dem dem Problem auf den Leib gerückt werden konnte. Die Ver- 
suche kamen aber des Krieges wegen nicht mehr zur Durchführung und haben bis vor kurzem 
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Bild1. Sich überschlagende Bild2. Schüttelapparat, Achse A und Berdfest, A angötrieben. Gefäße werden auf 
Wasserwelle Platte Z befestigt, r = Schüttelradius. CD bleibt bei der Bewegung parallel zu AB - 
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ser en 05 =) zeigten Ka daß a ee 
al ukelnde ‚Oberfläche 2 blieb, in en mit der ‚Theorie keine ö 


a un der Theorie Br endlichen Wellenamplituden, die Teilchen i in 'dön Wellenbergen x 
größeren Weg vorwärts machen, als sie in den ‚Wellentälern zurückschwingen) ; nach Dei 
utsaniem Abstellen war aber jede Umlaufbewegung verschwunden. Sowie aber auch nur eine 
Wellenkrone schäumte, blieb jedoch eine kleine Umlaufbewegung zurück, und bei starkem 
i Se konnte man beliebig starke Zirkulationsbewegungen erhalten. 
3m "Bei den im Vortrag vorgeführten Versuchen wurden erst mit großem Radius r (rd. 60 en 
ei erziwungene Schwingungen gezeigt und dann mit kleinem Radius (rd. 2 mm) Resonanzvorgänge. 
. Da die Frequenz der Eigenschwingungen der Wassermasse sowohl von der Wellenamplitude 
I "wie auch von einer etwa vorhandenen Umlaufbewegung des Wassers abhängt, ergaben sichaller- 
hand reizvolle Einzelheiten. ‚Bei überkritischem Betrieb z. B. bleibt die Anfachung aus. Wenn = N 
. „man aber auf eine Frequenz ein wenig über der Resonanzstelle der ruhenden Flüssigkeit einstellt 
_ und vorher für große Zirkulation in der umgekehrten Richtung sorgt, dann ergibt sich zunächst 
keine "Anfachung; wenn dann aber die Umlaufbewegung durch Reibung allmählich ermüdet, 
"kommt plötzlich eine hohe ‚Resonanz und damit on auch. Umkehr der Zirkulationsbewegung Me: 
zustande. - 
j Es mag bemerkt werden, daß sich die meist lJaminare Grenzschicht bei den Beobachtungen 4 
ni = nie irgendwie besonders bemerklich macht, abgesehen von ihrer bremsenden Wirkung auf etwaige 
Umlaufbewegungen. Beim Abstellen des Apparates zeigt sich durch die ‘hierbei auftretenden 
-  Druckfelder eine rasch wieder verschwindende Umlaufbewegung der Grenzschicht entgegen- 
ER der durch Schütteln erzeugten. 
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= Neuere Geräte für Strömungsmessungen & 
E : : | Von Alb. Betz in Göttingen : Be 
Er? Unter den Geräten für Strömungsmessungen sind die Mikromanometer von besonderer | 
E Bedeutung. Es gibt sehr viele Konstruktionen, die den besonderen Wünschen nach hoher Ge- 
nauigkeit, Verkürzung der Einstellzeit, Bequemlichkeit der Ablesung u.a. Rechnung tragen. ER 
Wegen seiner bequemen und genauen Ablesung auch aus größerer Entfernung hat sich das a 
Projektionsmanometer nach Betz gut bewährt. Dabei erscheint im Projektionsbild auf einer 
Skala jeweils eine bezifferte Marke. Die Ziffern geben den Druck in kg/m? und die Stellung der 
Marke auf der Skala die Bruchteile an. Der Meßbereich ist 400 kg/m?, die Ablesegenauigkeit 
nahezu 10-?kg/m2. Durch eine konstruktive Maßnahme konnte die gleiche optische Einrichtung 
auch für die 10fachen Drücke verwendbar gemacht werden. 
Ein anderes baulich besonders einfaches Manometer gestattet hohe Empfindlichkeit bei 
kleinen Drücken mit einem ausreichenden Gesamtmeßbereich (250 kg/m?) zu verbinden. Bei ihm 
> __wird der Meniskus in einer Kapillare durch Schwenken des Meßschenkels auf eine festgegebene 
Nullstellung eingestellt und die Druckhöhe aus dem SchwenkWinkel bestimmt. 

Für extrem kleine Druckdifferenzen bis herab zu etwa 10-°kg/m? ist ein Gerät von 
Reichardt geeignet, bei dem der Druck auf einen kleinen, an einem Torsionsfaden hängenden 
Drehkolben wirkt und ihn verdreht. 

Es besteht noch Bedarf nach einem bequemen Gerät ähnlich dem Projektionsmanometer 
für Drücke bis herab zu etwa 10-®kg/m?. Wichtiger wäre aber die allgemeinere Aufgabe, einen 
Druckoszillographen zu schaffen, der auf sehr kleine Drücke in sehr kurzer Zeit in registrier- 
fähiger Form anspricht. Die Arbeiten hieran sind leider durch die Zeitumstände immer wieder 
unterbrochen worden, sollen aber jetzt wieder aufgenommen werden. 

Vielfachmanometer halten gleichzeitig eine größere Anzahl von Drücken meist photo- 
graphisch fest. Um die große Auswertungsarbeit dieser Photogramme zu umgehen, wurde eine 


weiter verarbeitet wurden. Zur Registrierung rasch rander! 
war von Drescher ein rasch rogistrierendes ge 
Erfolg verwandt worden. EIN, NS CK 


Bei der Geschwindigkeitsmossung rfüllte Wi 
Messung kleiner Geschwindigkeiten, Das übliche ist, 
geeignet, für die Praxis aber meist zu unbequem. Für die Messung kleiner W 
keiten bis herab zu etwa d mm/s hat Ben mit gutem m rt Su armani 5 
Gerät entwickelt, das die Mängel Hitzdrahtes vermeidet. 
Für Richtungsmessungen wurde von Conrad ein ne chin Gert en 
"mit dem man auch in die Nähe von Wänden oder störenden Körpern 
Zum Schluß sei auf ein von Ludwieg entwickeltes Gerät zur a ar Mesung der Wi 
"spannung hingewiesen, das auf der Wärmeabgabe einer geheizten Fläche e % 
- ist bereits sehr erfolgreich für wichtige Forschungsarbeiten eingesetzt ver Vortrag. 
von u Ludwieg und W. Tillmann). nr 


4 Zur Theorie der Gefügebildung beim Erstarren reiner eaie) 
SR Von Hans Fromm in Alzey > Br Ba 


Bei Erstarrung unter Wärmeentzug entsteht ein kristallines Gefüge durch Diese von % 
- Kristallkeimen und deren Wachstum. 


= = ß 6; ° Phasenaustausch zwischen Kristall und Schmelze ergibt bei Überwiegen der Anlagerungen L 
Bi - „über Ablösungen ein Wachsen mit der mittleren Geschwindigkeit c (@), definiert mar, Sa 
Be ae Peer > REBEL |. 


(Vz = Volumen, F = Oberfläche eines Kristalls zur Zeit; © = absolute Temperatur). Nahe _ 
Br dem Schmelzpunkte (©;) lautet das „Wachstumsgesetz‘"etwa 


3 2 c(0) zo (0 — 9) mi SS <0.. num ann (2). 
2% Mit starker Unterkühlung (@—0) geht auch e—0. . 
2 Auch die „Keimung“ k (= Anzahl der in Einheiten des Schmelzevolumens und der Zeit 


gebildeten Keime) erfolgt in Wettbewerb von Bildung und Auflösung. Die Auflösung steigt 
mit der „Keimdichte“ » (= Keimzahl je Volumeneinheit). Bei n = 0 sei die Keimung &,- 
Mit v = Volumen eines Keimes und &y = Auflösungsfaktor sei gesetzt: 
„Keimungsgesetz‘“: | 
k=k,(9) 1 en)— (9) on hl) NOn.... 2... (3), a 
worin der „Keimdichtefaktor‘‘ 
N=[11—(1+&(9)k, (ON en). 


Keimungsgleichgewicht (k = 0, also N = 0) fordert zu jeder Temperatur eine Gleichgewichts- 
keimdichte n, (9). 


7 Die Zeit zwischen zwei Auflösungen ist ?, = U(ka - vn). Bei Gleichgewicht muß jeder Keim 
gerade n solcher Zeiten überleben. Also ist Ai „mittlere Lebensdauer“ eines Keimess 
u= nt, = 1/(kı t) a Be ra RENTE (4). 


In dem Ansatz 


v v x 


ko kget®.e® und u=ket.e® 


berücksichtigt e =. den Wärmetod bei tiefen Temperaturen. Anpassung an Versuche mit Eisen 
liefert am Schmelzpunkt etwa n, ='10tem-3; 2, =5000 Jahre. Überhitzung um 20° (50°) C 
senkt die Werte auf n, = 0,05 cm? (£4 = 10-*sec). 

Unterhalb © verschwinden Keime auch durch Wachstum. Die Zeit zum Wachsen um 
eın Atom heiße „mittlere Bestandzeit‘. Diese Zeit wird !) schon bei geringer Unterkühlung 


* Eine ausführliche Darstellung ist in der Zeitschrift „Optik“ 3. Ba. (1948), S.137—164 gegeben. 
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re Teegalglsichung für. a( Er stellt Er Bene setz“ a Sie Jäßt sich durch h wieerhoflen FE 


. Ableitung nach der Zeit t in eine Differentialgleichung. umwandeln, 


” Die Temperaturfunktionen kl ©) und c (9) folgen der Änderung von On mit den Zeiten ce 
ar bzw. t gemäß einer „Wärmebilanz‘: Bedeuten o, = Keimungswärme je Keim; o,, gg, 48 


‚so ‚lautet die Bilanz: he 
3 RER TEN OH a1 2 9 ; 
=. ee er + Metal: Ted 


ne Er ee ‘ 

5 x E Ey = 2 S > v 
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sind also auch als Funktionen des Ortsvektors 1(2, Y,2) Aa, und es gilt die Differential- 
as der „Wär meleitung“ (= Leitfähigkeit); 


> . [8290 820ı 20 
3% 11 | . (5,4). 


a = -4.49= EM TRrT 


bedingungen über die Wärmeabgabe bestimmt. 


Für Erstarrung bei konstanter Temperatur (0 < 6s), die nur Ba möglich ist, 
ergibt sich eine geschlossene Lösung: 


on = Co (pt) . cos (pt) mit .p=Y2nok 


T=aj/(2p); Kr = 0,79 (kje)?* De ee A 
: „Korngröße. “.Radiusr=0-(T—t1); Re rmx =eT= 0,99: e/k 
F „Größenverteilungsgesetz“: dKr/dr = — (kje) - «a (t) 


Br: Wegen der Durchführung bei gleichmäßiger Wärmeabfuhr und ihrer Ergebnisse muß auf 
DER die ausführliche Darstellung verwiesen werden !). 

es Die Korndichte K7 und die Größenverteilung beziehen sich auf die Volumeneinheit. Die 

"Beobachtung an ebenen Schliffen zeigt Schnitte durch die Körner, deren Größenverteilungs- 

. gesetz wesentlich abweicht. Während z.B. nach (6) das größte Korn (R) das häufigste ist, 

-ist es im Flächenbild. das Korn rz = 0,6: R. Zwischen „Flächendichte‘“ Zz und „Raum- 

dichte“ Kp besteht eine Beziehung Kr = 'Z7?", in der p von der Korngestalt abhängt; 

für die: hier angenommene kugelnahe. Gestalt ist g = 0,9. 
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se - Kristallisationswärme und spezifische Wärmen von Kristall bzw. Schmelze, alle auf die 4 
Sa 'Tumeneinheit bezogen, und. = Bea re = Wärmeabfuhr je Volumen- und Zeiteinheit, _ 


BR nen kann die Dushwarme nur durch Leitung abgeführt werden. Alle RR 


Damit ist Äas Problem formuliert und wird durch Anfangsbedingungen n, (rt); ©,(r) und Rand- 
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ist z. B. sicher nicht mehr streng anwendbar auf Strömungsvorgänge in nen um & 


treten in Gasmischungen mit inhomogener Druck- und Temperatu 
 gänge auf, die eng mit den Erscheinungen der sog, Drück- und hermoditinsion 


für Gase, leitet man exakte Bewegungsgleichungen ab für gewisse Mittelwerte, die die Rolle 
. der'hydrodynamischen Variablen übernehmen, außerdem,aber noch den Reibungsspannungen 


Der . Begriff der statistisch 
Von Max in 1 Horb 


| a, 
Du en Be Br = URS 
Kontinuumstheorie, ihre durch die molekulare Struktur der Materie I 


deren Dimensionen nicht sehr groß sind gegenüber der mittleren freien Khan Be, ER. 
7 ei 


hängen. Man kann sich die Frage vorlegen, ob es eine verallgemeinerte Sn bt 33 


mit wenigen, nur von den Koordinaten und der Zeit abhängigen Feldgrößen, in der diese Er 
scheinungen miterfaßbar sind. Diese Fragestellung ist z. B. von technischer Bedeutung für die 
Aerodynamik der hohen Atmosphäre, der sog. „„Superaerodynamik“‘. Die EN 
allgemeinerten Strömungslehre bildet die kinetische Theorie der Gase und "Flüssigkeiten, und 

ist daher eine statistische Disziplin. Sie kann „statistische Hydrodynamik“ genannt ee 
in Analogie zu der Schrödingerschen Bezeichnungsweise „statistische Thermodynamik“ für 
die Molekulartheorie der Gleichgewichtsvorgänge. Ausgehend von der Boltzmanngleichung 
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und den Komponenten der Wärmestromdichte entsprechen. Für die sog. Enskogschen Nor- 
mallösungen ’der- Boltzmanngleichung sind die letzteren auf die hydrodynamischen Größen 
und deren räumliche Differentialquotienten bis zu beliebig hoher. Ordnung zurückführbar. Für 
die Enskogschen Normallösungen bildet die Lösung der Bewegungsgleichungen für die hydro- 
dynamischen Variablen ein vom kinetischen Gesamtproblem abseparierbares Teilproblem. Der 
Gaskinetik bleibt das Problem der Bestimmung der Materialkonstanten in den Beziehungen 
für die Reibungsspannungen und die Wärmestromdichte aus speziellen Molekülmodellen. Die 
Enskogschen Normallösungen bilden nicht die allgemeinsten Lösungen der Boltzmann- 
gleichung. 

Ein schwieriges.Problem bildet die Frage der Randbedingungen an Grenzflächen. Die 
üblichen Randbedingungen der Hydrodynamik sind nicht ausreichend, da die Ordnung des 
‘Systems der hydrodynamischen Bewegungsgleichungen durch das Auftreten von höheren räum- 
lichen Differentialquotienten der hydrodynamischen Variablen größer aa ist. Dieses 
Problem ist z. Zt. noch nicht befriedigend gelöst, 
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Neue Rechnungen zur statistischen isotropen Turbulenz 
2 Von J. Rotta.in Göttingen 


In zwei Arbeiten, die demnächst im Druck erscheinen werden !), haben €. F. v. Weiz- 
säcker und W. Heisenberg eine Theorie der statistischen Turbulenz entwickelt. Ich habe 
versucht, diese Theorie etwas zu erweitern mit dem Ziel, auch für die größeren Turbulenzelemente 
eine brauchbare Annäherung an die Wirklichkeit zu bekommen. 

Das Spektrum der Turbulenz einer inkompressiblen Flüssigkeit wird durch eine Funktion 


an 


F (k) der Wellenzahl k = Wellenlänge - definiert, so daß die gesamte Turbulenzenergie je BT 


einheit 


beträgt. u® ist der quadratische Mittelwert des Vektors der turbulenten Bewegung. 

Den physikalischen Vorgang beim Abklingen der Turbulenz, bei dem die Turbulenzenergie 
durch turbulenten Austausch von den größeren Elementen in die kleineren Turbulenzelemente 
geleitet wird und in den kleinsten Elementen schließlich durch molekulare Reibung in Wärme 
umgesetzt wird, kann man nach Heisenberg, wenigstens näherungsweise, rechnerisch erfassen, 
Zur Bestimmung der spektralen Verteilung F'(k) fehlt aber noch eine Angabe, die die zeitlichen 
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“anftaßt, 2 ‚den ang (2) al, ie 6) läßt sich der en er Zwischen 
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Mit Hilfe der Gl. (2), (3) und (5) erhält man schließ- 
lich als Bedingung für das statistische Gleichgewicht 
der Turbulenz AR 


eine von k unabhängige 


2% 
t 


=> = Dabei == 

2 abei ist C E 

- Funktion der Zeit. Für hinreichend große Wellen- 
zahlen % nähert‘ sich E, dem konstanten Endwert 


Be von E nach Gl.(1), und k F(k) geht gegen Null. Bild 1 
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Damit geht Gl. (8) über in 


lim’. S,= konst . * 
k> oo = 


Für die großen k-Werte ist Gl. (8) also identisch mit der Bedingung, die den Untersuchungen 


von v. Weizsäcker und Heisenberg zugrunde liegt. x 
ET AT 

a de N h. v.Kärmän, V.Intern. Congr. Applied Mechanics (1938), 8.347. 

») L.Prandtl, Nachr.d. Akad.d. Wissensch. Göttingen (1945), 8. 6. 
R 4) Pür den Speziatfall & = 0,5 hat W.Heisenbergineiner weiterenim Druck befindlichen Arbeit auf ähnliche Weise 
die gleiche Beziehung gefunden. Die vorliegenden Überlegungen hat der Verfasser durchgeführt, ohne von der Existenz dieser 
Arbeit von Heisenb erg zu wissen. 
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Die Gl. (8) ergibt mit der von W. Heisenberg angegebenen ‚orme für 8, durch d en & 
die turbulenten Austauschvorgänge berücksichtigt werden, eine Integralgleichung zwischen den Bi 
Größen F(k) und k. Der Exponent & des Dissipationsgesetzes liegt nach bisherigen 
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zwischen 0,5 und 1. Die neuesten Messungen von G. K. Batchelor und A. A. Townsend 
die als die zuverlässigsten anzusehen sind, ergaben a = 0,5. Für diesen Fall & = 0,5 wurde 
die Gl. (8) schrittweise integriert. „en DER 
Die vorliegenden Versuche von L. F. G. Simmons’und C. Salter°) über die Ausmessung 
der spektralen Energieverteilung bei Windkanalturbulenz bieten die Möglichkeit, die Theorie 
‚mit der Erfahrung zu vergleichen. Die Rechnungen sind in sehr gute Übereinstimmung mit den 

- Messungen zu bringen, Bild 1, so daß einige Wahrscheinlichkeit vorhanden ist, daß die Theorie 

den physikalischen Gegebenheiten gerecht wird und in genügend großem Abstand vom turbu- 

lenzerzeugenden Gitter tatsächlich ein Gleichgewichtszustand in dem definierten Sinne besteht. 


») G.K.Batchglor, A.A.Townsend, Proc.Roy. Soc. A Vol.190 (194%), 8.534. 
e LEG et, ns, 0.Salter, Proc. Roy. Soc. A Vol.165 (1938), 8.73, 


Dreidimensionale Strömung um schlanke, nicht rotationssymmetrische 
Körper mit gerader Achse 


Von Fr. Riegels in Göttingen 


j Für den dreidimensionalen Fall der reibungslosen inkompressiblen Strömungen sind nicht 
so viele Lösungen der Differentialgleichungen bekannt wie im zweidimensionalen, so daß man 
aus bekannten Einzellösungen nur wenige Anhaltspunkte für die Strömung um Körper allge- 
meinerer Gestalt besitzt. Daneben birgt aber auch.die quantitative Behandlung solcher Strö- 
mungen noch erhebliche praktische Schwierigkeiten. Bekannt geworden sind bisher nur einige 
wenige Verfahren, die gestatten, wenigstens die Strömung um Rotationskörper zu berechnen. 
Ich nenne hier nur die Namen Rankine, Fuhrmann, v. Kärmän und Lotz. Die folgende 
Untersuchung setzt sich daher zum Ziel, quantitative Angaben für einen allgemeineren Fall zu 
gewinnen, indem sie die Strömung um schlanke, aber nicht mehr notwendig rotationssymmetrische 
Körper näherungsweise behandelt. 

Zu diesem Zweck ersetzt sie den beliebig vorgelegten Körper, von dem nur vorausgesetzt 
sei, daß seine räumliche Erstreckung in einer Koordinatenrichtung, etwa der z-Richtung, wesent- 
lich größer als in den beiden anderen dazu senkrechten Richtungen (die annähernd von derselben 
Größenordnung sein mögen) sei, zunächst dusch einen mittleren Rotationskörper, indem jeder 
Querschnitt des vorgegebenen Körpers durch einen Kreisquerschnitt mit gleicher Fläche an- 
genähert wird. Die Strömung um diesen Ersatzkörper kann mit Hilfe von Quellringen bestimmt 
werden und dient als erste Näherung. Weitere Näherungen werden erhalten, wenn man diesen 
Ersatzkörper mit weiteren Quellringen, jedoch periodisch veränderlicher Intensität belegt. 
Setzt man kleine Abweichungen von der rotationssymmetrischen Form ‚voraus -— eine Forde- 
rung; die in den meisten Fällen ohne weiteres erfüllt ist —, so lassen sich die an sich recht kom- H 
plizierten Rechnungen weitgehend linearisieren und man kann den Verlauf der Intensitäten i 
dieser Belegungen über & so bestimmen, daß bei der Überlagerung all dieser Quellströmungen 
über eine Parallelströmung schließlich die vorgegebene Körperform resultiert. 

Um die angedeuteten Rechnungen durchführen zu können, ist es notwendig, die Ge- ' 
schwindigkeitskomponenten von Quellringen mit periodisch über den Umfang veränderlicher 
Intensität zu kennen. Aus diesem Grunde wurde die dabei immer wieder auftretende Funktion 


u te a 


nf? 
cos2np 
@ (k? = a 1 1: II 5 ber a yasar ea 1 
n() Irre 6 
ß | 
für die Ordnungszahlen n=0 bis n=7 der ersten sieben Fourierglieder tabellarisiert. Die Er- | 


gebnisse der bisherigen noch nicht ganz abgeschlossenen numerischen Untersuchungen — über 
die später noch zu berichten sein wird — mögen dahin zusammengefaßt werden, daß sich an- 

. scheinend eine beiriedigende Näherung ergibt, wenn man zu der Tangentialgeschwindägkeit der 
Anströmung, die also z. B. in einem Schnitt & = konst über den Umfang veränderlich ist, noch 
die von der Verdrängung des ersetzenden Rotationskörpers herrührende, über den Umfang un- 
veränderliche Tangentialgeschwindigkeit hinzufügt. 


ma. = 7: . 


er Kenntäiise, de Verhaltens der Wendehubeh ren in 
a , in wir a, Biene: ee 


174 


chicht Grande ist, für die die u Bezhung 


= an iechubesanndng u = Geschwindigkeit i in der Helkingischiehe, A Were 
SE zuee) gilt, hängt die übertragene Schubspannung außer von der Stoffkonstanten | a 


du 
nur von dem Geschwindigkeitsgradienten ne an der Wand ab. Prinzipiell ‚könnte man also 


eh Messen des Geschwindigkeitsgradienten die Waschubsnnng bestimmen. Mit den.. 


x üblichen Hilfsmitteln, wie Pitotrohr oder Hitzdraht, ist dies bei der geringen Dicke der Unter- ie 
. schicht im allgemeinen versuchstechnisch nicht durchführbar. ‘Man kann aber den Geschwin- 


Kiaherkiärabenten mittels einer  Wärmeübergangsmessung erfassen, "Die Wärmeabgabe eines 


kleinen, a und wärmeisoliert in der Meßwand eingebauten, e Metallklötzchens ist 
3 


nämlich 2 © bzw. Tu proportional, wenn. die. Tiefe des Klötzchens in Strömungsrichtung 


. so klein gewählt Br daß die sich, bildende Wärmegrenzschicht innerhalb der laminaren 
* Unterschicht verbleibt. Ein auf dieser Grundlage beruhendes ‚Meßgerät wurde entwickelt. Aus 
mechanischen Gründen konnte die Tiefe des Klötzchens doch nicht so klein ausgeführt wer- 
‘ den, daß die "Wärmegrenzschicht vollkommen in’ der laminaren Unterschicht bleibt. Wie noch 
gezeigt wird, wird- dadurch der ‚eindeutige Zusammenhang von Wärmeahgabe und Schub- 
spannung nicht beeinträchtigt. 

Aus früheren Reibungsschichtmessungen ad theoretischen Betrachtungen wre ge-. 
schlossen, daß auch bei allgemeinen turbulenten Reibungsschichten mit Druckgradienten in 
- Wandnähe, d.h. in der laminaren Unterschicht, im Übergangsgebiet und im wandnahen Teil 
des voll turbulenten Bereiches, das gleiche universelle Geschwindigkeitsgesetz 


u .,(ury er 
7 el) ee 


(ut = V*- Schubspannungsgeschwindigkeit, 0 = Dichte, » = z- = kinematische Zähigkeit) 


wie bei der Plattenströmung (mit konstantem Druck) gilt. Im Bereich der laminaren Unter- 
schicht ist dieses Gesetz natürlich identisch mit Gl. (1). Unter der Voraussetzung der allge- 


r et F 5 T, 
meinen Gültigkeit der Gl. (2) wurde für den örtlichen Widerstandsbeiwert c; (= —--; U= Ge- 


% 
| — U: 
schwindigkeit außerhalb der Reibungsschicht‘ folgende Formel abgeleitet: e 
x 2 
er Re Sl re N ET REST ON 
2 A / ( ” 
Ran Bier 
er (re Y Nr 
Polar . 


Dabei ist re die mit der Impulsverlustdicke ö, er Reynoldssche Zahl’und y ein Profil- 


parameter, ‚der angenähert gleich dem Wert von 7 an der Stelle y= ö,ist. Bei der Platten- 


3) Enthencher Bericht über die folgenden beiden, demnächst im Ingenieur-Archiv erscheinenden Arbeiten:- 
dwieg, Ein Gerät zur Messung der Wandschubspannung turbulenter Reibungsschichten. — H. Ludwieg und 


H. Lu 
"Ww. a Imann 5 unsersuchingen über die Wandschubspannung in turbulenten Reibungsschichten. 


en ist dieser Profilparameter, mit yo re ı von on re 2 
_ jenige Funktion, die c; mit re verknüpft. Diese, beiden Funktionen Yg =, ol be 
für die Plattenströmung sind aus experimentellen Untersuchungen hinreichend } 
liefert obige Kettenfunktion Gl. (3), die übrigens gut konvergiert, den 
- für allgemeine Reibungsschichten als Funktion des Profilparameters y und der m 
Zahl re. Das Gesetz Gr. (3) wurde durch mehrere Versuchsreihen, die in n URE 
verzögerten Strömungen durchgeführt und bei denen die Wandschubspannungen mit | des 
- neuen Meßgerätes. unmittelbar gemessen wurden, gut bestätigt, Dieselben Messungen erfüllten 
natürlich auch die bei der Ableitung vorausgesetzte Gl. (2). Daraus folgt, daß das 7 
dann die Wandschubspannung eindeutig ie. wenn die a ea aus der a Bet 
naren Unterschicht herausragt. 
Für den praktischen Gebrauch wurde die abgeleitete Gl. ) durch folgende einfache 


| Potenzformel R Et. 
2 4 = 0,0580 - „17% - re=0208 u SE FE el N TE, Ar (4) x 3 de 


a chkrh Wegen der angenäherten Einparametrigkeit der turbulenten Geschwindigkeits- 
profile kann man an Stelle des Profilparameters 7 den bei Reibungsschichtuntersuchungen häufig 


ke = Verdrängungsdicke) einführen und erhält dann 


verwandten Profilparameter Au (=, 


aus Gl. (4) folgende Formel: { . 
| %,=0;246 - 10-0 Bas gg 208...) Du Se SE BE (;) 


Aus en Formeln für c; folgt in Übereinstimmung mit den ER} daß bei Rei- 
bungsschichten mit Druckabfall bzw. -anstieg größere bzw. kleinere c,;-Werte als bei der Platten- 
strömung gleicher re-Zahl auftreten. Es lassen sich also bei größeren Reynoldsschen Zahlen 
kleinere örtliche Widerstandsbeiwerte c; nicht nur durch Laminarerhaltung der Grenzschicht, 
sondern auch durch geeigneten Druckverlauf in turbulenten Reibungsschichten erreichen. Das 
von verschiedenen Autoren festgestellte Ansteigen des Reibungsbeiwertes auf ein Mehrfaches 
des Beiwertes der Plattenströmung in Reibungsschichten mit Druckanstieg ist durch diese 
Untersuchungen widerlegt. 


Zur Frage der. Schubspannungsmessung in turbulenter Strömung 
Von H. Reichardt in Göttingen 


In der Türbulnziorecking besteht das Bedürfnis, die durch die Schwankungsbewegungen 
hervorgerufenen Schubkräfte an jedem beliebigen Ort einer Strömung unmittelbar feststellen 
. zu können. Wir haben daher im Göttinger Strömungsinstitut die Entwicklung von Meßgeräten 
für die Schubspannung, insbesondere auch für die Wandschubspannung, seit etlichen Jahren 
betrieben. 


® 


_ Über das Problem des Messens der Wandschubspannung hat bereits Herr Ludwieg vor- 
getragen, Es wäre dazu noch manches zu sagen. Ich will mich aber aus Zeitgründen darauf be- 
schränken, die Frage des Messens im Innern eines turbulenten Stromes zu behandeln. 

Die Darstellung des turbulenten Schwankungsvektors auf dem Schirm eines Braunschen 
Rohres mit Hilfe einer aus drei Hitzdrähten bestehenden Sonde des Verfassers gehört seit 1930 
zu den Demonstrationsversuchen des Institutes. Mit dieser Hitzdrahtsonde ist es später auch 
gelungen, die turbulente Scheinreibung quantitativ zu messen!). Die Kompliziertheit dieser 
Methode erschwerte jedoch ihre praktische Anwendung. Es ist daher die Frage, ob die turbu- 
lente Schubspannung nicht in wesentlich einfacherer Weise unmittelbar mechanisch gemessen 
werden kann. 
ae. Ne ren. u die Geschwindigkeit i in der Hauptströmungsriehtung und ® die 

= oo + WE a a ee (). 


der zeitliche Mittelwert des gesamten in die v-Richtung überführten «- «-Impulsstromes je Flächeneinheit. 
euv bezieht sich auf den %-Impuls der stationären Querbewegung ?, und 0 w vv BeuE die durch die Schwan- 
kungen u und v’ hervorgerufene Schubspannung. 


| 
| 
| 


") Siehe H.Reichardt, Z.angew. Math. Mech. 18 (1938), 8. 358, Naturwiss. 26 (1938), 8.404, 


| | ges agu 
ei ät, d auf die Größe. oewv en Für ou 
9 ofuder als N ee Winkel zwischen Br 


Ben Gl ist: en: für Take Ale auch für. Sinkeleurpfiedliche Denen ee 
er im Flügel. bedeutet N die Normalkraft, bei Winkelsonden die Eee er Kan ® 
Adaher mit solchen Geräten grundsätzlich den ‚Impulstransport. messen. ° 
$ . Nehmen wir als Beispiel eines Flügels eine Windfahne (ebene Platte), die wegen a 
ehaklen Vorderkante auch ‚bei kleinen Re-Zahlen auf kleinste Winkeländerungen reagiert. Die 
ER Abmessungen -der Windfahne seien klein im Vergleich zu den Abmessungen der Wirbelballen. g> 
Wird ein solches Fähnchen in der u-Richtung. festgehalten, so entsteht eine mittlere Normal- 
‚kraft, die gemäß (1) und (2) proportional oe ad owv ist. In dem Sonderfall einer ausgebil- ze 
Ren . deten Kanalströmung, in der ® verschwindet, ist N now wo. Zur Zeit laufen Versuche, eine 
'solche Schubspannungsanzeige mit. ‚einem Plättchen zu verwirklichen, das an einem Torsions- ER 
 faden hängt und dessen Drehbewegung stark gedämpft ist. Z & 
Den größten Ausschlag %, der Fahne erhält man offenbar, wenn die Richtkraft ds Fädens fr 

zu vernachlässigen ist. Bei diesem Anstellwinkel ist die von der Hauptbewegung @ hervor- Hr 
bernlene: Normalkraft m ae im Gleichgewicht. N er mittleren Normalkraft der Schwan- 
kungen wowv. Somit ist. 


a, 
% 


| | er m = a 3). 
3 < er Spa die turbulenten Schubkräfte zwei bis drei Zehnerpotenzen kleiner sind als die Stau- 2 
-  drucke, sind die auftretenden Drehwinkel & < 0%, sehr klein (<=0,3°). Durch Verlegen der 


 ,.  Torsionsachse auf die Leeseite der Fahne kann 
| man jedoch den Ausschlag wesentlich vergrößern. 
Allerdings werden dabei auch unerwünschte Stör- 
=... effekte mit vergrößert. 
£ Die experimentellen Schwierigkeiten sind 
eben grundsätzlicher Art, unabhängig von der 
gewählten Meßmethode. Stets handelt es sich 
darum, die von der 'Scheinreibung herrührenden 
sehr kleinen Effekte von den großen Wirkungen 
der Hauptbewegung zu trennen und fernzuhalten. _ 
- Bei der Ermittlung der Schubspannung aus 
der Druckdifferenz einer Winkelsonde kommt 
es vor allem darauf an, daß die ‚‚wirksamen 
Achsen“ der Sonderöffnungen durch die gleiche Kirn. nuninke nel Aue er Drurkalhinene onen i 
“Tsotäche gehen. Andernfalls wird eine zusätzliche " ;Winkelsonae, Twantschnospannnne sa TA 74.107° aus BR. 
Druckdifferenz mitgemessen, die vom Unterschied er 
des Staudruckes an den Sondenöffnungen herrührt (dieser Druckunterschied kann erheblich $% 
B größer sein als der zu messende kleine Schubspannungseffekt). 5 
Bi Eine wesentliche Aufgabe besteht darin, die Lage der ‚wirksamen Achsen‘ zu ermitteln, 
‘die mit den geometrischen Achsen der Sonden im allgemeinen nicht übereinstimmen. Sämtliche 
von uns erprobte Sondenformen zeigten in zäher Strömung eine Abweichung der wirksamen 
"Achsen von den geometrischen Achsen und zwar in dem Sinne, daß eine kräftige Querbewegung 
in Richtung auf die Wand vorgetäuscht wird. 
Dieser noch ungeklärte Einfluß der Zähigkeit auf die Druckbildung der Sonden scheint ' 
in turbulenter Strömung nicht vorhanden zu sein. Wir haben jedenfalls in einer turbulenten 
Kanalströmung mit einer Winkelsonde, die nach einem besonderen Verfahren nach den geome- 
trischen Achsen ausgerichtet war 2), Druckdifferenzen gemessen, die mit den aus dem Druck- 


\ 


2) Beschreibung erfolgt an anderer Stelle. 


'sehen kann. 


 verlässigen Meßverfahren ausgebildet werden kann, wenn die ‚von ‚höherer RR 


TREND FREE ARE 


hängt. Die vorletzte Gleichung wird nach £, durch Iteration aufgelöst. 


abfall een Schubspannungsworten annähernd. 


‚Das vorliegende günstige Ergebnis darf jedoch. a Heaber in 
Methode, mit Winkelsonden die turbulente Schubspannung zu messen, erst er 


Effekte an und in den Sonden hinreichend erforscht sind. 


B. Elastizität ad Festigkeit: 


Die Biegung der. Rechteckplatte ohne die Bernoullischen. 
oder andere Annahmen ne % 

Von Konrad Ludwig in Hannover N = SE 
Die unverformte Mittelebene wird als {z, y}-Ebene gewählt. Die Versehen er Ba 


nenten eines beliebigen Plattenpunktes werden mit &, n und & bezeichnet und nach ERBE 
des Abstandes .z entwickelt. 


| 3 Rs 3 
= 3 & 2 n=3 m [= 4: 
k=0 k=0 k=0 
Aus diesen Potenzreihen werden die für die Verzerrungs- und Spannungskomponenten berechnet. 
Die Dichte wird mit e bezeichnet. In der Bewegungsgleichung 


° iu Tun. 02 E 0) _ 80, OTyy 02 ®n 


in Nr N ne Fr 


_ 


SET EITEeE 


u a ee ee ee ne en tan re ne 


werden die Koeffizienten der Blue für k= O und 1 verglichen. Diese Gleichungen werden 
nach &,, nz und £, für k = 2 und 3 aufgelöst. 
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Die Plattendicke wird mit A bezeichnet. Aus den Randbedingungen 


RER herr 
ENT [ 0 


ee ee he) ee ee 


folgt 
BR a m. 1» a, FAR: 
14n{} 2 Er Haltz Uta re ErrT, 
| a nn 9°N0 p 
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Vorausgesetzt wird, daß’die Dicke konstant ist und der äußere Druck von der Zeit nicht ab- 
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Aus den Randbedingungen 


(Tee}, LE: {tzu} 


Pr er 1ER 370 a et h2 a ie 03%, 
ESEL FRE 2 --. 0% E € Bu: 02 B- 8 Be L a Ex 
Ma N Bez ee en GE u SO ne 
ER EEE EBEN oy 7, =‘ og E 8° oy ei 2 
wird mit Hilte der Ditferentialgleichungen für & und 29 vereinfacht. ER 
en Re) nr Puh ein es 

Yy Baus 4 B. .>z Aik= ar Ri 2 
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$ ae N Erde: 
Die Bewegungsgleichung wird über © von — 5 bis nn = integriert. Bi ee 
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„Träger. mit wandernder Last“ (strenge Lösandy: 


ee; Fe Von E. Pestel in Hannover > BL 


5 Basar Pur die Berechnung des Lastweges unter einer über einen Träger fahrenden Last unter 
= gleichzeitiger Berücksichtigung der Lastmasse und der Trägermasse wurden zwei miteinander 
verwandte Verfahren entwickelt. Mi ihr 
: Beide Verfahren beruhen auf der Anschauung und behandeln das Problem in einer von a 
den bisher angewandten Lösungsmethoden abweichenden Weise diskontinuierlich, d.h. bei der 3 
Lösung wird die Last schrittweise auf ihrem Wege über den Träger verfolgt und die dabei auf- RR: 
tretenden physikalischen Vorgänge mit elementaren mathematischen Hilfsmitteln beschrieben. Bes. 
Bei mehreren. Lasten kann nur die Masse der Hauptlast berücksichtigt werden, während die 

Anzahl der wandernden Kräfte beliebig sein kann; ebenso dürfen letztere mit der Zeit verän- 

derlich 'sein, wie z.B. die vertikalen Fliehkräfte bei Lokomotiven. Furner kann sich die Ge- 

' schwindigkeit der Last und der Kräfte beim Überfahren des Trägers ändern, wie auch die Ge- 
 schwindigkeit der Last und der Kräfte untereinander nicht die gleiche zu sein braucht. 
Die Verfahren sind auf jedes Tragwerk anwendbar, für das bei beliebiger Laststellung 

die sin: und Eigenschwingungsform angegeben werden können. 
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- formungen durch Ausbohren und Messung der Stababmessungen bestimmt werde 
Bend soll durch eine theoretische Überlegung die Gesetzmäßigkeit ermittelt werden, die z 
‘den verhältnismäßigen Längenänderungen &, & und & in Zorgal radialer und t 


Die kigenspannungen | Kin rfläch ückte = St 
von Kreisquerschnitt" > Be we 

Von 0. Föppl in Braunschweig, Wöhler-Instibut ee 

1: Es sollen die durch das Oberflächendrücken bewirkten elastischen lastischen und 


> 
E45 


‚Richtung in Abhängigkeit von r besteht. SER 
2. Um die theoretische Durchführung zu erleichtern, wird die Poissonsche Konad Be 
m = ® gesetzt. Es wird ferner angenommen, daß die Länge ! des Stabes groß sein soll gegen- R 5 + 
über dem Durchmesser d, so daß 1:d— ® gesetzt werden kann. Die Störungseinflüsse, die 
an beiden Stirnflächen des Probestabes auftreten, fallen unter»dieser PRSMSEER, aus Gil Br | 
Betrachtung heraus. 
3, Es wird ferner ein einfaches Gedankenmodell vorausgesetzt, bei dem der Se in Fer, 
plastisch nicht verformten Kern a und in eine äußere plastisch gleichmäßig verformte Schichttdb 
4 unterteilt wird. Es wird ds 
elastische Spannungssystem, das ; 
zwischen beiden Teilen zum 
Zwecke der Vereinigung erzeugt 
werden muß, berechnet. 
4. Der Stab wird mit dem 
Radius r ausgebohrt und die 
„ verhältnismäßige Durchmesser- 
änderung z, bzw. die Spannung 
co; in Abhängigkeit vom Bohr- 
radius r aufgetragen. & = fı (r)- 
Die Spannungen o,, die mit der 
® Querschnittsfläche multipliziert 
: |" die Kraft in axialer Richtung 
ze geben, werden in Abhängigkeit 
Er l, von 7? aufgetragen. u = fa (r?). 
ER a a a oe Man erhält auf diese Weise 
, EEE £E: fr) ee zwei Hyperbeln für den plastisch 
4 RE 00 BES nicht verformten Kerna mit 2 
Mr EA. . aufeinander senkrecht stehenden  \ 
re EEE Asymptotenachsen X und Y. 
0 Lat OT un er BER REN Die eine Asymptotenachse = 0 
Pong ‚oder r=a-b gilt unter der 
hat} 2, Voraussetzung, daß die äußere 
0 / Schicht b unendlich dünn ist, für 
E die Stelle r=a--b, wobei a 
der Halbmesser des Kerns und b die Dicke der plastisch verformten Schicht ist. Die 2. Asymp- 
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en 


be 
„inne nur. 
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totenachse mit - = y= 0 bezieht sich auf den Fall, daß die Kernfläche «a durch außerhalb 
0 


beigefügten plastisch nicht verformten Werkstoff unendlich erweitert wird. (Siehe Mitteilungen 
des Wöhler-Inst., Heft 42.) 

5. Es sind die Meßpunkte für die verhältnismäßigen Längen und Durchmesser des Probe- 
stabes in Abhängigkeit vom Bohrradius r in Bild 1 aufgetragen worden. Zur = (0 (Stab ohne 
Bohrung) gehört der verhältnismäßige Längenunterschied y, zwischen dem oberflächengedrückten 
und ungedrückten Probestab. Zum Bohrradius r, wird der verhältnismäßige Längenunterschied 
Y„ gegenüber der Länge, zugehörend zur =0 bestimmt. Die verhältnismäßige Verlängerung %,. 
des Kerns infolge der plastischen Verformung der Schicht b ist also y, — YA Yo. 

6. Solange sich das Ausbohren auf den plastisch nicht verformten Kern a beschränkt, 
liegen die Meßpunkte auf den beiden Hyperbeln, die durch zwei Meßpunkte gelegt werden können. 


t 


190.4 


die ine M—= © weglassen und erhält die größten Druckspamnungen in der ober- 
_flächengedrückten Schicht zu: ; Re: 


JEp FAN ke: 2 A ° 


a (eu je ME); Tao = Orr, (end + m ea) 


v2 0 = me 


m diesen el nehmen die nen O0 Und oas0' in Karen aler Ba aafaier Rich- 
tung in der Schicht 5 bei Federstählen und entsprechenden Drücken Werte von etwa 90 kg/mm?® 
5 ‚bzw. von etwa 160 kg/mm? an. Die Werte stimmen mit den Werten ungefähr überein, dieman 
Be: beim Abätzen der äußeren Schicht durch Mare Längenmessungen bestimmen kann. A 


Über den en ieycland beim Koselen von Flä chenträgern 


Von A. Pflüger in Hannover | R TRAEE ee 5 


Br Bei Beuluntersuchungen an Platten -und Schalen ist ‚es FB und RT E 
Spannungen durch Integration über die Wandstärke zu 'Schnittkräften bzw. -momenten je 
e Panzsneinheit der entsprechenden Koordinatenlinien der Mittelfläche zusammenzufassen. 
#7 Man “ermittelt so zunächst — am unverformt gedachten — System die Schnittgrößen des 
B Grundzustandes, d.h. des Zustandes, der bis zum Ausbeulen gilt, und untersucht dann, für 
welche Belastung es einen brauchbaren Gleichgewichtszustand gibt, der dadurch zustande 
kommt, daß die Kräfte des Grundzustandes infolge der Verformung ihre Richtung ändern, 
- dabei aber mit den durch die Verformung zusätzlich auftretenden Schnittgrößen gerade im 
© Gleichgewicht sind. Diese Betrachtungsweise bringt jedoch zwei Schwierigkeiten mit sich, 
auf die nach Kenntnis des Verfassers im Schrifttum bisher noch nicht hingewiesen ist. 


Erstens betrifft dies die Tatsache, daß man sich nicht nur um die Richtungsänderung 
2 einer Schnittgröße als Ganzes kümmern muß, sondern auch darum, daß sich infolge der Ver- 
 formung ihre Zusammensetzung aus den Spannungen noch einmal wieder ändern kann, weil 
R die verschiedenen parallel zur.Mittelfläche verlaufenden Schichten verschiedene Verformungen 
haben können. Vernachlässigt man diesen Einfluß, wie es allgemein üblich ist, so wird das 
Ergebnis nicht genau richtig. Ohne auf Einzelheiten einer Ableitung einzugehen, sei als Bei- 
spiel angeführt, daß auf diese Weise in der bekannten Differentialgleichung der Beulung einer 
- auf Druck beanspruchten Rechteckplatte konstanter Dicke 
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= ; rt Itr* 2. 1 DR [2 
Be - = u + 2w U ne Ya 
(w Durchbiegung, t Wandstärke, o, Beulspannung, « Querkontraktion, E Elastizitätsmodul) 


auf der rechten Seite das Glied + d— u) w'" verlorengegangen ist, weil der Einfluß der 


Verwindung des Plattenelementes auf die Spannungen des Grundzustandes vernachlässigt 
worden ist. Der begangene Fehler ist jedoch nicht groß, da (1 — u?) = die Größenordnung 


der Dehnung des Grundzustandes hat, die gegenüber der Zahl 2 des Gliedes 2w”* stets klein 
ist. Es ist also schon berechtigt, den in Frage stehenden Einfluß zu vernachlässigen, wenn es 
auch unbefriedigend ist, ihn einfach zu ignorieren. 

Eine zweite Schwierigkeit taucht auf, wenn man bei Beuluntersuchungen an Schalen 
den Grundzustand als Membranspannungszustand voraussetzt, was fast immer getan wird 
bzw. getan werden muß, um die Rechnung mit ertragbarem Arbeitsaufwand durchführen 


.beulen auftretende Beanspruchung des gesamten Querschnitts ist dann nach der Biegetheorie 


den Grundzustand eine konstante Verteilung der Spannungen über die als hi eichend klei 2 
vorausgesetzte Wandstärke, wobei die Biege-, und Drillmomente verschwinden; in der 5 
nauen Biegetheorie, die zur Berechnung des Ausbeulens angewendet werden muß, gehören 7 
aber zu den konstanten Spannungen wegen der jetzt nicht mehr beliebig kleinen Wandstärke 
und der infolgedessen zu berücksichtigenden Trapezform der Schnittfläche eines Elementes er 
doch wieder Schnittmomente. Da diese jedoch lediglich durch den Widerspruch in der Rechen- 
methode zustande kommen, keinerlei mechanische Bedeutung haben und zahlenmäßig belang 
los sind, ist auch ihre Vernachlässigung berechtigt, was ebenfalls allgemein üblich ist. Er 
Die beiden genannten Schwierigkeiten lassen sich am ° 
einfachsten aus der Welt schaffen, wenn man nach dem neben- 
stehenden Bild statt des wirklichen Flächenträgers ein anderes 
statisches Gebilde von besonderem Aufbau und besonderen 
Eigenschaften voraussetzt. Man muß dazu eine Schale an- 
nehmen, die aus verschiedenen Schichten besteht, und zwar aus 
einer sehr dünnen Schicht von der Dicke t, und zwei weiteren 
_ außen liegenden Schichten von der Dicke {/2—t,/2,so daß wieder _ 
| - die Gesamtdicke t herauskommt. Die Kräfte des Grundzustandes 
sollen nun lediglich in der Schicht mit der Dicke t, wirken, während erst beim Ausbeulen auch 
der übrige Querschnitt in Anspruch genommen wird. «Man muß sich dazu vorstellen, daß 
im Grundzustand nur die Schale von der Dicke f, vorhanden ist. Diese Schale wird nach der 
Membrantheorie gerechnet und im übrigen ihre Verformung gleich Null gesetzt. Erst wenn 
sich der Grundzustand eingestellt hat, werden auch die beiden übrigen Schichten aufgebracht, 
so daß diese in der Tat spannungsfrei sind, solange die Schale nicht ausbeult. Die beim Aus- 


zu können. Es ergibt sich nämlich dann im Rahmen der 


zu rechnen, Dabei soll es aber durch Wahl eines hinreichend kleinen Wertes vom, berechtigt 
sein, nur in Z, lineare Glieder mitzunehmen, so daß man bei allen in der Mittelschicht wirken- 
den Kräften im Rahmen der Membrantheorie bleibt. 

Durch das besprochene Modell wird nun zweifellos .das Verlangte geräde erreicht: 
Erstens wirkt jetzt der Grundzustand in einer hinreichend dünnen Schicht, so daß die Unter- 
schiede in den Verformungen der einzelnen parallel zur Mittelfläche liegenden Schichten den 
Grundzustand nicht mehr beeinflussen können. Zweitens wird es möglich, Membrantheorie 
und genaue Biegetheorie gleichzeitig widerspruchsfrei nebeneinander zu benutzen. Durch 
Einführung des Modelles wird selbstverständlich an der üblichen Rechnung im Endeffekt 
nichts geändert. Es. wird nur klargestellt, was man eigentlich zu rechnen pflegt, welche Ver- 
nachlässigungen man dabei macht und inwiefern diese berechtigt sind. 
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Probleme der partiellen Differentialgleichungen der Mechanik 
der Kontinua 


Von E. Hölder in Leipzig 


Der Vortrag behandelt in Verallgemeinerung eines 1941 auf der DMV-Tagung in Jena 
vorgetragenen Variationsproblems ebenso die stationäre adiabatische Strömung im Raum ° 
eines Mediums mit dem Mollier-Diagramm für den Druck 


P=9l,s ww); 
die rechte Seite hängt unter der Annahme chemischen Gleichgewichts außer von dem Wärme- 
inhalt © und der Entropie s pro Masseneinheit noch von den k— 1-Molenbrüchen Y, der k-Kom- 
ponenten ab. Der Wärmeinhalt ist nach der Bernoullischen Gleichung 


1 
ß > 2 ee 
>) Ux + (= % 


durch die Geschwindigkeitskomponenten u, (a = 1,.2, 3) auszudrücken, diese wiederum durch 
zwei Stromfunktionen 9 (2,),X (%,) in folgender Weise: 

Ist die Geschwindigkeit im Unendlichen (U,,0,0), so drücken sich durch die dort 
asymptotisch angenommenen Lagrangeschen Koordinaten b,, b,, welche die zweiparametrige 


Zwar 29, 


’ $ 0: FR FE oainktionen; = v. 0. E8 Fa = Vermöge dieser Gleichungen Be, 
2 nen %, durch die ersten Ableitungen von 9, x ausgedrückt werden. Das  Variationsproblem R2 
% stationären u ist durch das Extremalintegral Era Meg 2 PR 


Er & Gag enden stationär 


+ 


Er De Ferteäui=r e Bpun): ‘ 
Y gegeben, | N EN ; 

0.9, x und x, selbst Aretenlein: das ei b,) auf jeder Stromlinie SR ist bis Auf etwaige. Re 
_ Verdichtungsstöße; i, kann außer von b,, b, beim Vorhandensein eines Außenfeldes auch 
noch von den x, abhängen. An einem Verdichtungsktoß hat man außer den bekannten aus 1 
 Energie- und Impulssatz und der Kontinuität der Gesamtmasse folgenden Relationen noch 
1 weitere Massenbilanzen, welche die Übergangsrelationen auch der y, festlegen. 

Zum Schluß wird das zweidimensionale Variationsproblem für den rotationssymmetri- 77 
schen Strömungsfall erwähnt, das im Weickmarn-Heft der Meteorologischen Rundschau 
erscheinen wird. Hier Herrscht in jedem Meridianschnitt derselbe Geschwindigkeitszustand 
u (2, Y), v (x, y); y Abstand von der Achse, x Koordinate in Richtung der Achse. Es besteht 
aber noch eine Geschwindigkeitskomponente w (z, 4) senkrecht zur Meridianebene, unab- ne, 
_  hängig- von der geographischen Breite. In diesem rotationssymmetrischen | Fall hat man 92; 
ebenso wie im ebenen und axialsymmetrischen Fall nur eine Stromfunktiony. Die Grund- 
funktion des zweidimensionalen Variationsproblems ö [[ F dx dy= — ( ist 


Br : = yß+e@+W%))=F (m, 9, 9 Ya yı)- 
. Dabei ist die Bermoullische Gleichung ER 
Er 1 22 (y 
En Zt ien N eh  Ue N) 
zu verwenden, wo 2= yw= Q(y) das nur von y abhängige Rotationsmoment und U(z, y) 
die. potentielle Energie pro Masseneinheit ist,-endlich | 
| Yyou—yy YOUV—=— Y%- 

Das Variationsprinzip liefert ebenso wie das in ihm enthaltene von Friedrichs (Math. 4: 
Ann. 109) für inkompressible Strömungen (ohne Rotationsmoment) äuch freie Strahlgrenzen, BE 
ganz allgemein die rotationssymmetrischen Diskontinuitätsflächen der atmosphärischen Dy- 


namik. Auch die Variationsbetrachtungen von Helmholtz (Wiss. Abh. III, S: 333) über Br 
„Die Energie der Wogen und des Windes“ finden ihre Verallgemeinerung aufs Kompressible. a 
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Bemerkung zur Lösung von Elastizitätsproblemen von vorgegebenen ee 
Randverschiebungen bei Scheiben | er 


Von Udo Wegner‘ in Heidelberg 


Zur Berechnung von Spannungen in Scheiben bei vorgegebenen Randverschiebungen, 
bzw. von Verschiebungen und Dehnungen benutzt man größtenteils die Methoden der Theorie 
'komplexer Funktionen, der Greenschen Funktion oder das Ritzsche Verfahren. Die Me- 
thoden der komplexen Funktionen, deren Verwendung bei mehrfach zusammenhängenden 

- Bereichen fast unumgänglich ist und deren Begründung auf Kolosoff!) und Muschilis- 
vili2) zurückgeht, eignet sich auch zur approximativen Berechnung der Komponenten des 
Verschiebungsvektors, bzw. Spannungsvektors, wie früher?) dargestellt wurde. 


1) @.Kolosoff, Über sinige Eigenschaften desebenen Problemes der Elastizitätstheorie. Zeitschr. f. Math. u. Phys. 62, 


1914. - 
2) N.Muscheli&vili, Praktische Lösung der fundamentalen Randwertaufgaben der Elastizitätstheorie. Z,angew, 


Math. Mech. 13 (1933). 
U,Wegner, Anwendung der Methoden der komplexen Funktionentheorie auftechn. Probleme, 2. Abschnitt, $3, 
a demnächst im Verlage Teubner, Leipzig. 


r 


mittels der Sätze über Laplacesche Integrale und deren U = a 


‚und der Spannungen längs der Kerbe. (Dies wurde auf der vorigen Versammlung der ( 


- recht kompliziert bei nicht einfach berandeten Bereichen sowie bei der Berücksichtigung von 


Bei unendlichen Ben die’ von BR Kreisbögen begrenzt sinc 


nungsprobleme in geschlosserier Form angeben, selbst GR dem ei ndbogen die 
Verschiebungen und auf dem anderen die Spannungen vorgegeben ‚sind, 2. B ä ür die 
kreisförmig gekerbte Halbebene bei Vorgabe der Verschiebungen längs des Geradk i 


in Karlsruhe im Frühjahr 1947 gezeigt.) 2 Re 
Die Methode der‘ Greenschen Funktion em: Treiftz u.a. )*) eignet sich E - 
tischen Durchführung — selbst bei Beschränkung auf approximative Sp = 
— fast nur bei einigen sehr einfachen Scheibenprofilen (z. B. Kreis-, Rechteck-, Ellipsenpı iR 
Letzten Endes greift man wieder auf das Ritzsche Verfahren zurück, das durch das Costa 
sche Konvergenzerzeugungsprinzip besonders verbessert wurde?). Die Annäherungsfunktionen 
werden aber bei dem Ritzschen Verfahren in den zu berechnenden Größen der Form nach 


eventuell vorhandenen Symmetrien, was zur raschen Konvergenz des Verfahrens sehr AREEBEN 
mäßig ist. Man kann aber auch mittels mehrfacher Anwendung der Friedrichsschen Trans- 
formation®) das zu den Grundgleichungen der ebenen Elastizitätstheorie gehörige Variations- 
problem so umformen, daß man mit Potentialfunktionen und anderen nicht den Rand- ni 
bedingungen genügenden willkürlichen Funktionen als Annäherungsfunktionen operieren en 
kann, um die Werte der Komponenten des Verschiebungsvektors zu erhalten. a 


Problem: RER ERBE Le 
1—v9x "4% = rechtwinklige RRERRAEEN 
1+v20 des Verschiebungsvektors - 
Av we 
1—v0y 
ou 9v er 
G= N — Querzahl, ® = einfach zusammenhängender Scheibenbereich, I’ sein 2 
Rand, 4%, %= Randwerte, s= Bogenlänge der Randkurve. 
vesigd.,. Be gle (f und g stetig) ar 
seien vorgegeben. > B 
Dann ist x 
| du 08 Du ed 
v= ®-+- const und En Br Ep Erf, 


NS he ‘ 
Dabei sind F(z, y) und — ro @(x, y) der Realteil und Imaginärteil einer in ® regu- 


lären analytischen Funktion (auf /’ stetig), die das folgende Integral zum Extremum 
machen: 


L. ee dedy 


+ else ir tete 


P (x, y) ist ein in ® stetig differentiierbare Funktion, die das Integral: 


1 oD12 x © 41 28 98 ob 
Ir, =] (5 ) —2(# TR 5)|azay - 0320 
B 


zum Extremum macht. 


(Man kann auch, was manchmal zweckmäßiger ist — s.u. — eine stetig differentiier- 
bare Funktion Y% et so daß: 


n) 6) Oy 
 ® +2 (r reifjner. Joe Extrem. ist, 


9 Handl Hanalun der Physik, Bd. VI (1928), 8. 122. 


»)R.Öourant, Über ein konvergenzerzeugendes Prinzip. .., Göttinger Nach 
*) K.Friedri chs, Ein Verfahren der Variationsrechnung. ""Adbtinger Nachrichten Uileey 
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| e a, 
und ee man in II D it: ein. © und y Sr man un Dh 


on n en Bes eventuell noch den durch das Problem T aegeBehen = 
a 3 N “ 


ngun | 

Ir Bei: = ir IT REN: “E nt 
DEREN TEE TE Fuss Rechteck; —a<ıza; ee BE IRH LTE SR 
Mey ne 


Fr en | 2a = = n a=H+a 2 m | % Fe So v 9) h 
“ EEE ee ISE vgl) für ER y=+b. - re 
i 2 
u =4g(s) hl s) genügt. Auf y=H4 b braucht sie keiner ‚Bedingung zu ‚genügen. Ebenso. 2 
ist Y auf I’ keiner Bedingung unterworfen. % wird aus II. ‘bestimmt. 


0... Auch in.der Theorie der Platten und der Zylinderschalen sind die oben genannten 'Approxi- 
 mationsmethoden sehr schlagkräftig, besonders auch bei den durch andere Methoden fast 


"ist. dann eine > reguläre analytische Funktion, die auf x = + a der Bedingung 


Br '  unzugänglichen sog. statisch unbestimmten Rechnungen, wenn z. B. die Zylinderschale an 

Be - ‚Randträger meet wird ?). | 

Di 3 3 

Ber; €. Angewandte Mathematik 

2 _ Über komplexe Nullstellen von Zylinderfunktionen 

ge: 2 Dr EN a -_ Von F. Emde in Stuttgart ' 

2. >,(Eine ausführliche Arbeit über diesen Gegenstand wird in einem späteren Heft dieser Zeit- E 
; es > r@gschrift erscheinen.) Se 


Quasistrenge geschlossene Horkiele für elliptische Normalintegrale, € 
Thetafunktionen und elliptische Funktionen = 


n 1 Von F. Tölke in Karlsruhe : Ja 
(Ein Referat ist nicht eingegangen.) 


Transformationen der Variationsrechnung und ihre Anwendungen 
auf technische Eigenwertprobleme 


— Von H. Schaefer in Braunschweig 


RR Eeetinieche, Eigenwertrechnungen führen gewöhnlich auf ein Variationsproblenr (V. P.) Be: * 
der Gestält 


2 
Be: [1 
an a [au nan a zn. sn Ams) 1e=0 en Ai 
e 


"U.Wegner, Ein neues Verfahren zur Berechnung der Spannungen in Scheiben. Forschung auf dem Gebiete des 
Ingenieurwesens, Bd. 13, Nr. 4 (1942), — U.Wegner, Eine neue Methode zur approxmativen Lösung von n Spannungsproblemen 
bei Platten und Scheiben. Forschungshefte aus dem Gebiete des Stahlbaues, Heit 6 (1943). 
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3 at dem E g nwer tpa; meter A. I 1 | Su 
 raden Balkens Laneeh, en { Si e| f >q 
27% hin, um An der folgenden Methoden erkennen re Hieraet. = 

ER Wir setzen aa ae Bee >. = 


a ee 
Ä und as Biegemoment pR= ame ein und bilden die Hamiltonsche Funktion | : 22 Ä : 
H (Pr Pas u 4 DSnitnd em wi 

- 445 1 (Brusgsinn)| 


‘ 2 Er . a 


Be Die ET er schen ge des ee Varsationsproblems" 
Mn - Fi \ £ ' 2 Mn ai ah '$; 
BE ar | msrantäpen Pa; 1» 92)} d=0..... RE - 8), = Pr 
Br; | in dem jetzt 9» 9» P,, 9, unabhängig voneinander variiert werden, sind die vier 
Bis kanonischen Gleichungen“ er ae 
er: Le 1 
nn 8% r EPA Iı j 2 
rstde) ERBE 
hr „= OEES2.:95, ER, Fan En 
:: 4, °P EJ | P, 09; En P: 4»4g 
2 (4 a) liefert den A an nEDnE zwischen inneren Kräften und Verschiebungen, (4 b) die 


Gleichgewichtsbedingungen or 
Wir fügen nun zu (3) als bee ante (4b) hinzu und erhalten EL 


+[B- #84 -Amat)de=0 BE N Re“ (5) 
ö 


mit den beiden Möglichkeiten, entweder mit (4b) in (5) q, und g, durch p, und 9, zu er-. 

, setzen und somit ein V.P. für p, und p, zu bekommen, o der durch Integration von (4b) p, 
durch g, und g, auszudrücken und in (5) einzuführen, um ein V.P. für q, und 9, zu erhalten. 
Beide V. Pe. erstreben die Erfüllung von (4 a). Im Verfolg der zweiten Möglichkeit wird {; M(z, r) 
Biegemoment durch Last ? in r): 


ı 
oM (z, 
me) = afarıe moqdrtu Partys 
. und 


I 
Pla) A ( M(x, r) m (r)gq, (r) ar) . ([ M (x, t) m (t) q,(t) a)+ USW,, 
also 


1 
[E5ge=% |I®« t) m (r) qı(r) m(t) q, (t) dr di + usw. 


) 


r) M (x, t) 
TE J (x) 


M (x 
mit dem symmetrischen Kern X (r, t ji %, 


so daß (5) die Gestalt 


2 2 
Rdn +2 Ahut au [Sat dem made) a N (6) 


‘) Zur Bigenwerttheorie derlinearen kanonischen Systeme vel, E.H ö1d er, 


und Stabilitätsproblem. Jahrbuch 1940 der deutschen Luftfahrttorschung. Dort weiter Banane als funktionale Verzweigung 


weitere Schrifttumshinweise, 


a er ? R 
(4 a) zu (3) als: Nebenbedingung a so erhält man’ a des ursprüng- 
NSD% da). oder ein Integralgleichungs- ‘V.P. für p,. Ein V.P., das g,’ = q, erstrebt, ent- 
, wenn die drei restlichen kanonischen Gleichungen als Nebenbedingung Er werden, | 
a au mit FR (Y) = ENT. Ba nd Br : 


offen er un): BE RER % 

; er P% 
Auf GE Möglichkeit, mit Hilfe, Kenonfacher Transformationen, ‘ wie ‚sie An. de an er 
RS _ Mechanik gebräuchlich sind, neue  Variationsprobleme, zu gewinnen, möge hier noch ie E “ 
Ki ‚ wiesen werden. Een - 


BC: DE. Über Abschä ätzungen von Lösungen Inhomogener 


Bi: = ‚Differentialgleichungen ° 

Eis ER Er ‚Von Hans Büchner in ‚Göttingen Br 

Br 923 = "Es Bench sich um die Differentialgleichungen ; 
N ee N uEacy 

BR 2 a RE ZEN, AR 

E23 > SR, EUR ds 0 R 

2 als konstanten Koeffizienten. ‚Für jedes Intervall a<x< b läßt sich mindestens . ein partiku- 
E . läres Integral finden, welches der Abschätzung 

en ee, IS C- Max fl 


u 


Sunterliegt- Für die a C darf gesetzt, den 


H FT. O% 
2 ir etgh en 


| wobei 

” - ' + 8% 
die Paare konjugiert Klier Wurzeln der charakteristischen Gleichung darstellen. Hinweis 
auf Differentialgleichungen stabiler Regelschaltungen und Servomechanismen. 


Einige Bemerkungen i über die Orthosonalpolynome zu einem endlichen 
Intervall mit einer vorgeschriebenen Belegung 


Von Hans @ebelein in Braunschweig 


Vielseitige Anwendungen der Orthogonalpolynome zu einer beliebigen Belegung w(x) >0, 
wie ich sie im Hinblick auf statistische Probleme in meinem Buche ‚Zahl und Wirklichkeit‘ 4 

dargelegt habe, legten mir den Wunsch nach einer Methode nahe, die es gestattet, unabhängig 

von der speziellen Wahl der Belegfunktion die SRtuszanelpol none vor allem hohen Grades in 

| ihrem gestaltlichen Verlauf zu überblicken. 

x Sukzessive Orthogonalisierung (siehe Courant-Hilbert, Method. d. Math. Phys., 
Bd. I S.34, 1. Aufl.) auf Grund der Beziehungen | u: 

Ja" mn( x), wix) dx = 0... - für re SE Re er ee 

2) Vgl.R. Grammel, Über die Lösung technischer Eigenwertprobleme, Forschungshefte aus dem Gebiete des Stahl- 


baues. Heft 6, Berlin 1943. Diein Sisasy Darstellung implizit benutzten V. Pe. lassen sich ebenfalls aus einem kanonischen V.P. 
herleiten. 


ee 
VE FREE AH | 
San 2 ee > 
FIR "h © 


Ara PR (ner 
a r [3 \ 
re a € . 


Rd - BE ya 


% N, 


scheidet für diesen Zweck ebenso aus wie der Gebrauch von R 
EN | Pur = (®— Ba) Pa — On Bn-ı 
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Be Do 30 OR, Ä ee 
0° dahier als Ausgangspunkt nicht eine rekursive, sondern eine independente Darstellung für 

A erforderlich ist. Anderseits ist die bekannte, auf Tschebyscheff zurückgehende Dets 

Be, minantendarstellung | | re Et ns 753, 

Bu F . L328 5,239 MT 2 : 
Br: \  Pu(e) = const he 1 PBr3 

I. = j M,„-ı M, Mal Ms .-ı Pr ar “ 
Be nicht geeignet, weil der Rechenapparat zu schwerfällig wird, sobald man über die allerersten - 
r : Polynome hinausgehen will. Als für unsere Zwecke brauchbar erweist sich hingegen die komp E 
ee, mierte Wiedergabe der Orthogonalitätsbedingungen (1) in Gestalt des auch von G. Szegö Be | 
P _ in seinem Buche „‚Orthogonal Polynomials“, New York 1939, angeführten Variationsproblems DS 

hr J=/p (2)-w(a)de= Minimum . . ...- 2... 0. (4), Se | 


wobei für p, alle Poltnarie n-ten Grades mit dem höchsten Koeffizienten Eins zur Konkurrenz 
stehen. Schreibt man diese Polynome in Produktdarstellung, so lautet das Problem 
e 


I=) I ia? wi) d& = Minimum a er RE (da), 
m=1 


und die Aufgabe läuft darauf hinaus, die n Wurzeln des in Rede stehenden Polynoms richtig 

zu legen. Gleichung (4a) eignet sich auch besonders gut dazu, in einfachster Weise zu zeigen, 
daß alle Wurzeln reell sind, keine außerhalb des Definitionsintervalls liegt und Doppel- oder 
Mehrfachwurzeln nicht vorkommen können. Ausgeschrieben lauten die Minimumsbedingungen 
bezüglich der Wurzeln (2J/d 2, = 0) 


Die im folgenden mitgeteilten Ergebnisse sind an die Voraussetzung eines endlichen Defi- 
nitionsbereichs gebunden, von dem ohne Einschränkung der Allgemeinheit angenommen werden 
kann, er erstrecke sich von —1 bis --1. (Für den Fall eines unendlichen Bereichs, bekannt 
durch die 

Hermiteschen Polynome mit uwu@)=c", -w<a<+tx 


und die Laguerreschen Y »’ BD =- Er, O<I<CH - 


konnte ich bisher keine wesentlichen Ergebnisse erzielen.) Auf Grund der Sätze über die Wurzeln 
läßt sich vorwegnehmen, daß der Verlauf des Orthogonalpolynoms n-ten Grades zwischen —1 
und -+-1 ein Wellenzug ist, bestehend aus n—1 Halbwellen im Intervallinnern und zwei Ansatz- 
stücken, deren Gestalt davon abhängt, wie die äußersten Nullstellen zum Intervallrand liegen. 

. Das Augenmerk ist daher zu richten erstens auf die Größe der Amplituden, zweitens auf die 
Verteilung der Nullstellen und drittens auf das besondere Verhalten in“der Nähe der Ränder, 
Für diese drei Befunde gelten asymptotische Gesetzmäßigkeiten, die um so genauer erfüllt sind 
je höher der Grad n ist. Jedoch zeigt die,rechnerische Nachprüfung bereits bei recht niedrigen 
Graden erstaunlich gute Übereinstimmung. 


I. Amplitudenverlauf: 
Das in Rede stehende Grenzgesetz für die Amplituden lautet 


Ya) = Amp (910) Yo TO) = 7 ante er SE nane *tb); 


I= [pe w)de= 0 ES (4b). 


Beispiel: Für das relativ niedrige Le gendresche Polynom P, schwankt die Kon- 


- stante C nach (6) noch zwischen 0,2910 in Intervallmitte und 0,2884 für das abseitigste Ex- 


tremum. Mit C=0,2910 gerechnet, folgt J=0,1330 gegenüber dem exakten Wert 2/15 = 0,1333. 
Über den Beweis sei vermerkt, daß die bekannte, aufdie C auch y schen Hauptwerte sich 
beziehende und für beliebige & zwischen —1 und 41 sen Integralformel 


TE Maar 


eine entscheidende Rolle spielt, wodurch es verständlich wird, warum der Gedankengang sich 


nicht ohne weiteres auf den Fall eines unendlichen Intervalls übertragen läßt. Die Abschätzung - 
des Integrals (4) geht so vor sich, daß für jede Halbwelle zwischen zwei aufeinanderfolgenden 


Nullstellen des Integranden von der für einen Sinusbogen streng geltenden Formel 


Lo —% 1 
j Y°(z) de = I Due ag >s I Yılaz 


Gebrauch gemacht wird. 

Die im Verlaufe der Rechnung auftretende Funktion Y1— 2? legt den Übergang von der 
Strecke —1,--1 zur Peripherie des darüber liegenden Halbkreises nahe. Durch diese .aus der 
Theorie der Legendreschen Polynome geläufige Transformation 2= cosp nimmt das 
Variationsproblem die Form 


J= mo. w(g) sinpgdgp = Minimum . . N. EN (4e) 


an, und das Amplitudengesetz sagt aus, daß die Funktion p,-yw-sing einen Wellenzug mit 
n Nullstellen im Intervall (0, x) und mit sehr nahezu konstanten Amplituden darstellt. 
II. Lage der Nullstellen: 


Je höher der Grad z, um so besser erfüllt das Polynom 7, das folgende, für die Anord- 
nung der Nullstellen maßgebliche Grenzgesetz Ä 


TE) sny = less (Ap HB) . : 2 af nenne (N. 


Aus dieser Beziehung folgt, daß die Nullstellen der Orthogonalpolynome, ausgedrückt ing mfit 


zunehmendem Grad mehr und mehr eine arithmetisehe Reihe bilden. 

Beispiel: Die Nullstellen des Legendre schen Polynoms P, liegen exakt bei 
op = (°, + 23056’, + 47052’ und + 71038’. Die Theorie liefert auf Grund des Wertes ö = 0,75 
(siehe unten Gl. 13) für die Nullstellen-C6°, 24°, 48° und 72°, mit Z = 0,761 nach Gl. (12) aber 
sogar 6°, 23055,8°, 47° 51,6’ und 71° 47,4. 

Es ist verhältnismäßig leicht einzusehen, daß die weniger weit gehende Aussage, wonach 
irgendeine herausgegriffene Nullstelle mit wachsendem Grade immer genauer mit dem arith- 
metischen Mittel der beiden Nachbarn übereinstimmt, zutreffen muß. Daß aber durchgängig 
im ganzen Definitionsbereich die Nullstellen gegen eine arithmetische Reihe streben, liegt an- 
scheinend tiefer. Hier sei nur angemerkt, daß jener akustische Vorgang; bei dem durch Über- 
lagerung wenig gegeneinander verstimmter Schwingungen hoher Frequenz Schwebungen vor- 
geschriebener Modulation erzeugt werden, ein Gegenstück zu dem vorliegenden Befund darstellt. 

Für die zeichnerische Darstellung der höheren Orthogonalpolynome ist es zweckmäßig, 
außer den Nullstellen und den Extremwerten auch noch die Neigungen p, (2,) an den Schnitt- 
punkten mit der Abszissenachse zu kennen. Für sie gilt, wiederum asymptotisch, das Gesetz 


Le Yolz,) = const .. .......82 


II. Das Gesetz für den Verlauf in Randnähe: 


Die oben besprochene Transformation vom Definitionsintervall auf die Peripherie des. 


darüber liegenden Halbkreises bewirkt, daß für randnahe Punkte zwischen dem Abstand & 


von einem der Ränder und dem zugehörigen Abstand des Winkels p von beire 
ein quadratischer Zusammenhang besteht. Im folgenden wird unter = - dieser rel 
ringe Abstand vom Rande verstanden und unter 2 der geeignet normierte Winkel, der ist 
falls vom Rande aus gerechnet wird, so daß «= 2? ist, und die Nullstellen des Polystar 

geschrieben, lauten: IE | DEE EN, 


2 a=C +1 1, Ve FE 


Es wird der Grenzfall eines Polynoms von sehr hohem Grade betrachtet; die Folge (9) erstreckt 
sich also ins Unendliche. Das in Rede stehende Polynom p(z) wird außerdem so normiert, daß 


»(0)=List. Ya 
Nach diesen Festsetzungen ist der einzige Parameter des herausgeschälten Problems Z, 

und das Ergebnis ist die transzendente Funktion 
(1 2 II: (£) x 
ee Pre | 
Sie stellt das allgemeingültige Grenzgesetz für den Verlauf der Orthogonalpolynome in Rand- 
nähe dar. An Stelle des Parameters © kann mai auch die Ordinate des ersten Minimums zur 

Kennzeichnung der einzelnen Funktionen der Schar benutzen. 


(23-m-1) Entscheidend ist aber für den Parameter £, 
m 2 daß er in engem Zusammenhang mit dem Ver- 
& Y\\ı83:75 halten der Belegfunktion w(z) in der Nähe des 
Randes steht. Es gelte für w(z) in Randnähe, also 

= 3=125 für kleine x das Potenzgesetz 
40 \ x w(x) = const. =*, (k<—1).. (ll). 
E Das Hauptergebnis, das dann aus (10) unter 
= \ Berücksichtigung der zuvor erörterten Gesetz- 
3, mäßigkeiten folgt, lautet: Eine exakte Zuordnung 
e \ zwischen den Funktionen »(z,£) einerseits und 
10 5 \ den Potenzgesetzen w(x) = z* anderseits würde 
«dann bestehen, wenn für konstante & und %k, aber 
() (\ beliebige ganzzahlige m eine Beziehung der Art 
e \\ Ve Re - (e* F für m-1,3,3, ... 
N 
gelten würde. Diese Beziehung ist allerdings 
z x N keine Identität. Nur für die beiden Wertepaare 

1 1 

/ En \ arm u, Ku und {/=1, un 


Ss Ss 
STR 


5,20 
. gilt (12) identisch für alle m. Es handelt sich 
hierbei um die Tschebyscheffschen Poly- 
nome erster und zweiter Art, deren Sonderstel- 
lung wie an vielen anderen Stellen der Theorie 


3:025 auch hier eindrucksvoll zutage tritt. 
and Obwohl nun zwar die Beziehung (12), von 
diesen beiden Ausnahmefällen abgesehen, nicht 
streng zutrifft, kommt sie doch der Wirklichkeit bemerkenswert nahe, wie man erkennt, wenn 
man für bestimmte und für variable ganzzahlige m die Binomialkoeffizienten 


? + BEN Si 


und die Größen 


* 


mM 


doppelt logarithmisch gegeneinander aufträgt (Bild 1). Es sind in der Tat die Kurven & = const, 


von einem Geradenbüschel kaum zu unterscheiden. Einfach wird der Zusammenhang zwischen ' 


& und k nach (12) für sehr große m; er lautet 


3 ; 
kd)=- az -oder Lk) “715 ee ehr“ (13). 


er > ee 


DE en) 
x 


cc. A u 


n man aus dio) Br M= 1 geninn, Er aibt, wie 
ellen erstaunlich gut Wieden. ya me I: 
AR : e. Verhalten i in Randnähe läßt sich RN zusammen- 
fassen: Der \ erlauf der Orthogonalpolynome in der Nähe des. Randes ist vom Grade des be- 
treffenden Polynoms im wesentlichen unabhängig bis auf eine Änderung . des Abszissenmaßstabs. 
ängt vielmehr davon ab, welches ‚Potenzgesetz für die Gewichtsfunktion w (z) bei An- 
äherung. an den Rand gilt, ‚Man kann aus der Lage der ersten. beiden Nullstellen zum Rande, 
oder auch aus dem Verhältnis zwischen der Ordinate des ersten Extremums und derjenigen. ee 
u des Bereichs auf den Potenzexponenten der Belegung schließen, und 


Lese 


Zusammenfassende Näherungstormele 


. Re 


Die Ergebnisse aller drei i Fragenkreie lassen sich zusammenfassen mittels der Formel 25 


Aser 
A IENEEE a . cos ar Be 4) IP ae } 
ER BER E x = Pr (P) 72. Te 


vw (p: ‚sinp 


“ 


Ener Derisuten n A Grad des durch » anesnaberk wiederzugebenden Grioenaltei e 
 .noms, w(g) die Belegfunktion, £, und £, aber die das Verhalten: der Belegung an den Rändern 
. charakterisierenden Parameterwerte. . GR 
Der explizite Ausdruck (14) ist, in = cos ) geschrieben, im allgemeinen eine BE, Er 
.dente Funktion mit Unendlichkeitsstellen an den beiden Rändern, und kann daher das verlangte 
Orthogonalpolynom n-ten Grades nicht exakt wiedergeben. ‚Reduktion auf exakte Polynome 
vn nur ein in es beiden Fällen % 
ey Br tt ee = ar h=eb-l, | > 
welche auf die Tscheb yscheffschen Polynome beiderlei Art führen. Im Falle der Le - 
gendreschen Salyacne liefert -Gl. (14) mit 22 = tL,= 3/4 die berühmte Lapl acesche Bes 
Formel. BE 
Durchgerechnete Beispiele zeigten, da man über die explizite Formel Gl. (14) hinaus u 
einer noch besseren Darstellung des Polynoms n-ten Grades gelangen kann, wenn man diese RB. 
Funktion durch Zurückgreifen auf (10) an den Rändern von den Singularitäten befreit, danach BER. 
als Fourierreihe entwickelt und diese mit dem Gliede cos n  abbricht. Die erzielte Übereinstim- A 
mung war so gut, daß bis zu niedrigsten Graden heruntergegangen werden mußte, damit über- 
haupt Unterschiede auftraten, die sich in-Abbildungen noch veranschaulichen lassen konnten. 
' Zam Schluß sei noch eine Folgerung der dargelegten Ergebnisse genannt, die sich auf de BE: 
. anfangs erwähnten Rekursionsformeln (2) bezieht. Wenn das Intervall von a bis d reicht, so gilt en x 
y y r = ” “ > : 3 t y 
(da)? | BE 
i = —— 15% 


% 


Liegen die Intervallränder symmetrisch zum Koordinatenursprung, so konvergieren die 
B,„ gegen Null; ist die Belegfunktion w(x) symmetrisch, so sind die B, sogar identisch Null. 
Die CO, betragen für die Tschebyscheffschen Polynome erster und zweiter Art überein- ; 
stimmend 1/4, während für die Legendreschen Polynome a, 


z N? ee y “ 1 A 
RE rn und daher ER = I ist. 
Für das Beispiel der Polynome zur Belegung w(x) = 1-+- x und zum Intervall (—1, —1) 


beginnt die Reihe der B,, mit 1/3, 1/15, 1/35, 1/63 und & Reihe der C,, mit 0, 2/9, 6/25, 108/441, 
so daß die in Gl. 15) E ecnte Konvergenz gegen 0 bzw. gegen 1/4 sehr eindrucksvoll sich zeigt. 


N De über das s Rosi ER 


das Ben R,(a) gewöhnlich in der OR 1) 


a Kari der Formeln zur n e 
"Von &. Weine in dena DE er EIEEENE 

#, Obwohl ee samktenkeng zwischen Interpe d Taylorschem ! 
„ bekannt gilt, scheint es nicht bemerkt worden zu sein, daß sich das Fehlerglied der ! | 
"in einer Form darstellen läßt, die dem Schlömilchschen bzw. Cauchyschen F 
entspricht und die im Grenzfall zusammenrückender Interpolationsstellen in das Sch 
milchsche bzw. Gauchysche Restglied der Taylorformel übergeht. SEE 

Den Zugang zu dieser neuen Form des Interpolationsfehlers gewinnt man aus einer In 
graldarstellung, die bisher anscheinend ebenfalls noch nicht verwendet worden ist. El, 
1. Bezeichnet man wie üblich die zu einem System 2, »=0,1.. m) von Interpolations- 

stellen gehörenden Lagrangeschen Polynome a L,(e) ” =0,1,.0427 Re BURN es % 

bs Interpolationsformel 


fa) = > fa) Le + Bla en BR Be 1 


2 le fr SE ——din ns a 2) 
v0 - 
oder in der Mittelwertform: 
Er n | 
x R„(«) = L(x) ar) a (3) 


angegeben, wobei & einen Zwischenwert aus dem durch die Stellen 2, 29 24... u besfimmienS 
Intervall und L(x) das ‚Produkt 


(an) re ER ur} (4) 

bedeutet. 
2, Die neue Integraldarstellung für a erhält man durch folgende Überlegung: 
‚Bei festem & stellt 


ern BR Ze ER RE pr nn 


or 


eine im oben genannten Intervall stetige und mindestens n-mal differenzierbare Funktion der 
Veränderlichen u dar. Beachtet man weiter, daß die Lagrangeschen Polynome die Form 


\L(e) 
L a A N 
haben, so geht Gl. (2) über ’'in 
ö g(%,) 5 
Rn (0) = L(8)7 7, Sm s 


Die rechts stehende Shrkane erkennt man als die n-te Steigung, der Funktion g(u), 
und somit ist wegen des Zusammenhangs zwischen Steigungen und Ableitungen 


R,(@) = L(a) 9" (d), 


wobei £ wiederum einen bestimmten Zwischenwert aus dem durch die Interpolationsstellen auf- 
gespannten Intervall bedeutet. 


Durch Differenzieren der Gl. (5) erhält man so die neue Integraldarstellung 


j f MH (Et Se 
Rn(a) = L(a) I em RE a). 


) G.Kowalewski, Interpolation und genäherte Quadratur. Leipzig 1982, 3.24. 
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erhält a das Bass Mia: de Integrakdastellung 2 27 n er ES 
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Hieraus folgt - — Fan der Voraussetzung, dad Sr (a) ei eine ntegierbare In 9)- -te © Ableitung besitzt 
= = durch. a Br | 
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et er 3 ie Se ni 4 ne ne = Be = Es 
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un! ‘ det > \ a EURE 
gi 


“Der Mittelwerlsatz der Inegrerechmung liefert die Formeln: 5 
EI N I(E +2) | 
re ar ER es vLiaf DE N) 


en = nen feel ar 
ee gem). 7 pirtaug) 
ee BE NEE A hose: Ly! L(a) ER RT (VII) 
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% - Harmonische Analyse von \ Kurven mit stark ausgeprägten Spitzen 
= Von H.-J. Dreyer in Darmstadt 


Bei der Harmonischen Analyse mit Instrumenten (z.B. dem Harmonischen Analysato!' 
_Mader-OÖtt) und erst recht mit numerischen Schemaverfahren (z.B. nach Runge oder 
-- .  Zippereroder Hußmann) bereitet es’ beträchtliche Schwierigkeiten, Kurven mit hohen + 
‚Spitzen zu bearbeiten, wie sie etwa in der. Technik bei Stoßvorgängen oder in der Medizin bei 

- Elektrokardiogrammen vorkommen. Die einem sonst vielleicht recht gleichmäßigen Kurven- 

_ verlauf überlagerten Spitzen bedingen hohe Harmonische, deren instrumentelle Bestimmung 

sehr zeitraubend ist, während die üblichen Schemaverfahren, wenn man nicht gerade ungewöhn- 

lich viele Ordinaten benutzt, sogar falsche Ergebnisse auch der niedrigen Harmonischen liefern 

können. Die-Schwierigkeiten sind um so größer, je ausgeprägter und j Je schmaler die Spitze ist. 

Hat man solche Analysen in größerer Zahl auszuführen, so ist es zweckmäßig, die Spitze 

_ durch ein Dreieck anzunähern, den dreieckförmigen Verlauf von dem übrigen Kurvenzug durch 
Subtraktion abzuspalten (Bild 1) und dann das Dreieck getrennt formelmäßig zu behandeln. . 

Die Spitze braucht dazu nicht streng dreieckförmig zu sein; die Abweichungen werden in der 

rue mitbearbeitet. 


Er ——— B ” - Y 3; 72 ;? 2 2 Baar & eyes rn Rex 6 
.. Die Fourierkoeffizienten c, — Ja?+ 5 der Dreieckkurve las en sich für 
0. Dreiecke sehr bequem bestimmen, wenn man ein für allemal die Kurve er 


2 1—eos® ec 
ga)= — 


— 


FA 2? 


aß: 


Bild ı > Bild 2 


genügend genau aufträgt (Bild 2). Ist nämlich die Dreiöckbreite 2b, die Höhe h und die Perioden- % 
n 3 
länge 7, so braucht man nur an den Stellen „—n"f=n'27,; die Ordinaten g„ abzu- 


T | 
lesen; sie ergeben nach Multiplikation mit dem konstanten Faktor A fdie Amplituden der Harmo- 


nischen 
— cosnf 


| 
„= an hf n2 f 
Der Phasenwinkel ist für alle Harmonischen 90°, wenn die Spitze am Periodenende,liegt. Dazu 


d h 
kommt ein weiterer Winkeln a’ = n' T 360°, wenn die Spitze um d nach links verschoben ist. 


Auch für den seltenen -Fall, daß eine zweckmäßige Annäherung der Spitze nur mit einem 
ungleichschenkligen Dreieck zu erreichen ist, lassen sich entsprechende Hilfskurven zeichnen 
und vorteilhaft anwenden. 


D. Mathematische Maschinen und Instrumente 


Vom ABACUS zum ENIAC : 
Von Heinz Wittke in Clausthal 


Das erste wissenschaftliche Rechengerät ist im Orient (Indien, China) erfunden; es gleicht 
dem Zählgerät unserer Kinder. Später haben es auch die Römer und die Russen übernommen. 


Jene nennen es „ABACUS“, diese „SCHTSCHOT“. Noch heute findet. man das Gerät in den 


russischen Amtsstuben und Lohnbüros. 

Als Neper (1614) die Logarithmen entdeckt, rechnet man mit diesen; Multiplikationen 
und Divisionen werden als Additionen und Subtraktionen gelöst. 

Pascal entwickelt 1642 die erste Addiermaschine, Leibniz 1672 die erste Multi- 
plikationsmaschine. Damit ist der Weg gewiesen. Andere Pioniere folgen: Polenus in 
Paduv (1709) mit seinem Sprossenrad; Babbage (1833) mit seiner unvollendeten ana- 
Iytischen Maschine; Wiberg (1863) mit seiner schreibenden Differenzenmaschine, mit der 
er seine im Jahre 1876 in Stockholm erschienene Logarithmentafel berechnet; Hollerith 
(1889) mit seinen Lochkartenmaschinen; Boll&e mit dem Einmaleinskörper (1889); Ha- 
mann (1905) mit seinen Proportionalhebeln’und (1925) mit seinem Schaltklinkenrad; Ellis 
(1906) mit seiner schreibenden Buchungsmaschine, die heute zur „‚National-accounting-machine‘“ 
weiterentwickelt ist. 

Aus der Vielzahl der heute vorliegenden Maschinen seien nur einige typische Beispiele 
herausgegriffen, die für wissenschaftliche Berechnungen wichtig sind: 

'Brunsviga20. Sie hat eine hohe Kapazität 11Xx12x20, Rückübertragung, teilbares 
Hauptzählwerk und ist für hochstellige wissenschaftliche Rechnungen hervorragend geeignet. 


Al Ze Da u Na A a Th rn Fe ae Ad u Da Ze > De Zee 


3 K paztät. mare Keine  Rückübertragung. 10 | 


BR = BR Sx8x18 Automatische Anika und Division, Kraei Yuliplikai 
ER Be ae ne 


ER "schine, ng ekoppelt eignet = ee im ee a Öi 
GUT hochstellige uitRetatische Berechnüngen. er reed 0 En a oder . ala RS 
Be ER EN BAR 
se © >:.Mauser- order. Universal: Den ersten Se Kasaen haben 
"die Mauserwerke. (Obefndort a am N geben, Dem ee don Schreibwerks eebr der, ‚hohe 
_Preis- entgegen. | Ba 
SZERemington- R an ei: N ew Yor "k: "Der Berpkten: Band Printing-Calewiator er 2 
- die erste schreibende. Kleinrechenmaschine, die ‚addiert, subtrahiert, multipliziert. und dividiert. 
- Vollautomatische Division und Multiplikation mit Druck der Eingangs- und Aa 
Ta Kleine ‚Kapazität 5x3x10. Keine Vorzeichenschaltung bei ,Multiplikation. ee Br 
N National- acco unting- machine, London: Buchungsmaschine für Addi- 2 
ER tion und Subtraktion.- ee für er N Untertaflung, numerische Inte- TE 
Sees gration usw. geeignet. SSR : Re. ar 
Fa Se Vs ollerithm aschine: Zahleneingang durch Lochkarten. Drückwerk ochkartenst 
RR speicher. ‚Anwendung für Differenzenrechnung, Untertaflung usw. In Verbindung, mit dem 
Rechenlocher lassen sich auch Multiplikations- und Divisionsaufgaben lösen, z. B. Koordinaten- 
aan ntegen ! in der Geodäsie, wobei die, Ergebnisse fortlaufend gedruckt werden. 22 


N ae Kuse»Hechenamlafe: ‘Sie beruht äuf dem Dualsystem. Einfachste Bauelemente. aR 
Steuerung der Rechenoperationen. Die Zuse-Anlage ist. eine „‚Klein“-Anlage In Vergleich zu 
den nachfolgend besprochenen amerikanischen und englischen Großanlagen. Die a, Ber. 
ist noch nicht abgeschlossen. R 


Elektronische Bechenänlsgen rc md zu nennen: a N 

_ > ASCC, Automatic Sequence Controlled Calculator (Aiken). Harvard University USA. > 

ä « ACE, Automatic Computing Engine; British National Physical Laboratories in Cambridge. 

EDSAC, "Electronic Delayed Sequence Automatic Computor; aneahes! aber) der? 

* University of Cambridge. BR 

Der z.Z. im Bau befindliche Rechenautomat der British National Physical Laboratories in BR 

Teddington bei London. A 

ENIAGC, Electronic Numerical Integrator and Computor; University of Dee in Phila- ms 

‚delphia, Moore School of electrical Engineering. EN 
EDVAC, Electronic Discrete Variable a MESRERENIORRUNE des ENIAC, ehöntalls >“ 

‚University of Pennsilvania. N 

R IBM-SSEGC; Selective Sequence Electronic Eaiesleron New York a ae 

MANIAG, gebaut‘ von Neumann, Princeton University. Be 

Bell Telephone Laboratory machines: Sie sind entwickelt von Stibite und Williams in 
den Bell Telephone Laboratories, New York. 


Diese elektronischen Rechenanlagen haben eigene Tafelwerke, eine hohe Speicherkapazität, 
steuern den Rechenablauf und haben eine ‚‚elektronische‘“ Rechengeschwindigkeit. So löst der 
IBM-SECC 3500 Additionen bzw. Subtraktionen von 19stelligen Zahlen in einer Sekunde. Die 

R. Anlagen sind für wissenschaftliche Rechnungen eigens erbaut und hervorragend geeignet. 

RI “ Nun noch ein kurzer Blick in die Zukunft. Wir benötigen eine kleine, vorzeichentreu 
rechnende, schreibende ‚‚Reiserechenmaschine‘“, billig. im Preis. Die Cordt-Universal und die 
Remington-Rand sind Ansätze dazu. Jedoch ist die erste zu teuer, und die letztere hat eine zu 
kleine Kapazität. Beide rechnen noch nicht algebraisch. Ferner benötigen wir eine wirtschaft- 
liche, mittelgroße Rechenanlage zu erschwinglichem Preis. Hier scheinen die von Walther 
und Zuse beschrittenen Wege Erfolg zu versprechen. Die Entwicklung ist noch nicht ab 
geschlossen. - 

In bezug auf die Rechenaufgaben können heute Probleme angefaßt werden, vor deren Lö- 
sung man bisher angesichts ihres ungeheuren Umfangs zurückgeschreckt ist. In der Geodäsie 
3* 
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‚36 D. Mathematische Maschinen und Instrumente * 


Man kann sich heute der Lösung kontinentaler und interkontinentaler Ausg) 


harren umfangreiche Tafelberechnungen ihrer Lösung, bedingt durch die, Umste 
360° auf 400°, von den Bessel- und Clarke -Erdellipsoiden zum 


zuwenden mit linearen Gleichungssystemen von 100, 500, 1000 Unbekannten. Die Erddimi 


sionen lassen sich durch Zusammenfassung geodätischer, geophysikalischer und astronomischer EM 24 


Messungen neu ableiten. In der Astronomie und bürgerlichen Zeitrechnung kann man dem 


Gedanken einer dezimalen Zeiteinteilung näher treten; der Einwand, daß viele | 
überlieferte Tabellen umzurechnen sind, wird hinfällig angesichts der uns nunmehr bescherten 
elektronischen Rechenanlagen, und die Instrumente wären in älnlicher Weise umzustellen, wie 
es beim dezimalen Längenmaß, dezimalen Gewichtsmaß und dezimalen Winkelmaß (400 g) erfolg- 


reich durchgeführt ist. 


Die mathematischen Voraussetzungen für die Entwicklung 
logistisch-kombinativer Rechenmaschinen 


Von K. Zuse in Hopferau . 


Nach einer kurzen Einführung über die bisherige Rechengeräteentwicklung von Herrn 
Zuse wird der Unterschied zwischen numerischen und den sog. ‚logistischen‘ Rechenmaschinen 
besprochen. r 

Während die bisher bekannten Rechenmaschinen der Zahlenrechnung dienen, ist das Ar- 
beitsgebiet der logistischen Rechenmaschinen die Lösung schematisch-kombinativer Aufgaben. 

Es läßt sich folgende Rangordnung in der Rechengeräteentwicklung aufstellen: 

1. Tischrechenmaschinen zur Lösung der arithmetischen Grundoperationen, - 
2. Programmgesteuerte automätische Rechengeräte mit Speicherwerken usw., 
3. Logistische Rechengeräte, 

ı 4. „„Künstliches Gehirn“. 

Bis zur Stufe 2 ist diese Entwicklung praktisch erprobt. Bei Stufe 3 ist der theoretische 
Ansatz klargelegt und die praktische Ausführung in absehbarer Zeit möglich. Stufe 4 ist vor- 


erst nicht aktuell. 


Die logistische Rechenmaschine arbeitet wie die numerischen Rechenautomaten der * 


Stufe 2 nach Anweisungen, die jedoch auf einem höheren mathematischen Niveau liegen. 

Der Unterschied der Arbeitsweise beider Gerätetypen wird an einem Beispiel klargelegt. 
(Aufgabe: Differentiation. Numerisches Rechengerät löst sie durch schrittweise Errech- 
nung der gesuchten Kurve. Die logistische Rechenmaschine ermittelt die Zeichenfolge 
der mathematischen Funktion der Ableitung.) x 

Die Anwendbarkeit des logistischen Gerätes erstreckt sich auf alle Operationen, welche 
sich in strenge Regeln fassen lassen. Auf mannigfachen Gebieten der Technik und der Wirt- 
schaft liegen derartige schematische Aufgaben vor. 

Die Entwicklung des logistischen Rechengerätes erfordert zunächst die Aufstellung eines 
allgemeinen rechnerischen Kalküls als Voraussetzung für die konstruktive Lösung. Die Grund- 
lagen hierzu liefert der Aussagen- und Prädikatenkalkül der Logistik. Jedoch müssen erhebliche 
Abänderungen und Erweiterungen mit diesen Kalkülen vorgenommen werden, um sie den 
Zwecken des praktischen Rechnens anzupassen; die Logistik hat zum Ziel, implizite Definitionen 
von mathematischen Sätzen aufzustellen. Rechengeräte hingegen benötigen explizite Anwei- 
sungen, welche für die in Frage kommenden Fälle explizit den Weg der Lösung angeben. Diese 
Forderung wird durch den von Herrn Zuse entwickelten „All gemeinen Plan- 
kalkül“ erfüllt. _ 

Im Allgemeinen Plankalkül wird der Begriff ‚‚Rechnen“ wie folgt definiert: „Rechnen 
heißt, aus gegebenen Angaben nach einer Vorschrift neue Angaben bilden“. Eine mathematische 
Analyse des Begriffs der „‚Angabe“ ergibt, daß sich grundsätzlich alle Angaben aüs Ja-Nein- 
Werten aufbauen lassen. Der Unterschied in Auffassung und Darstellung zwischen dem For- 
malismus der mathematischen Logik einerseits und dem Allgemeinen Plankalkül andererseits 
wird an dem Beispiel des Begriffs der ‚‚Kohärenz“ gezeigt. Der Allgemeine Plankalkül ist in 
den letzten Jahren von Zuse soweit durchgebildet worden, daß die Voraussetzung zur prakti- 
schen Entwicklung logistischer Rechengeräte geschaffen ist. 


\ 
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=: ‚3. Die mit Elektronenröhren arbeitende Relaistechnik. +» Ns 
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& B: ‚Da die weitere Entwicklung dieses Projektes in. mathematischer Hinsicht eine Se um- 
n _ fangreiche Aufgabe darstellt, wird eine Zusammenarbeit zwischen wissenschaftlichen und indu- 
E% ‚ striellen Ag EL Re \ 


IE ER Lösung gewöhnlicher Differentialgleichungen mit der EM 
EN FEN 8 Integrieranlage IPM-Ott A 
ee ER Von Alwin Walther in Darmstadt | | 


De Das Institut für Praktische "Mathematik (IPM) der Techn chen Fihchechere Darmstait 3 
a - hat in den letzten zehn Jahren zusammen mit der Firma A. Ott in Kempten (Allgäu) eine 
neue Integrieranlage zum instrumentellen Lösen gewöhnlicher Differentialgleichungen entwickelt. 
= ES "In dem zunächst geplanten Umfange ist sie jetzt fertiggestellt und erprobt; erste Anwendungen 
- haben ihre Brauchbarkeit. erwiesen. Wie ‚jede stetig arbeitende Integrieranlage (bekanntestes 
Be 2 - Beispiel: Differential Analyzer von Vannevar Bush am Massachusetts Institute of Tech- 2 
'nology; Gegensatz: ziffernmäßig arbeitende Rechenautomaten wie ENIAC mit finiter nume- BE 
_» . rischer Behandlung der Differentialgleichung) liefert sie eine partikuläre Lösung zu gegebenen 
| Anfangsbedingungen als gezeichnete Kurve oder als Wertefolge in einem Zählwerk. Um mög- 
. lichste Allgemeinheit zu erzielen, ist sie in einzelne Rechengeräte aufgelockert, die von Fall 
zu Fall gemäß der zu lösenden Differentialgleichung verkoppelt werden: Integriergeräte (deren 
bisher zwei vorhanden sind) lassen von der höchsten Ableitung zur Stammfunktion aufsteigen, 
Addiergeräte (bisher acht) und Multipliziergeräte (bisher eines) leisten die vier Grundrech-. 
nungsarten, Funktionstische (bisher drei, Zeichenfläche 620 x 460 mm) führen gegebene Funk- 
tionen mit Hilfe gezeichneter Kurven ein oder dienen wahlweise zum Aufzeichnen von Ergebnis- ERS 
kurven. Die Verkopplung geschieht mittels Steckschnüren an einer Kopplungstafel; dadurch B 2 % 
werden elektromagnetische Fernübertragungssysteme, die sich an den Ein- und Ausgängen der i 
 Rechengeräte befinden, zusammengeschaltet. Das implizite Gegebensein der Lösungsfunktion See 
durch die Differentialgleichung spiegelt sich in kreisartig geschlossener Verkopplung der Rechen- ° 
- geräte wieder; auf diesem von William Thomson!) 1876 klar formulierten „Rückkopplungs- | 
prinzip‘ beruht die Möglichkeit instrumenteller Behandlung von Differentialgleichungen. Die 
Anfangsbedingungen werden durch die Anfangseinstellungen der Rechengeräfe, die Maßstäbe 2 
durch Zahnrad-Wechselgetriebe berücksichtigt. 5 EN 
"Während V. Bush als Integrationselement die Planimetepoie (das Reibradgetriebe, den 2 


Gonellaschen Integrationsmechanismus) mit Drehmomentverstärkung verwendet, wirdim 
Integriergerät der Anlage IPM-Ott das im gewöhnlichen Integraphen bewährte, für Differential- / 
‚gleichungen erstmals durch Udo Knorr?) vorgeschlagene Schneidenrad benutzt. Stellt man | 
seine Richtung laufend gemäß der Ableitung y’ ein, so steuert es seine unter ihm mittels eines 
Unterwagens und eines Oberwagens in x- und y-Richtung bewegliche Unterlage nach der Funk- ° 

tion 9. Däbei wird es, um Rückwirkungsfreiheit herbeizuführen, durch eine lichtelektrische 
Nachlaufeinrichtung entlastet. Zum Addieren und Subtrahieren dient das übliche Differential- 
getriebe. Multipliziert wird nicht, wie sonst vielfach, mit‘zwei Integriergeräten auf Grund der u 
Teilintegrationsformel vv = f udv + a vdu, sondern mit einem besonderen Multiplizier- und 
Dividiergerät. Es stützt sich auf den elementargeometrischen Proportionensatz und ist ähnlich 

dem Integriergerät gebaut... Die Richtung eines Lineals wird laufend nach dem einen Faktor 


; 1) W.Thomson, Mechanical integration of-the linear differential equations of the second order with variable coeffi- 
eients. Proc.roy. Soc. London 924. (1876), 269—271; Mechanicalintegration of the general linear differentialequations ofany order 
with variable coefificients, 271—273; Addendum, 273—275 

2) Vgl. insbesondere U.Knorr, Apparat zur selbsttätigen Aufzeichnung des Fahrdiagramms,. Elektr. Kraftbetr.u. 
Bahnen 12 (1914), 310—313; Über einen Integraphen zur mechanischen Integrationeinersehrallgemeinen Gruppe von Differential- 
gleichungen. Diss. TH München 1921. 
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eingestellt, ein Unterwagen nach dem anderen Faktor verschoben; dann läßt s das Produkt RER 
‚als Verschiebung eines Oberwagens abgreifen. Klemmungen werden durch eine lichtelektrische "FRE u; 
Nachlaufeinrichtung verhütet. Beim Abgreifen einer Funktion auf dem Funktionstisch liegt 
ein besonderer Vorteil darin, daß die Funktionskurve nicht von Hand, sondern mit einer vn 
-H.-J. Dreyer entwickelten lichtelektrischen Abtasteinrichtung selbsttätig nachgefahreg wird. Y R 
Die ideale mathematische Kurve ist hierbei nicht die Mitte, sondern der eine Rand eines realen 
Striches; sie kann deshalb durch eine Photozelle nach dem Kriterium halb hellen, halb dunklen r 
Gesichtsfeldes verfolgt werden, und zwar mit einer Genauigkeit von +0,1mm. Im übrigen we | 
beträgt nach vielfältigen scharfen Genauigkeitsproben, deren Originalaufzeichnungen beim Vor- , 
trag. vorgezeigt wurden, der Fehler der einzelnen Rechengeräte weniger als 0,1 Prozent dr 
größten vorkommenden Ordinatenerstreckung, der Fehler der ganzen verkoppelten Anlage 
Bruchteile eines Prozents. Besonders empfindliche Proben lieferte die Lage von Nullstellen bei 
reinen oder gedämpften Schwingungen, bei der. Besselschen Funktion von der Ordnung 1/3 
‚und beiHermiteschen Polynomen 6. und 8. Ordnung. Jeder Durchlauf für eine Differentia- 


gleichung dauert etwa 15 bis 30 Minuten. 
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Die Anwendungsmöglichkeiten der neuen Integrieranlage können ‘hier nur durch einige 
Stichworte gekennzeichnet werden, etwa: Formelmäßig schwer angreifbare Differentialglei- 
chungen — Ersatz mühsamer numerischer Bearbeitung — Differentialgleichungen mit lediglich 
empirisch bekannten Daten, wofür als Beispiel eine Schwingung mit versuchsmäßig gegebener 
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mn Rückstellkraft gezeigt wurde. Besonders geeignet ist die Anlage für Serienuntersuchungen. 
f: Hierbei kann es sich darum handeln, im Sinne der,,Mosaikmethode“ einen Überblick zu gewinnen 
% oder Parametereinflüsse zu beurteilen. Vor allem aber gehören hierher Rand- und Eigenwert- 
3 probleme, bei denen‘ man in systematischem Probieren die Anfangsbedingungen so lange ab- 
E* wandelt, bis die Randbedingungen erfüllt sind. Als Beispiel für ein Randwertproblem, noch dazu 


mit Parametereinfluß, wurden die Aufzeichnungen für die Temperatur eines Glühfadens mit Ab- 
strahlung bei verschiedenen Umgebungstemperaturen vorgelegt ?). 
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Bericht über die Entwicklungsarbeiten an der Intesrieranlage 
der Rechenautomaten G.m.b.H. Göttingen 


Von Hans Büchner in Göttingen 


Es handelt sich um eine mit Reibradgetrieben ausgestattete moderne Integrieranlage. Die 
Verbindung der Elemente untereinander erfolgt fernelektrisch. Funktionskurven werden foto- 
elektrisch abgetastet. Besondere 10-stufige Übersetzungsgetriebe ermöglichen in Verbindung 
mit Summentrieben die Einstellung von über 600 Faktoren. Eine Hauptschalttafel ermög- 
licht die übersichtliche Herstellung topologischer Verbindungen, und zwar zugleich für drei 
Rechenschaltungen. Steigerung der Rechengeschwindigkeit auf durchschnittlich 5 Minuten für 
die Darstellung partikulärer Integrale. 
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Numerische Rechenmaschine mit Magnetophonspeicher ') 
Von H. Billing in Göttingen ?) 


ls wird ein Speicher für 192.20 stellige Dualzahlen beschrieben. Die Zahlen werden nach dem Magneto- 
phonprinzip auf einer schnell umlaufenden, mit Eisenpulver beklebten Trommel gespeichert. Pläne zum 
Bau einer numerischen Rechenmaschine unter Verwendung dieses Speichers werden erläutert. 


‚The report gives a description of a storage system for 192 numbers, each number having 20 digits in 
the binary system. Similar as in Ihe magnelophone numbers are stored on a cylinder rotating with high speed 
which is covered with iron powder. Desingns for the construction af a calculating machine are described where 
his storage system is applied. 


Le rapport donne la description d’une magasin pour 192 nombres, chaque nombre ayant 20 chiffres 
dans le systöme binaire. Comme dans le magnetophone les nombres sont emmagasines sur un cylindre tour- 
nant Ires vite qui est couverl de poudre de ferre. Des dessins pour la construction d’une machine & calculer 
sont deerits ou ce magasin est employe£. 


3) Vgl.auch A. Walther, Einige Anwend / Pü 
2. Okt. 1948,. Referat Arch. Math. 1 (1948), ats der Integrieranlage IPM-Ott. Vortrag Tübinger Mathematikertagung 


1) An Ste ‚* case A: e 
TE n Stelle des Vortragstreferates wurde hier ein am 26. 7.48 eingegangener Aufsatz des Verfassers 


°) Institut für Instrumentenkunde in Göttingen. ; y 


über das gleiche Thema 
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N“ ar Durchführung. ‚größerer, 'numerischer Rechnungen wie vor allem zur Duniartchen 
= sung von Differentialgleichungen sind seit einigen Jahren in ‚den USA und England Rechen- 

. maschinen in ‚Entwicklung, die unter Verwendung von Elektronenröhren die Rechnungen voll- 
ER automatisch mit großer Geschwindigkeit durchführen. Die erste derartige bereits seit Ende 1945 : R 
in Betrieb befindliche Maschine ist der amerikanische Rechenautomat ENIAC. Er ist im Prinzip 
aa eine Übertragung einer mechanischen Rechenmaschine etwa vom Typ Mercedes-Euclid auf Bau- 
- elemente der Hochfrequenztechnik. Er führt im wesentlichen die arithmetischen Grundope- 
Ser ‚ rationen aus. Um die hohe Rechengeschwindigkeit ausnützen zu können, ist er einmal miteiner 
. Kommandoanlage verbunden, die Art ünd Ablauf der einzelnen Rechenoperationen. befiehlt 

> und zweitens mit einem Zahlenspreicher, der Ausgangs- und Zwischenwerte zu ‘speichern ‚ge- NE = 
‚stattet. Obwohl die Maschine mit sehr großem Aufwand gebaut wurde — es wurden 18000. 

Radioröhren verwendet —, erwies sie sich für die Lösung allgemeinerer Aufgaben, vor allem 

hinsichtlich des Fassungsvermögens ihrer Zahlenspeicher, noch als zu beschränkt. ‘Es sind daher 

an verschiedenen Stellen Entwicklungen angelaufen, um den Aufwand der Maschine zu verklei- 
nern und das Fassungsvermögen des Zahlenspeichers wesentlich zu erweitern. Beides ist mög- Dr 
lich, wenn man von dem Vorbild der mechanischen Rechenmaschine weiter abrückt und de 
as bereits in den en mehr den elektro-technisch gegebenen MU En En 
ap t Te 
ER Ein wesentlicher. Schritt in dieser Fhtuge besteht darin, (daß man sämtliche Rechnungen ee 
3 ‘im Dosbsstem durchführt ?). Das Dualsystem enthält bekanntlich nur die beiden Ziffern 1 und 0. 
0. Dies begründet, daß man im ganzen Rechengang jeweils nur zwischen zwei Möglichkeiten zu 
„entscheiden hat; man kann die Rechnung aus einer Folge von ‚‚Entweder — Oder“- -Operationen 
‚ aufbauen. Da auch das wichtigste elektrische Schaltelement, der einfache Schalter, die zwei 
Möglichkeiten ‚‚offen‘‘ oder „geschlossen“ hat, leuchtet ein, daß das Dualsystem den elektro- 
technischen Gegebenheiten bestens angepaßt ist. Der Nachteil, daß eine Dualzahl mehr als 
‚dreimal soviel Stellen hat als eine gleich große Dezimalzahl, wird dadurch bei weitem aufgehoben. 
Entsprechend scheinen die modernen Entwicklungen elektronischer Rechenmaschinen sämtlich 
> das Dualsystem zu benutzen. 
= Der zweite Fortschritt besteht in der Entwicklung leistungsfähigerer Zahlenspeicher. An 
die Zahlenspeicher ist neben der eigentlichen Speicheraufgabe als zweite Forderung zu stellen, 
daß die Zahlen aus diesen binnen sehr kurzer Zeit entnehmbar sein müssen. Diese Entnahmezeit 
©, kann:sonst leicht das Rechentempo der ganzen Maschine bestimmen. Der ENIAC braucht zur 
- Darstellung einer Dezimalstelle 10 Radioröhren, die man sich nebeneinander angeordnet denken 
“möge. Ist z.B. nur die 7. Röhre stromführend und alle anderen gesperrt, so bedeutet das die 
- Ziffer 7. Zur Darstellung einer 10 stelligen Dezimalzahl werden mithin 100 Röhren benötigt. 
Beim Dualsystem braucht man zur Darstellung jeder Stelle nur eine Röhre. Da eine 33stellige 
. Dualzahl etwa einer 10stelligen Dezimalzahl entspricht, würde die Verwendung des Dual- : 
systems den Aufwand für den Speicher immerhin auf den 3. Teil herabsetzen. Da man für eine ie 
größere Bechnung einige Tausend Zahlenspeicher benötigt, bleibt der Aufwand trotzdem noch ua 
zu groß. | K 
Eine wesentliche Aufwandsverminderung bedeutet der in der‘ britischen Eitericklung in 
Cambridge verwendete ‚‚Schallwellenspeicher ®). Die Grundidee besteht darin, die Dualzahl 
"micht als eine Reihe statisch gleichzeitig nebeneinander bestehender Zustände (etwa.wie oben 
geöffnete oder gesperrte Radioröhren) darzustellen, sondern als eine Reihe zeitlich aufeinander- 
folgender Impulse. Die Zeit wird dabei in gleiche Intervalle A unterteilt. Wird in einem be- 
stimmten Intervall ein Impuls gegeben, so entspricht das der Ziffer 1, unterbleibt der Impuls, 
so entspricht das der Ziffer 0. Die zeitliche Aufeinänderfolge der Zeitintervalle entspricht der 
Aufeinanderfolge der Ziffern in den einzelnen Stellen einer Zahl. Zur Speicherung wird im „ 
Schallwellenspeicher die Zahl zunächst in eine Folge elektrischer Stromimpulse verwandelt. 
Mittels eines Schwingquarzes, der sich als Geber an einem Ende eines mit Quecksilber gefüllten 
Rohres befindet, wird die Impulsfolge in einen Schallwellenzug umgewandelt. Dieser läuft mit 
Schallgeschwindigkeit durch die Quecksilbersäule. Am Ende des Rohres befindet sich ein zweiter 
Quarzkristall, der Empfänger, der die ankommenden Schallimpulse wieder in elektrische Im- 
= pulse umwandelt. Diese werden unter Verstärkung und Herstellung der ursprünglichen Impuls- 
form unmittelbar wieder auf den Geber gegeben, womit sich das Spiel wiederholt. So’läuft die 
einmal eingegebene Zahl als Impulsfolge dauernd in dem aus Quecksilberröhre und Verstärker „ 


3) Eine im Dualsystem arbeitende Maschine wurdevin Deutschland bereits während des Krieges von Zuse gebaut. 
-. 4) Discovery, Februar :1948. 
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bestehenden Kreis um. In einer einzelnen ksilbe: e Können | 
gespeichert werden. Ein gewisser Nachteil derartiger dynamischer ner besteht 
satz zu den statischen Speichern darin, daß die Zahl nicht jederzeit zur Verfügung steht un« 
daß die Speicherung erlischt,.wenn der elektrische Strom abgeschaltet oder auch nur kurzzeitig 
unterbrochen ‘wird. Im britischen Schallwellenspeicher läuft die Zahl in wait einmal] 
um und kann mithin alle 10-? s entnommen werden. DR Be 2 
Noch kleineren Aufwand erfordert ein Dualzahlenspeicher, der nach dem Magneto 
prinzip arbeitet und im Institut für Instrumentenkunde in Göttingen vom Verfasser entwickelt 
wurde. Hierbei wird die Zahl ähnlich wie beim Magnetophon als Folge von magnetischen Stöo- 
rungen mittels eines Elektromagneten im Magnetpulver fixiert. Da die Zahl jedoch möglichst 
jederzeit aus dem Speicher entnehmbar sein muß, befindet sich das Magnetpulver nicht wie 
üblich auf einem aufspulbaren Magnetophonband. Vielmehr ist die Außenseite eines Aluminium- 
zylinders von 15cm Durchmesser mit dem Magnetpulver beklebt. Der Zylinder rotiert mit 
100 Umdrehungen/s um seine Achse. Seiner Oberfläche sind bis auf 0,2 mm nun 
chermagnete genähert, die im Abstand von 1 cm übereinander angeordnet sind. Die M 
ähneln in ihrer Bauart den üblichen Magnetophon-Gebemagneten. ‘Von einem als Geber verwen- 
deten Speichermagneten können rund um den Trommelumfang in der bisherigen Versuchsaus- 
führung 12 20stellige Dualzahlen = 240 Ziffern aufgeschrieben werden. Diese 12 auf dem glei- 
chen Umfang der Trommel liegenden. Speicher sollen im folgenden als die Umfangsspeicher be- 
zeichnet werden. Die einzelnen Ziffernimpulse haben voneinander einen Abstand von 1,5° 
Trommelumdrehung. Der gleiche Speichermagnet dient ebenfalls zur Entnahme der Zahl. Da 
bei Ausnützung einer Trommelhöhe von 16cm 16 Speichermagnete übereinander angebracht 
sind, können bisher 12 - 16 = 192 20stellige Dualzahlen gespeichert werden. Da die Trommel 
mit 100 Umdrehungen/s umläuft, stehen die Zahlen des Speichers nur jede 1/100s zur 
Verfügung. Dadurch wird die Arbeitsgeschwindigkeit einer mit diesen Speichern arbeitenden 
Rechenmaschine auf höchstens 100 Operationen/s (Addition, Multiplikation...) beschränkt. ° 
Unter Verwendung des obigen Magnetspeichers ist im Institut für Instrumentenkunde der 
Bau einer numerischen Rechenmaschine geplant. Die eigentliche Planung ist im wesentlichen 
abgeschlossen, eine große Reihe von Vorversuchen sind durchgeführt und die ersten Teilstücke 
sind bereits gebaut und in Erprobung. 
Bevor auf Einzelheiten der Planung eingegangen wird, soll an einem Beispiel gezeigt wer- 
den, welche Anforderungen an die Rechenmaschine zu stellen sind. Als Beispiel werde die Schrö- 


dinger- Gleichung 
y’+F&)y+T()=0 (1) 


gewählt. F(x) und T(z) sind numerisch vorgegebene Funktionen. Bei Annäherung durch ein 
Polynom dritten. Grades läßt sich die Schrödinger-Gleichung zur numerischen Berechnung um- 
formen in 
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+h?[—F (a) y(a)—10F (a+h) y(a+h)]+T(a)-+10 T(a+h)+T(a+2h)}+-Rest 


Zur Vorbereitung der Rechnung sind die Konstanten sowie die Funktionen F (x) und 7x) 
für die Argumente = a -+-N - h auf den Zahlenspeicher zu geben. Zu diesem Zweck werden 
die Funktionswerte zunächst maschinell in 20stellige Dualzahlen umgewandelt. Diese Dual- 
zahlen werden, wie es vom Hollerithsystem bekannt ist, in einen Lochstreifen eingeprägt, der in i 
jeder Reihe Platz für 40 Löcher hat. Die ersten 20 Löcher dienen zur Aufnahme der Dualzahlen, 
die nächsten 8 Löcher zur Wahl der verschiedenen Magnetophonspeicher, in die die Zahlen ein- - 5 
gegeben werden sollen, der Rest zur Kommandierung von Rechenoperationen. Das Lochband 
läuft über eine Kontaktwalze und ist mittels einer Zahnradübersetzung mit der Speichertrommel 
derart gekuppelt, daß bei jeder vollen Umdrehung der Speichertrommel eine neue Lochreihe ab- 
getastet werden kann. Im Lochband ist die Zahl wie in einem statischen Speicher durch 20 
nebeneinanderliegende Elemente (die Löcher!) gegeben. Der Eingeber muß diese statisch vor- 
liegende Zahl in eine zeitliche Impulsfolge verwandeln. Da die Ziffernimpulse auf der Speicher- 
trommel einen Abstand von 1,5° Trommelumdrehung haben sollen, ist der Eingeber mit der 
Trommel fest zu kuppeln auf eine derartige Weise, daß das Zeitelement At für eine Zifterlänge 
stets genau einer Trommeldrehung von 1,5° entspricht. Auf diese Weise ist das richtige Arbeiten 
des Eingebers unabhängig von der Drehzahl der Trommel. 

Da die Ziffernstelle in der Zahl durch den Ort definiert ist, an welchem der entsprechende 
Magnetimpuls auf der Speichertrommel liegt, muß einmal die Lage der einzelnen Speicher auf 
der Trommel genau definiert sein. Zweitens ist dafür zu sorgen, daß die Eingabe der Zahlin den 
gewählten Speicher genau in dem Moment beginnt, wenn der Anfang des betreffenden Speichers 


(2). 
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Nachdem die Ausgangswerte in die en eingebracht sind, a die ee 
E se gemäß Gleichung (2) beginnen. Die Addition zweier Dualzahlen, die als zeitliche Im- 
"pul folgen ‘vorgegeben sind, geschieht verhältnismäßig einfach in einer Anordnung von Radio- 
-  röhren. In diese ‚Anordnung laufen, mit der kleinsten Stelle beginnend, die beiden Impulsfolgen 
‚gleichzeitig hinein. Die Röhrenanordnung hat nichts weiter zu tun, als die einzelnen Ziffern 
"nacheinander von hinten her beginnend zu addieren, wie man es bei schriftlichem Addieren ge- 
"wöhnt ist. Sie muß die Operationen 0-+0,0-+1, 1-1 ausführen. Das Resultat erscheint 
Sr ‚Stelle für ‚Stelle gleichzeitig mit dem Einführen der beiden Summanden am Ausgang. Größere 
Schwierigkeiten als die eigentliche Addition macht es, die beiden Summanden gleichzeitig zur 
Addition zur Verfügung zu stellen. Zweiin verschiedenen Umfangsspeichern gespeicherte Zahlen 
- laufen ja nicht gleichzeitig unter den Speichermagneten hindurch. Um die beiden Zahlen trotz- 
_ dem gleichzeitig zur Verfügung zu haben, werden sie auf zwei Summandenwerke gegeben. Ein 
- Summandenwerk besteht aus einem Gebemagneten @ (Bild 1), der ebenso wie ein Speicher- Be 
_ magnet die Zahl in dem Magnetpulver der Speichertrommel « ein- 
BIS prägt. Auf dem gleichen Trommelumfang, im Abstand von genau 
einer Speicherlänge = 30°, befindet sich ein zweiter als Empfänger 
 wirkender Magnet E, der eine Speicherlänge nach. der Eingabe die 
 Ziffernimpulse der Zahl abnimmt und sie nach Verstärkung und 
-__  Impulsformung im Verstärker V wieder auf den Geber gibt. So 
läuft die Zahlin dem durch Geber, Empfänger und Verstärker ge- 
bildeten Kreis dauernd um — ähnlich, wie es oben beim Schall- 
ES wellenspeicher geschildert wurde — und steht somit beim Beginn BildT 
eines jeden Umfangsspeichers zur Verfügung. Der Löschmagnet L Be 
R löscht die eingeprägten Magnetimpulse, so daß sie nicht nach einer vollen Umdrehung der Trans e 
melin den Summandenspeicher nochmals gelangen können. NL 
-  » Die Multiplikation zweier Dualzahlen geschieht wie bei der schriftlichen Multiplikation 
RR durch wiederholte Addition. Um die Rechengenauigkeit nicht herabzusetzen, hat das Multi- ” 
.  plikationswerk die doppelte Stellenzahl des Additionswerkes. .Die notwendige Stellenverschie- 
- bung des Multiplikanden um je eine Stelle pro Addition läßt sich leicht bewerkstelligen, indem 
- man den Multiplikanden in einen Umlaufsspeicher bringt, der dem obigen Summandenspeicher 
ähnelt. Der "Abstand zwischen Geber. und Empfänger‘ beträgt jedoch eine Ziffernlänge 
mehr als zwei volle Speicherlängen. Die eigentliche Addition dauert nur eine Speicherlänge 


= ,: 10-3 5, die ‘Multiplikation dagegen zweier 20stelliger Dualzahlen beansprucht die 


A0fache Zeit = !/,, s. Dazu kommt noch je %/oos für die Herausnahme der zu multipi; 
zierenden Zahlen aus dem Zahlenspeicher. 

2 Während die Subtraktion durch eine kleine Zusatzeinrichtung vom Additionswerk mit 
erledigt werden kann, soll auf den Einbau eines Divisionswerkes bei dieser ersten Entwicklung 
verzichtet werden. Man kann einmal den Rechengang meist so wählen, daß man Divisionen 
weitgehend vermeidet. Andererseits läßt sich die Division durch eine iterative Methode zur Er- 
mittlung der Reziproken ersetzen, ohne eine Division auszuführen. 

Das eigentliche. Kommandowerk besteht aus einem Lochbandgeber, wie er schon oben 
beim Zahleneingeben erwähnt wurde. Mit seiner Hilfe werden die. Speicher ausgewählt und 
Befehle gegeben, welche Operationen mit den einzelnen Speicherinhalten auszuführen sind. 
Das Lochband gibt alle !/,os einen Befehl. Zu jedem Befehl wird eine neue Lochreihe ab- 
getastet. Das Lochband für eine z.B. 5 Minuten währende Rechnung hätte demnach 30 000 , 
Lochreihen. Damit würde es nicht nur unangenehm lang, seine Herstellung wäre auch äußerst 
zeitraubend. Die meisten größeren mathemätischen Rechnungen führen jedoch irgendwie auf 
Iterationen oder sich immer wiederholende gleichartige Rechenschritte, wofür Gleichung (2) 
ein Beispiel gibt. In dıesem Fall ist es anzustreben, daß das Kommandoband lediglich die Be- 
fehle für einen einzelnen Rechenschritt zu enthalten braucht. Durch eine zusätzliche Schalt- 
einrichtung müsseu von Rechenschritt zu Rechenschritt die nötigen kleineren Umschaltungen 
‘vorgenommen werden, so daß ein Rechenschritt nach dem anderen durchgeführt wird. Das 
Kommandoband wird dazu zu einem Ring zusammengeklebt. Jedem vollen Umlauf des Kom- 
‘mandobandes entspricht dann je ein Rechenschritt. Die für derartige Aufgaben vorgesehene 


Pr a" 


e- Zusnehrerhuun erfüllt im ee die Aufgae, ne‘ W: 
 Kommandoband in Wirklichkeit der Speicher N + M geöffnet wird. M muß sic 
 Rechenschritt um 1 erhöhen. Nur ein von der jeweiligen mathematischen Aufgabe a 
0 Teil der Speicher wird über diese ‚zusätzliche Einrichtung ko kommandjert, der a 

* kommt seine Befehle weiterhin direkt vom Lochband. 


bedürftig ist. Vor allem muß die Zahl der Speicher wesentlich erhöht werden. Unt 


' Ein wangeihnliches Gerät für haryonneh Analyse und Synthese = 


Ein 


Es ist uns klar, daß die geplante, von uns nur als Modell Aiigeschene Mäschine 1 


von 2000 Speichern für 30stellige Dualzahlen auf einer Trommel erscheint ohne üb ormäßigen 
Aufwand als möglich. Da zweifelsohne für die verschiedenen mathematischen Aufgaben noch 
weitere Zusatzeinrichtungen benötigt werden, veröffentlichen wir bereits jetzt unsere Fe % 
der ER aus Kreisen der angewandten BAR ARRONR Ei ee = £ Di 


Von A. Walz in Teningen (Baden) ar FE Br ER 
Die üblichen. harmonischen as ee; werten die Fourier- IR 


ye „470 rar TE 2: a) j 2 
21 eines EEE 3 5. (2) | 
ö # S , Fi 

1 | re | 
= 5, En er (3) # 
M) g 


durch eine Planimetrierung aus. Es ist dabei erforderlich, dem graphischen Bild der zu nr # 
sierenden Funktion f(x) mit einem Fahrstift nachzufahren. = 

Das hier zu beschreibende Gerät unterscheidet sich von diesen Analysatoren grundsätzlich ' 
dadurch, daß dieFourier (1), (2), (3) mit Hilfe der Besse1lschen Summenformeln - 


2n—1 £ i j 
a,= e SR eHs (v5) She rer (4) E| 
1 Bam = | 
a = SE) SERRIUER  FERR E (5) 
ee v=1l,3,.n—1 
= 5, De (6) 
=U 
1 we 
=; Def) Ten e  e (7) 
e=V0 


ausgewertet werden. Es wird also bei diesem Gerät nicht der ganze Verlauf der zu analysierenden 
Funktion für die Berechnung der Fourier -Koeffizienten verwendet, sondern nur eine be- 
schränkte Anzahl äquidistanter Funktionswerte fo, fi fa: ++ -fo- 

Die Anzahl dieser ausgewählten Funktionswerte pro Periode sei 2n, Dann kann das Glei- 
chungssystem (4), (5), (6), (7) n BE -Koeffizienten, n— 1 sin-Koeffizienten, sowie das Glied a, (den 


„ Flächeninhalt der Funktion f(x) von = 0 bis 2) liefern, 


Das Ergebnis der-Be.ss€ schen Summenformeln stimmt dann exakt mit dem Ergebnis 
der Fourier-Integrale (1), (2), (3) überein, wenn’sich die Funktion f(x) nur aus Harmonischen 
bis zur n-ten Ordnung zusammensetzt. Enthält f(x) auch noch Harmonische höherer als n-ter 
Ordnung, so kann man mit den Koeffizienten aus den Besselschen Summenformeln zwar 
immer noch eine endliche trigonometrische Reihe 


Fix) = 4,4 4,008 2 + A,cos2r +..... + A„cosn« 
+ 2, sin + B,sin2®..... + B,.,sin(n— 1) & 


!ı) Ein ausführlicher Bericht wird im Archiv der Mathematik erscheinen. 


ae ee 


En: aufbauen, die an den alien "Stellen ve a der. vorgegehähen Funktion ft (x) = A a, 
608 c4.. =D, a ‘übereinstimmt; an den Zwischenstellen ergeben sich jedoch 
Pe Abweichungen zwischen F(x) und f (x) und die Koeffizienten A,, B, stimmen nicht mehr genau 
mit den wahren Fourie r-Koeffizienten U; 0 überein, vor allem bei den höheren Ordnungszahlen. 
EN Dies muß als Nachteil der Besselschen Formeln und damit auch unseres Gerätes ge- 
wertet werden. In der- Regel hat man es jedoch i in der‘Hand, sich einer gegebenen Aufgaben- 
“ stellung durch Wahl einer hinreichend großen Zahl von Funktionsstellen 2n anzupassen. Über- 
dies interessieren oft gar nicht die wahren Fourier- -Koeffizienten a,, by, sondern nur die 
‚.Besselschen Koeffizienten A, B, der endlichen trigono- 
Red) 
' metrischen Reihe, die eine nur durch eine diskrete Zahl von 1.3 
- Werten gegebene Funktion interpoliert. ' | 
00. Das neue Gerät bildet auf ‚rein mechänichem Wege 
die Produkte f(x) cos(vx) und f(x) sin (rx) und addiert 
diese Produkte nach der Vorschriftder Bessel schen Formeln. 
Es wird dabei von einem sehr einfachen mechanischen Prinzip 
"Gebrauch gemacht, nämlich dem statischen Moment, das ein 
Hebel von der Länge f(x) mit Einheitsgewicht an seiner 
Drehachse erzeugt, wenn & um den Winkel vx gegenüber 
. der Lotrichtung ausgelenkt wird (Bild 1). Werden 2n sol- . 
cher Hebel auf einer gemeinsamen Achse angeordnet, deren 
Länge den Funktionswerten fu fr» Ja ----- fo proportional 
ist und deren Auslenkung gleich dem Winkel £—»x bzw. 
vx ist, so ist das an dieser Achse angreifende resultierende 
‚Drehmoment bis auf einen ‚Proportionalitätsfaktor unmittel- 
bar gleich den Koeffizienten A, bzw. B,. Um dieses Dreh- 
moment bequem und genau messen zu können, wird die 
Drehachse an beiden Enden in Schneiden gelagert.. Diese Bild ı 
Achse mit den darauf fixierten Hebeln bildet dann das 
schwingende System einer Waage, die z. B. mit einem Sehe an einem ron 
‚abgeglichen werden kann. = 
Zur konstruktiven Ausbildung des Gerätes ist im.einzelnen noch Be zu bemerken. 


Für eine bestimmte Ordnungszahl n kommen nur ganz bestimmte Winkelauslenkungen 
der Hebel vor (o = 05 00,1, 2.20 — ı) . Um das Einstellen der einzelnen Hebel in 
die vonder Besselschen Formel geforderte Winkelauslenkung und das Festhalten des Hebels 
in dieser Stellung zu erleichtern, ist jeder Hebel mit einer halbkreisförmigen- Scheibe fest ver- 
bunden, in der die Winkelteilung durch äquidistante Nuten oder Langlöcher festgelegt ist. Durch 
diese Löcher oder Nuten wird, nachdem die Hebel in die richtige Auslenkung gebracht sind, ein 
abgeflachter Bolzen hindurchgesteckt, der in zwei fest mit der Drehachse verbundenen Stützen 
gelagert ist. Durch diesen Bolzen ‚werden dann die Hebel in ihrer Auslenkung gegenüber der 
Drehachse festgehalten. Die endgültige Justierung erhalten die Hebel durch eine Drehung dieses 
Haltebolzens um 90 Grad. Erst nach dieser Drehung steht der Bolzen in Berührung mit beiden 
Nuten- oder Loch-Flanken. 

Bei großer Zahl:von Hebeln wäre das Herstellen der Auslenkung der Hebel ebenso um- 
ständlich wie das Nachfahren der zu analysierenden Funktionskurve mit einem Fahrstift. Da die 
Auslenkung der Hebel für eine bestimmte Harmonische unabhängig von der speziellen Aufgabe 
stets die gleiche bleibt, kann man für jede Harmonische eine treppenförmig gestufte Schablone 
vorsehen, mit der die gleichzeitige und irrtumsfreie Auslenkung aller Hebel vollzogen werden 
kann. Bemerkenswert ist, daß diese Schablonen keine Präzisionsarbeit erfordern, sie brauchen 
die einzelnen Hebel nur in die ungefähr richtige Auslenkung zu bringen. Die genaue Justierung 
erfolgt, wie schön erwähnt, durch den Haltebolzen. 

Das Gerät hat die angenehme Eigenschaft, ohne weiteres auch für die Durchführung von 
harmonischen Synthesen verwendbar zu sein. An Stelle der Funktionswerte werden in diesem 
Fall an den Hebeln die Fourier-Koeffizienten eingestellt. Die anzuwendenden Hebelauslen- 
kungen sind die gleichen wie bei der Analyse. Es können daher auch die gleichen Schablonen 
benutzt werden. Jede Wägung liefert die Summe der cos- bzw. sin-Glieder der trigonometrischen 
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' gleichmäßigkeiten in der Massenverteilung der einzelnen Hebel ergeben bei jeder Ha mo ischen 
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Durch die Verwendung des Meßprinzips einer Schneidenwaage ist mit ve hi Itnism 
geringer Präzisionsarbeit eine drei- bis vierstellige Genauigkeit der Ergebnisse erre 


eihe. Durch Addition dieser beiden Summen hat man ( 


Pal) 


- 


nur eine feste additive Korrektur. ? a EBEN 
; Das Konstruktionsprinzip des Gerätes erlaubt eine Weiterentwicklung zur vollautoma- 
tischen Durchführung von harmonischen Analysen und Synthesen. 3 
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HAUPTAU FSÄTZE 


Konvergenztypen und das Verhalten von Approximationen 
in der Nähe einer mehrfachen Wurzel einer Gleichung 


Von E. Bodewig in Den Haag, Holland az 


| 
' 

Zunächst wird ein Maß für die Konvergenz einer Folge aufgestellt, und der „Grad der Konvergenz‘ 
definiert. Damit ist es dann len zwei Lücken in der Theorie der Approzimationsmethoden zu schließen. 4 
Dies wird angewandt auf die Approximation von Wurzeln einer Gleichung, und es wird ges, daß die Kon- \ 
vergenz fast jeder Methode abhängt von der Vielfachheit-der zu approzimierenden Nullstelle. Z. B. gilt 
dies für die Newtonsohe Methode und ähnliche Methoden, für die Laguerresche Methode 
und für eine allgemeine Methode zur Konvergenz n-ten Grades. i 

First a measure of the convergence of a sequence is established and a new notion, the "degree of con- 
wergence‘‘ is defined. By this it is possible to shut two gaps in the theory of all approzimative methods. This 
is applied to the approximation of roots of an equation, and it is proved that the convergence of nearly all t 
methods will depend on the multiplieity of the root in question. This holds f.i. fr Newton’s method and‘ 
for similar methods, fr Laguerre’s meihod and for a general method of a convergence of degree n. nn 

D’abord' une mesure de la convergence d’une suite est #tablie et le « degr& de convergence » est defini. 
Par cette motion il est possible de remplir deux lacunes de la theorie de toutes les möthodes d’ zimations 
suttessives. Cela est appligue & l!’approximation des racines d’une Equation, et il est prouve que la convergence 
de presque toute methode depend de la multipliplicite de la racine & approximer. O’est valable p. ex. pour la ' 
methode de Newton et des möthodes semblables, pour la möthode de Laguerre et pour une 
generale qui fournit une convergence du degre n. 

CHayana YCTAHAB-IUBAeTCH Mepa MIA CXONHMOCTH HEKOTOPOH MOCHEANOBATENIBHOCTH | 
H ONpepengeTca ‘,,cTeNeHb CXOAHMOCTH“. HTO AACT BOSMOSKHOCTB BANOAHHTB 184 mpodesa } 
B_ TeOPUH &ANIPOKCHManna. JBaTeM 9TO HPHMEeHseTcH IPH IIPHÖNMKEeHHOM ONpeneneHun 1 
KOpHeH YypaBHeHHA, IIPHYeM HOKASEIBACTCH, YTO CXONHMOCTB HOYTH BCEX METONOB ANLPORCH- 
MaNHH 3AaBHCHT OT KPATHOCTH KOPHA. ITO MeÄCTBHTENBHO, HANp., AAIA MeTona HusowToHa u 


IA AHAJIOTHYHBIX METONOB, MIA MeToNa Jlareppa H MIA OÖINeTO MEeTOoNA, KOTOPHH NacT CXOAH- 
MOCTB n-HOA CTeIeHn. 4 


Ein Hauptproblem der praktischen Analysis ist das Aufsuchen von Nullstellen von Funk- 
tionen. Dazu benutzt man sukzessive Approximationen, wobei man von einem Schätzungswert 
der Nullstelle ausgeht und diesen sukzessive verbessert. Es gibt nun Verfahren, wo man von 
einem beliebigen Schätzungswert ausgehen kann und stets die nächstbenachbarte Nullstelle 
als Grenzwert einer Folge von Approximationen bekommt. Das sind aber nur ganz wenige 
Verfahren, wie z.B. die Laguerre’sche Methode, auf die wir weiter unten zu sprechen 
kommen’ werden. Sonst ist es fast stets so, daß der erste Schätzungswert gewisse Bedingungen 
erfüllen muß, damit die Folge der Approximationswerte konvergiert, und zwar gegen die benach- 
barte Nullstelle konvergiert. Wie man einen brauchbaren ersten Schätzungswert bekommt, 
ist eine. Sache für. sich und hat nichts mit der Approximationsmethode selber zu tun. 

Nun bestehen aber in der Theorie aller Approximationsmethoden 
zwei Lücken. Man beweist nämlich erstens nur, daß die Methode konvergiert, aber 
es fehlt die Möglichkeit, verschiedene Methoden zur Lösung desselben Problems miteinander 
hinsichtlich ihrer Wirksamkeit zu vergleichen. Man hat sich bisher vielmehr damit begnügt, 
an einem oder mehreren Beispielen die Verfahren zu demonstrieren, und zog aus diesen speziellen 
Parken. Beispielen Rückschlüsse auf die Methode im allgemeinen. Dieser Schluß ist natür- 
ich falsch. 

Zweitens nahm man an, daß die Konvergenz einundderselben Methode unter allen 
Umständen ähnlich verlaufe, soweit man überhaupt instinktiv von einem „ähnlichen Verlauf“ 
sprechen konnte. Denn ob die Konvergenz z. B. abhing von der Vielfachheit der zu approxi- 
mierenden Nullstelle, wurde nicht einmal als Problem gestellt. Man nahm vielmehr stillschwei- 
gend an, daß dies nicht der Fall sei, was: wiederum falsch war. 


h BY Ni SEE 
ER TRERER E ER (n—X) | 
5 (a kann GEBE irgerideine positive. Zahl sein. ) Führt man die Abveichingen der, Glieder ı von. on 
hrem EL ein: In — - = An, so, wird. 


=e+0. für. a ag = 
- girr. "% EL: : 


k En 


BE Re m I, d I ee ee 
TE RE EN ee UN, 00% ' 
ur Br ER 7 . N j ® Zr N E j Fe Bern 7 g” PR Es ! EN N 
Slip: den Grad der Konvergenz zu Bean hat man. , t dus 1 nach De von m © 
zu ‚entwickeln. 'Beginnt diese Entwicklung mit dem Me ch. so ist die ee vom. 
9, und der zugehörige Koeffizient ist, c. 


- Je größer gist, um so schneller konvergiert die F olge. Und) inan kann sagen, daß, wen n. 
e Me Anzahl der genauen Dezimalstellen von ı, gleich k ist, das 
nächste Glied derReiheetwakggenaue Dezimalen hat. 
. Ein Beispiel für zwei konvergente Folgen verschiedenen Grades werden wir unten geben. 
Ich halte aus den eingangs erwähnten Gründen diesen Angriff des a > für 
sah a in der ganzen numerischen Analysis. | ae 
- Wir geben zunächst eine allgemeines Beispiel. Be : 3: a 


Die Reziproke eines linearen Operators. Wir suchen einen "Prozeß, welcher die 
Er Reziproke eines linearen Ber durch eine im- ‚Grade 'k konvergierende Babe sukzessiver 
Er ee un liefert. 

Sn Dazu suchen wir zunächst einen <Äteken Prozeß für die BR einer Zahl a; Dann. 
ST 2 und gemäß unserer Forderung muß | 


dnzı = = cd en 


©. Die einfachste Annahme ist offenbar 
er 3 | & re 


also | RAR x 
Ent — 1ja= c(&n — 1ja)*. 


Die Konstante c ist noch verfügbar. Sie muß so, gewählt werden, daß alle Potenzen des un- 
- bekannten 1/a sich auf beiden Seiten herausheben. Dies ist der Fall für c=at4, Wi 
Damit wird ae 
BE = En Enkel) ante) 2ı .. N u WO Zn AXn» SE 


_ Für k=2 bekommt man die einfache Formel 


ER I VL) a ER N Ge Sek ig 
. Die Folge x,, &,, %p ... ist also quadratisch konvergent, falls sie überhaupt-konvergiert. Leiz- 
teres ist aber,. wie man leicht sieht, der Fall, wenn Be 
S Bern Aldi, & ee 


Dieselbe Formel (1) gilt aber auch für jeden linearen Operator, da bei ihrer Ableitung kapde 
‚anderen Gesetze vorausgesetzt wurden als sie für jeden solchen Operator gültig sind. Die Kon- 
vergenzbedingungen selber müssen allerdings für jeden Operator gesondert festgestellt werden. 


ö 1) E. Bodewig, Sur la möthode de Laguerre. Akademie van Wetenschappen, Amsterdam, Proceedings, 
vol. 49, 1946, p. 911—921. ‘ 
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„endlich‘‘ ist. Daher ist dort mit großer Annäherung: F 
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Die Formel gilt also z. B. auch, wenn «eine Matrix ist. Dann h 
lim = at, Für Matrizen wurde diese Formel zuerst aufgestellt durch G 
und wurde seitdem mehrfach von neuem entdeckt. In einer früheren Arbeit ®) a 
Ableitung der allgemeinen Formel gegeben, die etwas anders verläuft Sl obige. re 


Sr, pre 
‚ıiall 


II. Newtonsche Formel und Verallgemeinerungen ER ER: A A 
Wir geben zunächst eine Verallgemeinerung der Newtönschen Formel, zugleich mit einer 


neuen Ableitung der gewöhnlichen Formel. Während man diese nämlich bisher mittels 47 


Taylorschen Satzes bekam, gehen wir hier anders vor. x 3% 
f(x) habe die p-fache Wurzel X: | S= Ve ‘ 
fa)=(e—Xp-g(e), wo 1g(M#+0 ...:....( 
l 
“= fa) = ple — Xp g+ (erg), a Be 
d.h. | B 3 
Sal) =Pla— + ll) nennen =. 2) 


In der Nähe von x = X ist das erste Glied der rechten Seite sehr groß, während das zwei 


ala, | 2 
X=2—pfe)lf'(e). 


Damit haben wirdie allgemeine Newtonsche Formel in der Nähe einer | 


p-fachen Nullstelle: ; 
5 Mn In+1= an -Ppflan)/f (2) ee Are a ee Se a TE Fe (3). B 


en a a ee a ee ER 20 
und 'bestimmen den Charakter der Konvergenz der Folge: 

TR räkees de Ar - 
Dazu haben wir &.11— X = d„ıı nach Potenzen von 2, — X —=d, zu entwickeln. Die Ent- 
wicklung von g(x) an der Stelle X liefert nach (1): - 

Fan) = (in — AP I (DH + (an -APgR) +... : 1 
= AR) +EUR) +..., | 


also 


Wir setzen aber allgemein 


‚ k \ 
\ Y h ! ar } . Pr L gi 
2 } & . fi 
sn rn nn a nn nn 


also IE } 
Fan paTg KH +E+NBR) +... 7 

Somit nach (3): R 
Er a, ag(X) j 
Hm ut —d +... Ei 

| n-+1 n ) » N Hp g(X) Pr ' 

also.nach Subtraktion von X auf beiden Seiten: e 
a a g(X) 3 

d, - (1-2)a — ——du—.... 3 

in nen A 

Das lineare Glied von d, fällt also nur- dann fort, wenn @a= p. Damit haben wir das et 
Ergebnis. ? 


1. In der Umgebung einer p-fachen Wurzel liefert die Formel (3’) eine nach der Wurzel 
quadratisch konvergierende Folge (fall$ die Folge überhaupt konvergiert), wenn «= p. + 

2. Hingegen ist für a=# p die Folge nur linear konvergent, und zwar umso schlechter, je kleiner Y 
das Verhältnis @«/pist. Hier unterscheidet sich die Newtonsche Formel in ihrer Konvergenz j 
also nicht wesentlich von z. B. der Regula Falsi, die ebenfalls nur linear konvergiert und ' 
oft einfacher anzuwenden ist, vor allem bei wiederholter Iteration,, S 

3. Die gewöhnliche Newtonsche Formel, d.h. diejenige mit «= 1 konvergiert daher nur in °_ 
der Nähe einer einfachen Wurzel quadratisch, in der Nähe einer mehrfachen hingegen 
nur linear, 


4. Weiß man, daß die Wurzel mehrfach ist, kennt aber ihre genaue Multiplizität nicht, so 
setze man a=2. 


Eine stets quadratisch konvergente Formel. Man kann aber auch 
eine quadratisch konvergente Formel erhalten, ohne die Multiplizität zu kennen. Dazu hat 


2) G. Schulz, 
p. 57—59. 


®) E. Bodewig, Bericht über die verschiedenen Methoden zur Lösung eines Systems linearer Gleichungen mit 
reellen Koeffizienten. Akademie van Wetenschappen, Amsterdam, Proceedings vol. 50 und 51, 1947 und 1948, — Wieder 
abgedruckt in Indigationes mathematicae, Amsterdam, vol. 9 und 10, 1947 und 1948, . e 


Iterative Berechnung einer reziproken Matrix, Z. angew. Math. Mech., Bd.*13, 1933, 


“ 
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i f(x) en It, aber nur ein ich. 
mit. a=l liefert für sie Ss ER, a ee 


wo a der u Seite, es &n der einfacheren Schreibweise halber wenielassen. R: 
wurde. we zu rechnen hatı man, wenn man in in der: ‚Form schreibt; BEE a 


EL 


: ee = m. 
Br un dr Charakter der Konvergenz allgemeiner zu Fe rsuchen, Se wir wieder 5) ; 
x aut der rechten Seite ae ‚Botenzen, von .d . Wir bekommen dann für Re 4, 

SS 37 5 % S | | 3 “ A er Sr nz fon) f( (&n) = = D- (p “ & = Ei - v8 a j A » } ® E%7 f: 2 
a Ger > k VE a ar 
EEE NT pi ae 2 (gg? Be ne 
PR wo 2 g, = für: ‚das: Krane Cr gebildet sind. Damit ei man aus (A: ER Er Rn 
= GE a ES x Ba "dn+ı = = .cdh+. u ; e=—g( Xlpg(X)+ 0. ar ee se 5 ee 


Se Auch für eine Einfache ‚Wurzel konvergiert das Velen also "quadratisch, allerdings wegen 
| “des Nennners p hier umgekehrt - weniger gut als für ‚eine ‚ mehrfache Wurzel. 


Eine ähnliche N wie we Is BASE NER LE I ER A ei a 
Hier bekommt man für f( a ei x ar . = EN N > ; 
ME RS rg: An = cd nt 02d at. SE 
0 WO 

Zu PNlp+a—ap) | = —latap— 2) Ba, 


Er verschwindet nun wegen : seiner Symmetrie i in a und» in den Fällen: a—1 und p haichig) 
sowie a beliebig und p = 1. Der erste Fall liefert Formel (4). Im zweiten Falle kann man noch 
= 6, zum Verschwinden bringen er ul 1% Somit haben wir 
En. Ergebnis. AHN 

Test 1. Formel (4) liefert eine stets quadratische Konvagah, Die Koreidee ist im Falle: einer . 
einfachen Wurzel am schlechtesten, und: wird besser, wenn die Vielfachheit der Wurzel 3 
wächst. E 

- 2. Formel (5) ist im Falle einer einfachen Wurzel bei heliebiesm a quadratisch konvergent, ist 
also nicht ‚besser als die Re Newtonsche Formel, die keine zweite lung: be- 


nötigt. > - £ 

3. Formel (5) gibt für a = 12 i im Falle einer einfachen Wurzel eine ku bi ische Kon-. SE 
> DE, ON mean a 
ef vergenz mt = — gs im Falle einer mehrfachen Wurzel jedoch 35} 
f g 3 2 g' rt 9 (X) gr 18 | a 


- nur eine lineare Konvergenz mit c, = (p — 1)/(p+1). : 

4. Eine stets kubische Konvergenz bekommt man, wenn. man von den Resultaten der For- 
meln (3) und (4) jedesmal den Mittelwert nimmt. Denn die Koeffizienten c sind in beiden 
Formeln entgegengesetzt gleich. 


Die vereinfachte Newtonsche Methode. Bekanntlich hat man an 
Stelle der gewöhnlichen Newtonschen Formel 


In+1 = un — fl (F- (Xn) 
oft eine andere verwendet, indem man auf der rechten Seite in dem Bruch nur den empfind- 
licheren Teil f(x„) stets wieder berechnete, hingegen den unempfindlicheren Teil f’(&,) konstant 
hielt, also ihn etwa durch f’(x,) erstzte. Die so entstehende Formel 


ET, | i mn nf) rennen (6) 

wird von vielen als beinahe ebenso wirkungsvoll. angesehen wie die ursprüngliche Formel, so 

z.B. von Whittaker-Robinson‘d). Dies ist aber offenbar nicht richtig. Denn wenn 

. die Folge &.4+1 = %n— cf (&„) überhaupt konvergiert — was wir hier, wo essich nur um den 

"Charakter der Konvergenz handelt, annehmen —, so ist natürlich dnzı = dn — cf (kn) 

— dn—cdhg (X) —..., wenn wieder f(x) bei = X eine p-fache Wurzel hat. Somit haben wir 
Ergebnis. Formel (6) konvergiert nur linear. 


2) Whi ttaker-Robinson; "Caleulus of observations, p. 91, 1926. 


ei er RN ” 
0. Das gleiche Ergebnis hatte schon Ost owsi?) beissen 
wesentlich einfacher. in muß hinzugefügt werden, os 2 
vergenzbedingungen angibt. a er RE 
Beispiel. Wir wollen den Unterschied eines quadratisch und 2 
vergenten Verfahrens und gleichzeitig ebensolcher Folgen an einem Zahlenbeispie! Sc 
berechnen wir eine Quadratwurzel einmal mittels der Newtonschen Formel, sod: 
mita = 1/2. Dann haben wir f (x) = # —b. Dann ist nach Newton 


= 


a 


y Ba ae: U in 
2 R- © % i En}1 = (zn + b/xn)/2, mit.o, = 1/2 yb a 
Br,‘ während die andere Formel TE ER Dr 
Br: N Er 41 = (ai + 2b)2,/2 +6) mi seljab i 
E ir gibt. Für d= 2 und von ,=1 ausgehend bekommt man die Folgen: 
A - quadratische Konverg.. kubische Konverg. 
en; | Be \ el s 
D- "5 Burr, i Xp 1.5 , 1,4 £ 
Be.‘ Eye Sl 1,4142132 - 
Be 2: 1,414216 . 1,414213562373095048796 de 
Bi, BR Kr 1,414213562375 RR 
Sa . Beide Verfahren konvergieren hier sogar noch, wenn man als Anfangswert 2, = 10 wählt. 
Pe Dann wird nach Newton 2,= 1,444, während die kubische Formel 9 genaue Dezimalen 
er liefert x, — 1,414213562. | | ag: 
FR ; II. Die Laguerre’sche Methode für algebraische Gleichungen und Verwandtes 
e Be >. Während die. Newtonsche Methode und die Methoden, die man als Verallgemeinerungen 
er von ihr auffassen kann, für beliebige Gleichungen und für reelle und komplexe Wurzeln gelten, 
Bi hat Laguerre®) eine elegante Methode angegeben, die nur für algebraische Gleichungen 
E: und nur für solche mit reellen Wurzeln gilt.. Ich habe diese Methode ausführlich an anderer 
2 N Stelle analysiert”), worauf ich hier nur: verweisen kann. - DER 
RE ER Die Folge ist definiert durch ! „ir 
| | i er nf) 


ER a ra » 

Fa) EVA) - 
wo n der Grad der Gleichung f(x)= O ist und H (x) = (n — 1)? f?— rn (n — 1)ff”' = Hesse’sche 
Kovariante von f. Es gilt sowohl das + wie das — Zeichen bei der Quadratwurzel. Die 
Methode approximiert nämlich zwei Wurzeln gleichzeitig. 

: Die Hauptmerkmale der Methode sind: E 
Be, 1. Wie Laguerre bewies, konvergiert die Methode stets, also bei beliebigem Ausgangs- 
wert x, (im Gegensatz zur Newtonschen Methode, wo x, gewisse Bedingungen erfüllen muß). 

2. Sie liefert, wie Laguerre ebenfalls bewies, für beide Vorzeichen der Quadratwurzel je \ 


eine konvergente Folge: R 
BE EN ah 
1} 
Re Wlan e D 


wobei X’ der zu x, nächstgelegene linke Wurzelwert von f(x)= 0 ist und X” der zu x, nächst- 
gelegene rechte Wurzelwert ist. — Dieser Satz gilt immer, wenn man sich die gerade Linie im 
Unendlichen geschlossen denkt. 


3. Wie ich dann bewies, ist die Konvergenz im Falle einer eirfachen Wurzel kubisch, im Falle 
einer mehrfachen Wurzel hingegen nur linear. 

4. Um im Falle einer p-fachen Wurzel ebenfalls eine kubische Konvergenz zu bekommen, hat 
man, wie ich bewies, in obiger Formel 4 (x;) zu ersetzen durch aH (z;), woa— (rn — p)/p(n—1) 
ist. 

9. Im Anschluß daran bewies dann van der Corput®), daß, falls die Gleichung wirklich 
mindestens eine Wurzel der Vielfachheit > p besitzt, die in (4) erwähnte Formel nur die 
beiden zu x, links und rechts benachbarten Wurzeln der Vielfachheit > p approximiert, 

x BR m w j = ‚ Über eine Modifikation des Newtonschen Näherungsverfahrens. Travaux de Y’Institut mathem, 


‘) Laguerre, Sur une methode pour obtenir par approximation les racines d’une &quation alg&bri i 
> h ue 
toutes ses racines reelles. Nouvelles Annales Math., 2e serie, XIX, 1880. — Oeuvres de ee 1898, Se D.87 Ios). 
x ’) E. Bo da ewi 8, Sur la meöthode de Laguerre pour l’approximation des racines de certaines @quations alge- 
briques et sur la critique d’Hermite. Akademie van Wetenschappen Amsterdam, Proceedings vol.49 , 1946, p. 911—921). 
B) Van d er Corput, Sur l’approximation de Laguerre des racines d’une equation qui a toutes ses racines 
r6elles. Akademie van Wetenschappen, Amsterdam,, Proceedings, vol. 49, 1946, p. 922-929. e j 


aı gen wichtigen Beleblelen nr Unterschied der | 
und einer AueBrEneneHL Nullstelle ERET Wir wollen un- 


wir an, wir a N Prozeß F @) wie früher, welcher zu u einem 
besseren Wert Br liefert: ERSTER RE En are $ 


Ei) + == en Bu = = en We er . ; . a {* 
ne a &%y, ... ansich und mit a Bedingungen, unter hen sie 


beschäftigen wir uns nicht, nehmen. vielmehr an, daß diese vorhanden ist. Wir such 
enige allgemeine Funktion a! welche eine im n-ten. konvergierende. For Er 


vn 


dr: = et. 
ee ee 
| Bedeutet Di= dijde x die i-te Ableitung, SO-isE Be ae GE EEE 
0 DFa)= Dian=ern—1).. ee 
Setzt ı ‚man hier = = X, ‚so verschwindet die rechte Seite, wenn i<.n. Somit ist. 


er BER URS x DF(X)=0 für wol 28 RE i 
er ist die Bedingung dafür, daß die Folge &,, &,... im Grade n konvergiert. Bevor wir daraus.” 
: AR (2) een wir.nöch Austern. . z 


. Fi ver 


m x Br F (lim >= F(X), 
ats = F (X). Diese Bedingung gilt für jede Wurzel X von f(«)=0. Daher ist ee 
Ex ee Enz +@)ge) BE a ee 
Sr wo ae lebt beliebig ist. Er 
IE R Soll aber (3) gelten fürn =2,d.h; ‚soll die Folse der % quadratisch konvergieren, s so 
SEE  muß=wegen 
RUE = | EOSSURK) IHR): Ref (X), 
SITES 98) =— 1 ()+FE)91 le), 
und die Funktion (4); | | es 
BEN re ER ETLDIELEFHR) EL INT as rel 
—- liefert für beliebiges g, (x) eine quadratisch konvergierende Folge der xy » = 
Sall jetzt (3) gelten für n = 3, d.h. soll die Folge kubisch sein, so muß 
0= DF(X)= —De(fif')+ De(fegn). 
Der Ausdruck f?g, wird nun nach der Leibnizschen Regel differenziert und im Resultat wird _ % 
xz= X gesetzt. Dabei verschwinden alle Summanden bis auf den letzten: g,:DX(fP) . 
| o—g,(X)- 2f’*(X). Andererseits ist an def Stelle X wegen f(X)=0: _ en 
Be 2.0. Desif)= Di(fr)=2Df: Dr+ Def. sfr =—fr | 
- wo zur übersichtlicheren Schreibweise r = 1 f’ er wurde. 
Somit ist 


ee 0= DRK) fr). 
Daraus folgt 2 = 
| (= — zer la)+fte) 2) = 


also 


[ 


® Fa=2—r Pfr 4rg, ee. (6). 5 


" Dieser Prozeß liefert, falls er konyergieri; eine zubach Ru eüte Folge 2, 2. 0 zbeirsbe- 
_ liebigem 95 ‘ x). 


" Um A eine im Grade n konyergente Folge zu bekommen, 
und (6) mittels des Operators 


o0= =r' JEM wo r = 1/f'(«) AR D = d/dx | 5 an ar = > 


zu: | 
> — f"(&) 9n(®) 1 EN ER EIT 
ER | Be In 
ul) FH er 
Nun ist wegen D(fr) = DD ae Or: v2 | 


DF,@)=1—- (+ On+(f Orts: RE TE Or+ff-Or)+.. 
re! Un VE an) Ir m Or. 
Aus Dr [f(X)m = m! (f(X))” folgt dann: 
„D"Fn(X) = D'1DF,(X) = Re FAT: BEER 


Andererseits ist 


DL (X) (X = nl’ (Krm(R), 
also 
| Hl = on. 


Somit hat F(x), wenn seine n ersten Ableitungen an der Stelle X verschwinden sollen, d.h., 


wenn die Folge x,, &, ... im Graden +1 konvergieren soll, die Form von Fy-1 (x). Die | 
"Form (7) gilt also allgemein, da sie fürn=1 gilt. 


Nach ihrer ganzen Konstruktion konvergiert die durch (7) gegliederte Folge der z,, & 
%a,... im Grade n (falls sie überhaupt konvergiert). Die hier gegebene Ableitung hat aber einen. 
gewissen symbolischen Charakter. Wie sich die Sache genau verhält, also in der Nähe 
einer Wurzel der Vielfachheit p muß noch besonders untersucht werden. 

Nun ist die vorstehende Ableitung nur gültig, wenn 


0=DF(N=1+f(X)g(X) 
gemacht werden kann, d.h. wenn f(X) #0, also X eine einfache Wurzel ist. Andernfalls wird 


eben die rechte Seite der vorigen Gleichung gleich 1. Nur bei einer einfachen 


Wurzel ist also die Konvergenz vom Grade n, bei einer mehr- 
fachen Wurzel ist sie wegen DF(X)#0O nur linear. 
Schließlich formen wir Formel (7) um und geben eine andere Begründung für sie. In 


seiner Differentialrechnung hat Euler °) folgende Verallgemeinerung des Newtonschen Ver- 


fahrens gegeben !P). 
Es soll von der Funktion y = f(x) derjenige Wert.X gesucht werden, der y= 0 macht. 


"Wäre die Umkehrfunktion & = u (y) also bekannt, so wäre die Lösung X gegeben durch X°= (0). 


Da dies aber nicht der Fall ist, vielmehr nur ein Näherungswert a vön X mit dem zugehörigen 
Funktionswert f(a) bekannt ist, so daß demnach X=a+h=u(0), so entwickelt Euler 
den Wert «4 (0) um die Stelle A = f(a) in eine Taylorsche Reihe: 


X = u(d—4) = ud) Au (A)+ 5 Aw (A) AA)... a . 


Die Ableitungen von u an der Stelle A sind aber bekannt, denn es ist 
uA)=a, w(A)=deldy=1:dylde= l/f(a)=r(a) 
u (A) = —fz(a)rr(a), W’(A)= [IfHa)—fı(a)fs(a)]rs(a)  usw., 
wo zur Abkürzung f, (x) an Stelle von fW (x Si gesetzt wurde. — Setzt man allgemein 
| u (y) = d" zldyr = (— 1)" 1 U, (a) rer (x), 
so ergibt sich für U,„ das rekurrente Gesetz 
On+ı= (An—1)f, Un—f'- dU,ldx. 


%) E u le N a Institutiones Caleuli Differentialis aL, Cap. ix. I Ziet Opera Omnia Ser. 15 vol. X, — 
2) Ich habe darüber berichtet in den Commentarii, Mathematiei Helvetiei vol, 8 1935, 36, pP. 1-4 und einige An- 
wendungen davon gegeben. & ; 


/ 


I ann Ba rn eh ae nn ten nun en nn tn nn 


des 
.gentenverfahrens 
gentenverfahrens mit anschließender 

.  Iteration (TVT/), die Frage offen geblieben, wie 
- sich diese drei Verfahren, die sämtlich Näherungsver- 
fahren zweiter Ordnung sind, hinsichtlich ihrer Ge- 


= 0017 - = On - — f Dun =pr: 


ormel (3% 
' Somit haben wir das Erge bnis. 


Om rk 
en [U nen = 
Ned. RE EN 


ie Neihenentwiektung von Euler Fr also identisch RT 


2 2 Die Euler sche Entwicklung ist eh, mit Formel 
* = Formel (7) liefert eine im n-ten Grade konvergierende folge: % Ba a RER 


denn PR S 7 


a “gs 


die Wurzel X, gegen welche die Felge konvergiert, einfach ist, ee eine lineare Falaea 


. wenn X mehrfach ist. 


Pe 


3. Bei einer mehrfachen Wurzel wendet man also Besser Harn! (4) von Seite 17 an. 
-4. Das Vorliegen einer mehrfachen Wurzel darf man a wenn f’ (21) an ist. 


» Eingegangen: 29.5.48. 


es KLEINE MITTEILUNGEN 


eranekalksrrafeleh: für die Halbschnitt- 


_ verfahren der graphischen Integration. PR 


‚In einer früheren Mitteilung!) „Bemerkungen zur 


graphischen Integration‘ ist bei der Erörterung der 
‚bekannten Halbschnittverfahren der ne Inte- 
einer BSR eEUE 2: 


=f(e, Y),. 


Sehnenverfahrens (SPY), des Tan- 
(TV) und des Tan- 


nauigkeit zueinander verhalten. Diese Frage soll im 


nn beantwortet werden. 


.a. OÖ. sind bereits die bei den einzelnen Markahten 


Se Zuwächse %j, kr und kıır mit den Anfängen ° 


ihrer Reihenentwicklungen nach Potenzen der Schnitt- 


länge h’ angegeben worden: 


sv: kr hf + S, %r En 
HdR +EDNH 
>: kır en Harz +h, yo +hfo 


= ; h2 h3 
hf +5 Dh +7 Plot Se 
TVI: ku=gfotgStoth, Yo + kım) 


; ’ h2 8 . 
Tr (D’f+pDf)+ --- 


- 


= RL: | 
(Dabei bedeutet D den Operator 5, +f öy: und der 
- Index o bezieht sich auf die Ausgangsstelle x4, Yo). 


Von diesen Entwicklungen aus erhalten wir durch 
Subtraktion des wahren Zuwachses 


’ h2 Ro, 
Ley, taNnNtgN F-- 
i en N 

=hfp+tg Dh + PIthDMot--- 


1) Diese Zeitschrift; Bd.20 (1940), S. 121—123. 


- 


die Fehler 


x -u-ierl- a re =; 


Be : ap 
dr -nok-m mar-4 6 Bor)ar me. ”. 


dur = Br Sie (5 er 5 Dr), + M(...) 


Die Frage ist jetzt, wie sich 61, örr und örrr zuein- 


ander verhalten. Wir beantworten sie durch Unter- 
suchung des sich beim Grenzübergang ho 0 Rabe 
‚den Verhältnisses 


| Bajörl: [6x] : |dmml]) = I(D*f + 4, Dal 
(2 Da Dyl:|(@ DPF +21, DNN 
A) Dort, wo ist, ergibt sich Er 
i 1 
la! : [rel : öl) = 5::1:1. 


Das heißt: dort ist das SY doppelt so genau wie die 
beiden anderen Verfahren, und diese sind untereinander 
gleich genau. — Bei der Oma dratur, d.h. bei 
der Differentialgleichung :y’ f(z), ist dies wegen 
f» = 0 identisch der Fall, so , daß bei ihr das SV in 


jedem Falle das genauere Halbschnittverfahren ist 
Ein TVI gibt es in diesem - 


(vgl. Zahlenbeispiel 1). 
Falle nicht. 


B) Ist [525 #0, | so bilden wir den Kennwert 


ae Df fr t+2ffry+S? Behr 
BIT DE EtL). 


Dann ist 
lim (|ö1}: |6tzl 
h—0 


und.es gibt folgende sechs Möglichkeiten: Es ist (im 


: öl) =I|Ro+41:]2R, —4|: 2 Ko +2), 


‚Limesfür k — 0) für 


ee |örl < lörrıl < lOr! 


— 2 <K,< 0 lörır! < dr | < I$inl 
0 <Ko< 05 [örml <|örr |< lörlı 
05 <K, <2 [öır| <Idıml < löil 
2 <K,<8 lörl <löil.  < lörml 
8 <K, <<» Iörl <Törıl: <lönıl 


4* 


wenn 


ir 


fahren gebietweise verschieden, End es ist in dem 
Gebiet, in dem 


K <—2oder K > 8ist, das SP 
£ VRR <R ist, das 7V 
2<K<05 ist, das TVI 
Halbschnittverfahren. 
Zahlenbeispiele (auf Grund numerischer 
Berechnung der Zuwächse): 
1. Zu A (Quadratur): 


1 
y- Ir mit der geschlossenen Lösung 


y=arctgx + const, 
o=1h,h=0,2 


SV: kı =h 


das genaueste 


a) 


1. 0,090498 
EUEe 
1+ (20 + 2) 


h 1 
Ma eier aseıoee ln 
7 
k=arctg(@,+h)—aretge,=aretgl2 7 
= 0,090660 
dr = — 0,000162; örr = 0,000322 


1 a y2 
y’ = xy mit: der geschlossenen Lösung y = yo Kan e 
Hierist K(z,y) = ? 2 und daher stets 0<KX < 2. 
Also ist 
a) im Gebiet 2> K, na > 0,5, 
d.h.für- je <V3: Be <lörıl <lörl, 
b) im Gebiet 0,5> RK, = 


d.h. für 
In Übereinstimmung damit: ist 
a) für, =1, h= 0,2: 
K,=1,k= 049216, 
kı = 0,48400, ör = — 0,00816 
kıı = 0,48800, örr = — 0.00416 
kyrı= 0,50000, örm= 0,00784 
b) für 20 = 5, h = 0,02: 
Ko = 0,077, k = 0,21078 
kı = 0,21042, ör = — 0,00036 
kır = 0,21044, ör = — 0,0003 4 
kııı = 0,21099, örm = 0,00021 


H. Heinrich. 
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a 


Kujbyschew. 


Halbraum mit halbkugelförmiger Schub- 
belastung'). 
I. Aufgabenstellung 

Drückt eine Kugel auf eine Platte gleichen Werk- 
stoffes, die hier als Halbraum aufzufassen ist, wobei 
die Elastizitätsgrenze nieht überschritten werden soll, 
so stellt sich nach Hertz [1]? in der Druckfläche 
eine halbkugelförmige Druckverteilung ein. Der Ra- 
dius A der Druckfläche bestimmt sich zu 


A Van, Yo ar N 


Hure Für die Anregung zu dieser Arbeit möchte ich meinem hoch- 
verehrten Lehrer, Herın 0, Prof. Pr. phil.L.Föpp1, meinen ver- 
bindlichsten Dank aussprechen. 

°) Die Zahlen in eckigen Klammern beziehen sich auf den 
Literaturnachweis am Ende der Arbeit. 


Dr, Bari da Fr ori ie 


lzl > V3:] örml = |örıl <lörl 


as mit a 


Halbkugel über 
ga=- ty Ron mit @=uPp ++ Q) 


Die Schubkräfte g sind alle in Richtung des Schlüpfens 
2 er die aus der Schubbelastung 
resultierende ue in der Ober- 
fläche des Halbraume, In er Umgebung der relativ 


a a ee 
ungen wie für die Platte, nur mit um- 
E oe rzeichen. Der Halbraum mit Hertz-. 
untersucht; die Lösung dient in Verbind Weber[2] 


I. Lösungsdurchführung 

Unter Benutzung eines rechtwinkeli z, 9, 2 
Koordinatensystems, dessen x, y-Ebene Peosakenne 
ebene des Halbraumes ist, lautet für eine einzelne 
'Schubkraft@ nach Bild 1 der Spannungszustand [3] in 
der Oberfläche in Zylinderkoordinaten mit = r : c08g, 
y=r-sin® 

(0:0, [3 cos®p 
—(1—2»)eosp (3cos?9— 2)] 

(0y),=0=5 a le cos p sin? @] (4). 


(Fry)z= 0= 3 3 cos? p sin o 


— (1— 2») sing (3 cos? —1)] 


ar? 


Bild ı Bild 2 


In Zukunft wird auf die Zeiger x = 0 verzichtet, da 
es sich immer um die Spannung .n in der Oberfläche des 
Halbraumes handelt. Zunächst stellen wir uns eine 
beliebige Schubbelastung q vor, die palallel zur X- 
Richtung und über eine Kreisfläche vom Radius A 

°) Der „ligemeine Ansatz, der vom Schlupfgrenzfall stetig zu 


völligem Haften führt, liegt ineinernoch unveröffentlichten Arbeit 
des Verfassers vor, 


En mr et ee PL 


nsgrenzen für y wurden ‚noch nicht 
en unterschiedlich sind, und 
y=n, und für = 


a 


Bild 3 


FAR Ka 


“ In den Grenzen für oin (6) bedeutet % die halbe 
“ Sehnenlänge. 2 


PN ‚Um die unbequemen Grenzen für 0 zu vereinfachen, 
‚setzen wir (vgl. Bild 3) 


3 m e=rcosy— sc08Y, de=ssiny-dy, 
: damit wird 
E % 3 
3 a q:‘ssiny : 
| ET Ye Por: dy B). 


y=0 


Für die Hertzsche Schubkraftverteilung gilt, wie 
- man aus‘(3) und Bild 3 abliest, 


R Io et A Fu 7 
eg | BEE (9) 
und aus (9) folgt $ 
ag sony BI, (10)., 


Wenn man g aus 1s (10); in (8) einsetzt und cos y durch 7 
_ substituiert, folgtmitdn/dy= —siny 


+1 
2 = % 
ee Dart zii ; 


ne 


yı n? 


5 in: Do,ay A) 
3 cosy —n 


{5° vu = —afesin (an 
2 
la sie für Punkte innerhalb und außer- 


y=- — are sin (dl, und ee sin 


" und -bietet keine Schwierigkeiten. An Stelle. der 


De UBR 


\ di I Inner: 
: ze SL. «dor, Kreisfläch, 


RD \2\ En aner SIR3 
ee, (eon'zı 21 der . ‚Rreisfläche, Er 


= Einsetzen von Do, = f(p = y) nach (8) in B 


bietet die Integration über y keine Schwierigkeiten. 
Es folgt für 0% innerhalb. der ea : 


mtr, B-u-en] dir rzala) 


Da für r > A die Tntegrationsgtenzen für M höchstens 


8 erreichen, ist durch den "Faktor a 
Vorzeichenwechsel während: der. Integration zu er- 
ray so daß sie in einem Schritt anseelahl werden 3 

ann, 


Es folgt mit Mi I =& ‚außorhalb der " Bolastungsfläch 


co% 


Ur c0sP en are sine 


re 
Fee Ze a-2n|Karc 


® He | (for) Fur. 
Die’ Berechnung von o, und 7,, geht analog vor sich. 


Funktion 9, =f(p+y) ist Po, bzw. De ein- 


zusetzen, eh: 

Dog, = 6vcos (p-+ y)sin? (p En 

D,,, = 3. (p+psin@+Y) 2 65 
—(1— 2») sin (P + y) (308? (p #V)- 1) 

- bedeutet. 

Die , liefert mit der Ars e= A/rfür- 


el =" 34 " 608p rc I “ . (15a) 
rei gy=g, ee | x 
1. (15h), © 2 
ee 
r <A Ty= go esinp 7; + 2v)] “>. (16a) = 
BEL ER 
rZ A Ty= q8inY "17, aTesine 
EM 1 3 
ee 
+/1 &(& (200% 2) r) 
(16b) 


4e 
+ zele(20p—;) + a) 


III. Ergebnis: 


Der Spannungsverlauf nach Bild 5 erinnert stark 
an das ebene Problem von Rad und Schiene [5]. 


— (1-2 - an arcsine 


Di, 


Sale “ 


Außerhalb der Belastungsfläche, d.h. für r>A4, 
klingen alle Spannnngen schnell ab. Bild 5 zeigt den 
Spannungsverlauf von o, längs der x-Achse. Während 
beim ebenen Problem o, unabhängig ’von » ist, schwankt 
hier entsprechend »=0 bis r = 0,5 die maximale 
Randspannung 0,0 von 1,57, bis 1,7799. Da mit 


Bild 5 


r-cosp=.x nach (13a) o, und nach (15a) o, für 
r < A unabhängig von y sind, geben o, und o, in der 
Berührungsfläche ein . hufförmiges Spannungsbild. 
o, verschwindet für »=0, d.h. wenn keine Quer- 
dehnung auftritt. Der maximale Wert von o, wird 
fürz=+A und v=0,5 d.h. für einen vo umen- 
beständigen Körper, zu 0,0 = 0,589 q, erreicht. Die 
Schubspannungen verschwinden, wie aus Symmetrie- 
gründen zu erwarten war, für y=rsing= (0 und 
sind innerhalb der Berührungsfläche unabhängig von x. 
Der Höchstwert liegt für y= + A zwischen T,,o = 
0,589 g, und 0,785 g, entsprechend v= 0,5 undv =. 
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Vohburg/Donau. G. Sonntag. 


Quadratische Interpolation bei großen Dif- 
ferenzen. 


Die gewöhnlichen Interpolationsformeln versagen 
bei Funktionswerten mit großen Differenzen; nur 
Thieles Kettenbruchentwicklung !) mit ‚„Reziproken 
Differenzen‘‘ ist dann brauchbar, worauf bereits 
Milne-Thomson?) hingewiesen hat. 

Für die praktische Anwendung ist nun nicht nur 
das Rechnen mit Reziproken Differenzen umständlich 
(wegen der vielen Quotienten), sondern es existieren 
hier aurch keine Anwendbarkeitskriterien wie bei den 
gewöhnlichen Interpolationsformeln. In der Praxis 
nützt aber ein Formelergebnis wenig, dessen Genauig- 
keit nicht in einfacher Weise festgestellt werden kann ®), 
Nachstehend sollen deshalb sowohl für die direkte wie 
für die inverse quadratische Interpolation ‚Formel- 
paare angegeben werden, die den gesuchten Wert ein- 
schließen und außerdem nur die Funktionswerte ent- 
halten, 

Die Reziproken Differenzen werden definiert durch 
die Ausdrücke?) 


*) Thiele: Interpolationsrechnung. Leipzig1909, 8.139. 

®?) Milne-Thomson: The Caleulus of Finite Differences. 
London 1933, p. 108. . 

®) Das betrifft auch die numerische Quadratur; der Verfasser 
wird hierauf bei nächster Gelegenheit ausführlich eingehen. 

*) Die Darstellung der Reziproken Differenzen durch spitze 
Klammern wird hiermit vom Verfasser in Vorschlag gebracht; 
die von Thiele eingeführte und von Nörlund über- 
nommene Schreibweisehalteichnicht nur für umständlich, son- 
dern auch für unzweckmäßig, 


ORTE 


2, —% TRIER EEG 
<YoYıyao?) = REES, +Yı Ben: % 
. Hieraus ergeben sich die Identitäten 


zz —%g 
TUT ya 


Pine A 

“u = ey? 7 we u 80° 
L—2 Kr 

KYYoHı? = Hoya ya) — Zyyayayd pa 
Beschränken wir uns auf eine Annäherung zweiter 


Ordnung, so erhalten wir nach Einsetzen und leichter _ 
Umformung die Interpolationsformel 


fla+M=yots7 Tuer Fr 


-Ya—2Yı+%o Kia 
= m— 04, 7 
er s Y2%0 j a 


und u = h/o (» Argumentschritt der Tafel). 2 
Gehen wir anderseits von y_, aus, so ergeben sich 
die Identitäten 


1. z — lg 
gg 


i 1 TUT j " ° 
f ö h ” 
En ee ba nn nn 


2 —%o 
Wy-r= YAIr Ft yg yo ya’ 


7ı — 
YY-190 = Y-1Yoyı) — Zyy ya ya pn Yod 


USW. 


Beschränken wir uns auch hier auf eine Annäherung 
zweiter Ordnung, so erhalten wir nach Einsetzen und 
Umformung die Interpolationsformel 


en 


—Y-4 


Sa+m=-yatttW | * 


wo 


A ee ee ei 


Da die Identitäten der zweiten Zeilen verschiedene 
Vorzeichen haben, schließen dig abgeleiteten Näherun- 
gen den gesuchten Wert ein. 


Als Beispiel soll tg 1,5685 berechnet werden. Mit 
den Funktionswerten 


tg 1,567 — 263,41126 tg 1,569 =: 556,69098 
tg 1,568 — 357,61106 tg 1,570 — 1255,76559 


ergeben die Formeln (1/2) bzw. (3/4) die Werte 
435,47730 435,47734 
und hieraus folgt durch Mittelbildung 
tg 1,5685 = 435,47732. 


bzw. 


Die Formeln für die inverse Interpolation finden 
sich sehr einfach durch Auflösung der Glen, (1) und (3) 
nach x; wir erhalten so 


(5) 


und 


wo x und A nach (2) und (4) zu berechnen sind. 
Göttingen, F,K. Rubbert, 


ionsverlaufes gewisse rechnersiche 


gliches Maß gebracht. 


Beil 
nn ne Funktion en 


. und zwar springt der Funktionswert, wenn man sich, 
- von positiven x -Werten kommend, dem. Wetz= ö 
nähert, im ‚Augenblick des Durchgangs durch Null 

 vony=0aufy=I1, um sich dann in Richtung der 


Zu u nähern. 


E7 | 0506 1 0 


Merkwürdigerweise ist bei vielen -Autoren, die diese 
E: - Funktion angeben, das dazugehörige Funktionsbild 
’ nicht immer einwandfrei dargestellt. Dies mag z.T. 
daher kommen, daß man den z-Wert zwischen« = — 1 
und = +1 vielleicht aus Gründen des erforderlichen 
größeren Rechenaufwandes nicht die nötige Beachtung 
schenkt. Nimmt man jedoch die Berechnung dieser 
Werte auf graphischem Wege vor, so läßt: sich inner- 
halb kurzer Zeit für beliebige Werte a zwischen 1 ‚002 
„und 50 der y-Wert bestimmen, womit man den Ver- 
lauf der Kurve zeichnen kann. 

Wenn auch die Gleichung in der vorliegenden „Form 
nicht unmittelbar in ein Nomogramm übersetzt werden 
kann, so liefert’eine Umformung jedoch sehr bald einen 
brauchbaren Ansatz, wenn man wie folgt vorgeht; 


en 1_1_ 008 


26) 


1 
=-Iga 


fordert: a nt mung 

Diese werden durch graphische. 
hr: Erklärung des Bogrikfen der Unstetigkeit 
der -Ditferentialrechnung findet man als Beuepiel, 


- gleichungen der eihzelnen Leitern hinschreiben nn 
“ lauten: \ 


BR rh T 2 "| x 
Be En Ira geh 


ie areen hat in 2 =0 eine Unstetigkeitsstelle 
£ Di danach gezeichnete Noirogrhane zeigt Bildl, 


a3 negativen z-Werte. dem ee y=05 re, 


Ir © = sr lg I _ ))- = lg 1g a u) VER, 08 
w 6% . (Grundform) EN. 


(2) ist ber eine elschmeklorn, die man als Leikör 
tafel mit drei parallelen Teilungsträgern ansetzen kann, 
wie die daruntergeschriebene Grundform zeigt!).. 


_ Durch Zuordnung kann man unmittelbar die Träger- _ + 


=2u 


Em, j (en ) = 


mittels dessen man, wie oben gesagt, den Verlauf der 
Funktion genügend genau feststellen kann, hir 
nach dieser Methode gefundene Kurvenbilder für 
s=Ll; 1,5; 2; 5 und 10 a 


095 


497 
19093 


Bild 3 SSH 


Wünscht man den Verlauf der Kurve bis zum 
Werte x = + 0 zu verfolgen, so würde sich’die Leiter- 
tafel bis ins Unendliche erstrecken müssen. Durch 


2) P.Luckey, Nomographie. Mathem.-Navurwiss. Bi-- 
klioth. Reihe I, Bd.59/60. B.G.Teubner, Leipzig 1937. — 
H.Schwerd t, Lehrbuch der Nomographie, $34.J. Springer, 
Berlin 1924. —M.Zühl ke, Rechentafeln und Sonderrechen- 
stäbe. III. Aufl. 1941. RKW-Schrift 116. B. G. Teubner, Leipzig: + 


ojektive Verzerrung kann man aber diesen Punkt 
des Zeichenfeldes abbilden und ist, damit 
in der Läge, die des Wertes y bis zum 
Erreichen des Wertes y = 0,5 für = +0 zu ver- 
wie Bild 3 zeigt . 

azu ist noch folgendes zu sagen: 

Die Gleichung (2) wurde mit positiven x-Werten ent- 
wickelt, wir können jedoch das gleiche Nomogramm 
auch zur Berechnung der y-Werte benutzen, wenn # 
d! negativ ist und zwar unmittelbar durch Verzifferung 
: der x- und y-Leitern. Es ist nämlich 


1 Le: 
er a Te 
Ir aUR 1 WEN (3) 


d.h. die Teilung der x-Leiter gilt sowohl für positive 
wie auch für negative Werte, nur die Werte der y- 
Leiter ändern sich und zwar so, daß die Summe der 
an den Teilstrichen angeschriebenen Werte gleich 1 ist. 
Im Nomogramm Bild 3 ist diese Tatsache bereits be- 
rücksichtigt. 1 

Den: Verlauf der Funktion nach = + ® kann 
man auf gleiche Art und Weise bestimmen, doch sei 
wegen des grundsätzlich gleichen Weges nicht näher 
darauf eingegangen. 


Zusammenfassung: 


Der Kurvenverlauf, besonders in.der Nähe von Null 
ist nur durch umfangreiche Rechnungen zu bestimmen. 
Ersetzt man diese jedoch durch eine en Me- 
thode, so ist der Kurvenverlauf sehr schnell Baden, 
wobei die Annäherung an den Grenzwert besonders 
deutlich wird. Durch projektive Verzerrung kann man 
den Wert x= +0 in das Zeichenfeld hereinholen. 


Berlin-Treptow. Johannes Fischer, 


Der „natürliche Maßstab“ von Schaubildern 
empirisch gefundener Funktionen. 


(Ein Beitrag zur Bestimmung der 
ausgleichenden Geraden.) 


Hat man, um einen gesetzmäßigen Zusammenhang 
zwischen zwei Größen / und /’ zu bestimmen, durch 
Versucheeine Reihe von Wertepaaren (2,14), (l,, 23), --- 
(ln, In) ermittelt, so werden in einem Schaubild mit den 
Abszissen 17 und den Ordinaten / die ‚beobachteten 
Punkte“ (1,7), i=1,2,...., n, nicht streng auf der 
gesuchten Funktionskurve liegen, sondern sie werden 
eine „fehlerzeigende Punktreihe‘‘!) bilden. 

Man hat es nun völlig in der Hand, im Schaubild 
durch eine entsprechende Wahl der Achsen-Maßstäbe 
oder der „Überhöhung‘‘ der Ordinaten gegenüber den 
Abszissen die vorhandenen Widerspruchs- oder Fehler- 
beträge deutlich aufzuzeigen oder zu unterdrücken. 

Im folgenden soll nun nachgewiesen werden, daß 
es unter der Vielheit aller diesbzgl. Maßstabsverhält- 
nisse ein ausgezeichnetes gibt, das mit den mittleren 
Fehlern #’ von 1’, bzw. « von ! im Zusammenhang 
steht. Es soll dieses als das ‚natürliche‘ oder auch 
als das „isometrische‘‘ Maßstabsverhältnis bezeichnet 
werden. 

Im folgenden soll nun angenommen werden, daß 
1. die der Messung zugänglichen Größen I und I’ im 

Grunde eine lineare Gesetzmäßigkeit befolgen (Bei- 

spiel: 7 = die Länge eines Metalldrahtes, I’ = die 

Draht-Temperatur), 

2, der Ausgangspunkt für die Zählung der I und /’- 

Werte in den geometrischen Mittelpunkt (= Schwer- 

punkt) des Punktsystems verlegt worden sei, so daß 


M=0, mM=0 (M. 


Infolge der,unvermeidlichen Messungsfehler sind nun 
sowohl an den I! die Verbesserungen v als auch an den 
I’ die Verbesserungen v’ anzubringen. Dann muß nach 
erfolgter Verbesserung oder „Ausgleichung‘‘ der 


‘) Vgl.P.Werkmeister: Ermittlung derplausibelsten Ge- 
raden einer fehlerzeigenden Punktreihe, Z,angew. Math. Mech. 
1921. 8.4911, 


talk ala 1 
Er Bru=eirt 
i=1,2,...,n;tund ssind die noch z 
ist r3 
der Fur treih 
dieser Geraden auf der l-Achse b GI 
indie Der Neigungswinkel 9 dieser Geraden ist dann 3 
g=artgti.....- le + (3). 
einen Näherungswert t, ein, so daß 
t=4+Ä4t Pag ee (4), 
so ergeben sich aus (2) die Fehlergleichungen 3% 
b+v=lW+ANlÜHV)+s 


mte 


Führt man für ı 


oder: 
4=v-hvi=l-At+s4+hi—h+At: Vi (). 
Zunächst wollen wir das letzte Glied als kleine Größe 
II. Ordnung vorü wir werden 
es jedoch im Ver. der Rechnung wieder in Rück- 
sicht "ziehen. . 

Will man die Unbekannten A i und s nach der Me- 
thode der kleinsten Quadrate bestimmen, so man 
die Größen 4; als fingierte Verbesserungen?) anzu- 
sehen. Ihnen sind die Gewichte 
; „ont 

pi = EFRn® Fe I 
beizugeben ?). Über die willkürlich wählbare Gewichts- 
konstante wird hier zweckmäßig so verfügt, daß 

= Pi = 1 (6b). 

Dann hat man bekanntlich die Normalgleichungen 
aufzustellen; sie lauten: 

} m. 


[7]-At+[7-s + FI) —RrI=0 
(8), 


(6a) 


ya. 


W-At+ ns +4) —M =0 
Im Hinblick auf (1) wird aber: 


pen, 


1%] 


t= 


d.h. se 


At=—b+ 


sowie 
(9); 


d.h. die ausgleichende Gerade muß durch den ge- 
wählten Ausgangspunkt, den Schwerpunkt des Punkt. 
systems, verlaufen (Lambertscher Satz %)). Die fin- 
gierten Verbesserungen A, lassen sich dann berechnen zu: 


ae 
1 
Spaltet man diese nach Maßgabe der entspr. mittleren 


Fehlerquadrate auf?), so gelangt man zu den Ver- 
besserungen v; und v': 


yell—l Kl: 73,70). 


FF ER zus tu’? 
und ® = Fr yadi (11). 
Nunmehr sind wir in der Lage, das in (5) vernach- 
lässigte Glied zu berücksichtigen und die Ausgleichung 
von neuem zu beginnen. Dann hat man, falls wir den - 
zuletzt erhaltenen Wert für t als neuen Näherungs- 
wert i, ansehen, so daß 


ti=W+At 
») W.Jordanu. O,Eg 
kunde. 1. Band. Stuttgart 1 
mehreren Beobachtungsgrößen. 
°») H.Wolf: Boitrag zur Frage der strengen oder näherungs- 
weisen Ausgleichung von verschiedenartigen Beobachtungen. 
Berlin 19836, 
‘) E.Hammer: Bestimmungeiner Geraden aus gemessenen 


Koordinaten einer Anzahl ihrer Punkte, Zeitschritt für Ver- 
messungswesen 1925, 8, 52ff, 


FR. 
ERFRRRE 


gert:Handbuchder Vermessungs- 
935, 8.183, Fehlergleichungen mit 


Pr} 


+ür u 


3 Er (io wer = @ ä su — [1] — @ 7 un) =0 
Man könnte nun ‚hieraus At errechnen, diesen Wert 
= t, anbringen, so daß man einen verbesserten Wert 
= für 2 ‚erhielte, diesen wiederum als neuen Näherungs- 
REN Werb: betrachten und ihn erneut in die Ausgleichung 
er; 2 Se bis die Rechnung in der gewünschten De- 
 ..... zimale zum Stehen gekommen i ist. Dieser Endzustand: 
ist dadurch charakterisiert, daß At, t, bzw. At gleich 
DR = - Null wird. Der zugehörige Wert i,für den Steigungs- 
Be ‚faktor muß dann notwendig mit: dem definitiven Wert 
' für i zusammenfallen. Aus (15) ergibt sich dann als 


de use 


Pl] + Kur — un) — un =0... (16), 
ie a 
ee . wr] 

Be = 210 a BR ] 
ER | [217] vg em (17a). 
. Erweitert man mit y 

RE _ WI 71 TE 
mE 7 nym- =) —. = 


so erhält man auch 
| 211 
me 1m F ya —euny + sg ij ir 


Die en für die. Quadratwurzel hat 
K. Friedrich?) in dem Sinne beantwortet, daß 
stets das positive Vorzeichen gilt. Nur in dem zu- 
gehörigen Ausdruck für [A A] steht das negative Zeichen 
vor einer entsprechenden Wurzel:. Aus der bekannten 
"Kontrollformel für die kleinste Quadratsumme der 
verbleibenden Verbesserungen 


[44] = [il] — [27° 


a folgt nämlich: 
- gr 2) 7 = 7 
Bee Bilder + Syn EN + T)as) 


2eZ Wir wollen nun einen Sonderfall näher betrachten, 
nämlich wenn i 
u’=y wird, 


N 


d,h. q=1. 


Dann wird: 
FT 
= 1 = ao | 
— [77] + YEN— TIP + EL) 9) 


5) K.Friedrich: Mechanische Deutung und Bestimmung 
der auszleichenden Geraden. Zeitschrift für Vermessungswesen 
„1981, 8.129ff. : 


(19a), 


W-Yer ENDE in diesem Fälle 1a EN ‚noch a weite 
ey formung Rn ! et 


zu unterziehen, daß die veränderten Koordinaten- 


a 13 

A u “ R Si 

aditr + N Sr an 
€ ne Iautet, somit: 

ul [7 gi u 2 ty g2fll a 

en a li 14 g { 1+ ei D+ ns 


zu Baleen. Der zugehörige Steigungsfaktor tist jedoch 


g für den ee, dieser. Annäherungs- 


wenn man die Ordinaten im Verhältnisq : 1, d.h. B Ki 


 „natürliche‘“ 


ae 


Mit ig= DZ 


wird närnkeh 7.77 
El ran 
tg2y= —.. ——  —— . 
u ET TEIGEn- llesn ne 
v ER a der Neigungswinkel in einem Koordinaten E 
system mit gleichgenauen Koordinaten. 
Es ist nun naheliegend, die gegebenen Voreikieh 
genauen Koordinaten einer solchen Transformation 


werte, nachstehend mit z und y bezeichnet, die Eigen. 
schaften von gleich-genauen Koordinaten Ba 
Man braucht hierzu nur er 


ls und» v4: I 


nicht der Er wie im I-/’-System, sondern es gilt; 
Sr Eee (216,8 
Hat man die I und /’-Koordinaten nach (21) in das ° 
x-y-System übertragen, so braucht man nunmehr, 
‚wie in (19a) und en nur Be den Sg a Be 
rechnen: ; 


= IE (ty) 1a + IF Fee ea, 


bzw. 


SE 2 (& y) BE 

BI—[2] + Ye) —[y?))? + leyl? 
oder bequemer: - 
291, 
| war gi 


Überträgt man den Vorgang ins Geometrische, so folgt _ R 
offenbar: 
Die Koordinaten werden gleichwertig oder „isomelrisch““, 


. (23). 


gegenüber den Abszissen überhöht, bzw. die Abzsissen . 
im Verhältnis uw: gegenüber den Ordinaten reduziert. 

Der Maßstab, in dem dann die gemessenen Koordi- 
naten: zur Darstellung gelangen, kann daher als der 
eder der ‚‚isometrische‘‘ Maßstab be- 
zeichnet werden. 

Man findet auch leicht die folgende Formulierung: 

Das Schaubild ist’so anzulegen, daß die Beobachtungs- 
fehler der Abszissen und Ordinaten linear als gleich 


große Strecken erscheinen. : 2 
Nur nach einem so ausgeführten Transformations- Er 
vorgang darf man nach der klassischen Schumann- rt 


schen Formel (23) rechnen, wobei man dann auf die 
mittleren Beobachtungsfehler keinerlei Rücksicht 
mehr zu nehmen braucht. Man hat lediglich zum 
Schluß noch den Übergang von =tgy auf i nach ; 
(21b) zu vollziehen. (Selbstverständlich könnte man . 
auch unter Außerachtlassung dieser geometrischen 
Verhältnisse nach 


1 2gll). 
a een 


bzw. nach N ‘oder (17b) rechnen.) 

Man erkennt übrigens, wie hier das Geometrische 
zu einem bloßen Spezialfall eines allgemeineren 
Problems wird, das nur auf den funktionalen Zu- 
sammenhang zwischen zwei Größen hinzielt. 


% (24), 


°%) R.Schumann: Bestimmung einer Geraden durch Aus- 
gleichung der beobachteten Koordinaten ihrer Punkte nach der 
Methode der kleinsten Quadrate. Sitzungsberichte der Kaiser. 
lichen Akademie der Wissenschaften, mathem.-natuıwissersch 
Klasse. Wien 1916. 
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Es.erscheint nun noch notwendig, einen Blick auf 


- die beiden Grenzfälle mit u” = 0 und a= 0 oder g = 0 


und g= ® zu werfen. Dann sind nämlich nicht die 
Abszissen, bzw. Ordinaten unendlich groß aufzutragen, 
sondern die Formeln schrumpfen vielmehr zusammen 
auf: 


a) u =0: Ku'=0) = Ip "2. + (258) 
b) u 0 u) = 0... (206). 


Sodann besteht in diesen beiden Grenzfällen keinerlei 
Notwendigkeit mehr, die ursprünglichen Koordinaten 
1, U’ einer Transformation zu unterwerfen, und es kann 
dann die „Überhöhung‘‘ beliebig groß gewählt werden. 
R.Schumann) hat erstmalig bewiesen, daß im 
Falle eines echt-linearen Zusammenhanges 7 = f(l’) 
der Differenzbetrag 
At=ta=0)—tu=0) (26) 


nur von der Größenordnung der kleinsten Quadrat- 
summe [44], d.h. von der Größenordnung kleiner 


Größen II. Ordnung ist. 


Es soll nicht unerwähnt bleiben, daß die in der 
Korrelationsrechnung benutzten „normalen‘‘ Koordi- 
naten ?) sich als ein Sonderfall der hier angegebenen 
„natürlichen‘‘ oder „isometrischen‘‘ Koordinaten dar- 
stellen lassen, nämlich insofern, als dort von der An- 
nahme ausgegangen wird, daß — entsprechend der hier 
verwendeten Bezeichnungsweise — 


ya 
n 


Wie wenig die hier dargelegten Beziehungen nu 
theoretisches Interesse besitzen, sondern auch praktisch 
bedeutsam sind, zeigt ein von E.Pinkwart?) 
angegebenes Beispiel: Er hat sich die Aufgabe ge- 


“stellt, durch vier genau auf einem Kreis 


gelegene Punkte die ausgleichende Gerade zu be- 
stimmen. Da diese 4 Punkte achsensymmetrisch liegen, 
liefert die Anwendung der Gleichung (23), bzw. (22a) 
oder (22b) keine eindeutige Lösung, sondern sowohl 
die senkrechte, wie auch die waagerechte Gerade 
durch den Kreismittelpunkt als Lösung, was auch 
durchaus plausibel erscheint; da es nach vorstehendem 
nämlich ein leichtes ist, durch eine entsprechende 
Maßstabswahl für die Ordinaten, bzw. Abszissen die 
auf einem Kreis gelegenen 4 beobachteten Punkte 
einmal zur waagerechten und zum. anderen zur senk- 
rechten Geraden nahezu „zusammenzuquetschen‘‘, 

Wären jedoch a.a.O. gleichzeitig auch noch die 
mittleren Koordinatenfehler u, bzw. u’ mitangegeben 
worden, so wäre — falls einer der Grenzfälle (25a) 
oder (25b) gemeint gewesen war — klar ersichtlich, 
daß es nur eine einzige Lösung gegeben haben würde, 
nämlich (25a) oder (25b) °)/ 

Wenn jedoch von E. Pinkwart ein Fehler- 


verhältnis 
0<1<» 


gemeint ist, so ist eine eindeutige Lösung gar nicht zu 
erwarten: Einmal scheitert es an der grundsätzlichen 
Unmöglichkeit, eine achsensymmetrische Punktreihe, 
diestrengeiner Kreisgleichung genügt, durchein lineares 
Gesetz annähern zu wollen (eine gerade Funktion, nicht 
eine ungerade, hätte hierbei zugrunde gelegt werden 
müssen!), zum andern ist es das Fehlen der mittleren 
Fehlerwerte % und uw’; da nach Gl. (17a) oder (17b) 
oder (24) eine einfach-unendliche Mannigfaltigkeit von 


‚ Lösungen existiert, — je nach Wahl des Fehlerverhält- 


nisses q = : Durchläuft q alle Werte von Null bis 


’, F. Bauer; Korrelationsrechnung. Leipzig—Berlin 1928. 

?, E. Pinkwart: Nochmals die Bestimmung einer Geraden 
aus den gemessenen Koordinaten ihrer Punkte, Zeitschrift für 
Vermessungswesen 1942, 8. 2171f. 

°) Vgl.H.Wolf: Beitrag zur Bestimmung der Gleichung der 
plausibelsten Geraden einer fehlerzeigenden Punktreihe, Zeit- 
schrift für Vermessungswesen 1041, $8.4111f, 


‚--1/R ,.. (en. 
Be en. 


Lin in- 


Hioser Monniglaltigkeit, Berausgeguiffenen. Wert dat. 


leren Beobachtungsfehler unerläßlich ist, um so mehr, 
wenn die urs ichen ee gr 
von verschi Dimension sind (z. B. Längen- und 


Temperaturwerte, Winkel- und j m 
so daß gleichsam über den „natürlichen“ oder „iso- 
metrischen‘‘ Maßstab des Schaubildes eindeutig ver- 
fügt werden kann. 
Bamberg. H. Wolf. 


») E,Lindinger: Über die a Gerade. Zeit- 
schrift für Vermessungswesen 1943, 8.135ff. 


Nomogramme für die Gleichung "4. Grades 
mit reellen oder komplexen Wurzeln. 


Für dienumerische Lösung einer Gleichung 4. Grades _ 


wurde von H.Heinricht) ein praktisches Ver- 
fahren angegeben, das die reellen und komplexen 
Wurzeln mit, beliebiger Genauigkeit zu berechnen ge: 
stattet. Die eigentliche Rechenarbeit besteht da 
in der Bestimmung einer — stets positiven — reellen 
Wurzel der kubischen Resolvente der Ausgangs- 
gleichung. Auch bei dem Verfahren von Ferrari?) 
ergeben sich die Wurzeln der ‘Gleichung 4. Grades 
dadurch, daß man zunächst eine reelle Wurzel einer 
bestimmten kubischen Gleichung ermittelt. 

Wenn — wie es z.B. bei Stabilitätsuntersuchungen 
häufig vorkommt — sehr viele.ähnliche Gleichungen 
4. Grades (bei systematischer Variation der Koeffi- 
zienten) zu lösen sind, kann man die erforderliche 
Rechenarbeit durch Benutzung von Nomogrammen 
vermindern. Als Nachteil muß man dabei eine be- 
schränkte Genauigkeit in Kauf nehmen; als Vorteil 
gewinnt man’ aber die Möglichkeit, ohne weitere Rech- 
nung aus dem Nomogramm das Verhalten der Wurzeln 
in der Umgebung des untersuchten Punktes qualitativ 
zu erkennen. Mir sind bisher in der Literatur nur 
Nomogramme für die Gleichung 4. Grades bekannt, die 
entweder die reellen Wurzeln?) oder die komplexen *) 
liefern. Deshalb seien im folgenden Nomogramme 
angegeben, die in gleicher Weise zu den reellen wie zu 
den komplexen Wurzeln führen. 

Die Gleichung 4. Grades ji 


+ Hr Her +d=0 ..- ll) 


wird zunächst durch y=x + 4 a auf die reduzierte 


Form 


„+b®+ay+d=0.:...0) 


gebracht und dann durch x = er weiter vereinfacht 
6 
zu 
A+s2+rz+ß=0.,....% (3). 


In einer Netztafel mit der Abszisse x und der Ordinate ß 
erhält man für jede reelle Wurzel eine Gerade. Das 
so gewonnene Nomogramm (Bild 1) zeigt deutlich die 
drei Bereiche, in denen vier, zwei oder keine reellen 
Wurzeln existieren. Das Ablesen der Wurzeln ist aber 
in der Nähe der Grenzen dieser Bereiche und in dem 


. )H.Heinrich: Zurrechnerischen Auflösung einer Gleichung 
vierten Grades. Z. angew. Math. Mech. 1941, S. 304. 
?®) Vgl.F,C. Haus ; Stabilit6 et Maniabilit6 des Avions. (Paris 
1930), $. 310, 
°)H. Schwerdt: Die Anwendungen der Nomographie in der 
Mathematik. (Berlin 1931), 8. 74. 
*) F.A pt: Nomogramme zur Bestimmung der komplexem Wur- 


zeln einer beliebigen Gleichung vierten Grades, Z. angew. Math, 
Mech, 1932, S, 382, 
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t De TERN er A Re a n ‚ Mi 
. Hierin ist 'w stets reell, während ‚v, und v, entweder 
‚reell oder rein imaginär sein können. Falls zwei oder 


2 ar 


essere 


Ei haslüge ‚Erfassung der reellen und komplexen 
Wurzeln erhalten wir auf folgende Weise: Wir setzen 


‚die vier Wurzeln, wie es schon De scartes3)getan 
hat, in folgender Form an; "I 0, 


Ad a A 6 
DE et \ A) 


2, a=-—u4r% 


vier reelle Wurzeln vorhanden sind, werden also auch 
diese durch den Ansatz (4) paarweise zusammengefaßt. 


ch mit vier reellen Wur- ‘oder 3 allein auskommen, da 


Ablesemöglichkeiten und außerdem a die 


. praktisch. 5 1 or 
‚In dem Bereich mit vier'reellen Wurzeln erhält man 
‚durch jeden. Punkt drei Kurven u = const., dadesdeii 
ı verschiedene Möglichkeiten gibt, die vier Wurzeln 
zu Paaren zusammenzufassen. Für die Lösung der Gei- 
chung 4. Grades genügt es aber, auchin diesem Bereich 


OR 


H 1 
2—pHt- j 
Er Sn EIER a, 
ist. Da jedoch die Ablesegenauigkeiten dieser beiden 


Netztafeln für v? und v} sich gegenseitig gut ergänzen, 


ist die Verwendung beider Tafeln nebeneinander recht 


mit nur einem Wert uzurechnen. In den Bildern 2und 
3ist daher nur eine Schar dieser Kurven ausgezogen. _ 
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Für die reellen Größen u, v? und v3 erhält man durch 


. Koeffizientenvergleich von Gleichung (3) mit 


| (? — 21) (— 2) @— 2)@ 2) =0 


folgende Beziehungen: 


1 

4B = (a ART 
EP Fe 45) 

z 4u 

1 r 1 

er ET 


Eine Vertauschung von u mit —u zieht auch einen 
Austausch von v? und v3 nach sich. Positive Werte 
von v? und v2 ergeben reelle Wurzeln z, negative Werte 
dagegen komplexe. - 

Die Gleichungen (5) stellen wir nün in einer Netz- 
tafel dar, deren Abszisse x und deren Ordinate ß ist. 
Um die Ablesung zu erleichtern, zeichnen wir für 
jedes Wurzelpaar ein gesondertes Nomogramm (Bild 2 
und 3). Man könnte natürlich auch gut mit. Bild 2 


5) Vgl. E. Netto: Vorlesungen über Algebra, (Leipzig 1896), 
90. IR 
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Bei. der Herstellung der Nomogramme ist zu. be- 
achten, daß die. Kurven u = const. sämtlich von der 
gleichen Parabel gebildet werden, deren Scheitelpunkt 
auf einer Hyperbel wandert; die Verwendung einer 
Schablone für diese Parabel erleichtert daher die 
Zeichenarbeit erheblich. Für das Zeichnen der Kurven 
v? = const. mit konstantem Av? nutzt man den Um- 
stand aus, daß ihre Schnittpunkte mit den Kurven. 
u = const. aequidistante Abszissen haben. 

Statt der Kurven v2? = const. hätte man auch Kurven 
v = const. zeichnen können... Ich habe die ersteren 
vorgezogen, da die Kurven v = const. ungleichmäßigere 
Abstände bekommen, die das Interpolieren erschweren. 

Die Bezifferung in den Bildern 2 und 3 ist der be- 
quemeren Handhabung wegen so gewählt, daß die 
gesuchten Wurzeln in beiden Fällen durch u + v ge- 
geben sind. Mit dem Vorzeichenwechsel in u ist der 
Übergang von v, zu v, oder umgekehrt verknüpft. 

Der erforderliche Rechenaufwand sei an einem 
Beispiel gezeigt, das auch an anderer Stelle 6) zu dem 
gleichen Zweck benutzt wurde: 

a—42°+52° +52 —10=0. 
Die Anwendung des Hornerschen Schemas ergibt 
mit y=r—1: 


®), F.A. W illers: Methoden der praktischen Analysis (Berlin 
1928), 8.209 und H. Heinrich, loc. eit, - 
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Bei der Anfertigung der ee ar Bild? 


und 3 ist es notwendig, den Bereich der &- und ß-Werte 


‚den.zu erwartenden Werten der Koeffizienten. der 
. Ausgangsgleichung anzupassen. 

- Falls die Nachbarschaft des Brenagchich rien: 
- halb der negativen a-Achse häufig gebraucht wird, 
_ empfiehlt es sich, hierfür ein besonderes Nomogramm 
"in anderem Maßstab herzustellen, in dem außerdem 


die von unten kommenden Kurven u = const. in das 7 
kleine dreieckförmige Gebiet hinein durchgezeichnet 


sind. Dann fallen auch hier die in der Nähe von Ein- 
a immer vorhandenen Ableseschwierigkeiten 
o 

- Falls die Anfertigung der Nomogramme nach Bild 2 
> ne 3 zu mühsam erscheint oder nur. verhältnismäßig 
wenige Gleichungen gelöst werden sollen, besteht noch 
"folgende Möglichkeit einer teils" graphischen, teils 
numerischen Lösung: Man zeichnet auf einem Blatt 
die Parabel 16 (u?) über u? als Abszisse und auf einem 


zweiten transparenten Blatt die Hyperbel en über 
u? in gleichem Maßstab wie die Parabel. Beingb man 


jetzt durch Parallelverschiebung den Nullpunkt des 


Hyperbel-Blattes über den Punkt I: 4ß des Parabel- 
- Blattes (Bild 4), so-gilt en für den Schnitt- 
punkt (oder die ee der Parabel mit der 
 Hyperbel die Beziehung: - 


* . r 


Ra u En der Kbadine des Ra Blakten, Br 


zuleeR ist. Ns ee DB mit. = ersten 


v? und v2. Da man auf dem Parabel-Blatt keine Be- 


zifferung (mit u und 16 u*) benötigt, kann man Ab- 


szisse und Ordinate mit den Werten & und ß beziffern 


und sich dadurch die praktische Handhabung erleich- _ = 


tern. Man erkennt leicht, daß für negative & auch drei 
Schnittpunkte zwischen Parabel und Hyperbel auf- 
treten können; in diesem Falle ergeben sich wie oben 
vier reelle Wurzeln der Gleichung. 

Bei der Reduzierung der allgemeinen Gleichung 
4. Grades auf die Form (3) behalten die Wurzeln ihre 
Eigenschaft, reell oder komplex zu sein; nicht erhalten 
bleibt dagegen i.a. das Vorzeichen des 'Realteils. 


Braunschweig. =D H.Blenk. 


BUCHBESPRECHUNGEN 


Dr. Siegfried Flügge (o. Prof. a. d. Univ. Marburg) 
unter Mitarbeit von Dr. Hans Marschall, Rechen- 
methoden der Quantentheorie, dar- 
gestellt in Aufgaben und Lösungen. 
I. Teil. Elementare Quantenmechanik. (Die Grund- 
lehren der. mathematischen Wissenschaften in Einzel- 
darstellungen. Band LI). X.+ 240 S. mit 18 Abb. 
Berlin u. Göttingen 1947. Springer-Verlag. Preis 
brosch. 18,— DM, 


Das vorliegende Buch bezweckt, wie der Verfasser 
schon im Vorwort. betont, den Lernenden mit den 
komplizierten mathematischen Hilfsmitteln der Quan- 
tentheorie vertraut zu machen und die Scheu vor ihrer 
Handhabung zu beseitigen. Dies ist den beiden Auto- 
ren auch in ausgezeichneter Weise ‚gelungen. 

Heute mangelt es nicht an einfachen Darstellungen 
der Quantentheorie; es ist aber für den Lernenden oft 
. unbefriedigend, sich mit Plausibilitätsbetrachtungen 
begnügen zu müssen. Auf der anderen Seite gibt es 
eine Anzahl strenger Lehrbücher und Handbuch- 
artikel über Quantentheorie, deren Rechnungen nur 
mit Mühe von einigen Eingeweihten verstanden 
werden. 

Dieses Buch zeigt nun in Form einer Aufgaben- 
sammlung, wie ein Problem angepackt und durch- 
gerechnet werden soll. Gleich am Anfang. kann z. B. 


der Lernende die 4. Aufgabe „Eindimensionaler Po- 
tentialtopf“ ohne große Mühe nachrechnen, und es 


. wird weiter an der.Aufgäbe 16 der Behandlung der 


periodischen Potentials "ohne Anstrengung folgen 
können. ‚Das gleiche gilt auch für die anderen Auf- 
gaben. Das Problem der gekoppelten linearen Os- 
zillatoren, das quantentheoretisch de van der 
Waalsche Bindung verständlich macht, ist in ele- 
ganter Weise dargestellt. Ebenfalls wird die che- 
mische Bindung (Wasserstoffproblem) bis zum Ende 
durchgerechnet, so daß der Leser sich nicht nur mit 
einer allgemeinen Lösung abfinden muß. 


Die Auswahl der gewählten Aufgaben ist vortreff- 
lich. Man findet praktisch jedes wichtige quanten- 
theoretische Problem gelöst, soweit das mit einfachen 
Hilfsmitteln möglich ist. Ein Wunsch möge noch ge- 
äußert werden. Es sei zugegeben, daß das Buch den 
Fehler vermeidet, den Leser mit einer Flut von physi- 
kalischen Begriffen zu überschütten, da ja die Ver- 
fasser das Hauptgewicht auf die Rechenmethoden 
legen. Andererseits wäre es auch für einen mehr ma- 
thematisch orientierten Leser begrüßenswert, wenn 
z. B. in Aufgabe 55 die zugrundeliegende Fragestellung 
etwas deutlicher auseinandergesetzt würde. Der Ler- 
nende würde erfahren, wozu das Rechnen gutist, und 
würde noch mehr Freude beim Lesen des Buches haben. 


der Gleichungen (5) überein. Durch Verschiebung des 

Hyperbel-Blattes über dem Parabel-Blatt kann man 
also zu jedem Wertepaar (a, ß) das zugehörige u?finden. 
Daraus ergibt sich dann mit Hilfe der Gleichungen (6) 


2 


Das Buch kann jedem, der beim Lesen der quanten- 
theoretischen Lehrbücher nicht weiterkommt, weil 
ihm das Rechnen zu große Schwierigkeiten bereitet, 
an erster Stelle zur Ergänzung empfohlen werden. 


Dresden. 0, Stasiw. ° 


Dr. M. (zerny (o. ö. Prof. a. d. Univ. Frankfurt) un- 
ter Mitarbeit von Dr. H. Müser (Wiss. Assistent.a. d. 
Univ. Frankfurt), Anweisungen zum Physi- 
kalischen Anfängerpraktikum. VIlL + 
1148, mit 63 Abbildungen, Frankfurt a. M.1948. Verla, 
von Dr. Dietrich Steinkopff. Preis brosch. 8,— DM. 

Die „Anweisungen für das Anfängerpraktikum‘‘ 
umfassen als Hauptinhalt die Beschreibung von 52 Ver- 
suchen. Durchweg handelt es sich dabei um Aufgaben, 
die in jedem Anfängerpraktikum stehen, so daß das 
kleine Buch wohl überall gut verwendet werden kann. 
Die Anweisungen sind so abgefaßt, daß die Studenten 
die Versuche möglichst selbständig ausführen können, 
Unnötig sind vielleicht die Zahlenangaben, die sich 
auf die Frankfurter Versuchsapparaturen beziehen, 
die also an anderen Stellen nicht. verwendet werden 
können Deneigentlichen Versuchsbeschreibungen sind 
allgemeine Bemerkungen vorangestellt über die Ge- 
nauigkeit von Messungen, Fehlerabschätzung, Dar- 
stellung der Versuchsergebnisse usw. Ohne Zweifel 
stellt das Buch gerade bei dem jetzigen Lehrbuch- 
mangel hinsichtlich Umfang und Inhalt eine sehr will- 
kommene Bereicherung der Lehrbuchliteratur dar. 


BRecknagel, 


Dr.-Ing. W. Oppelt, Grundgesetze der Re- 
gelung. (Bücher der Technik. Herausgeber: Dr.- 
Ing. Alfred Kuhlenkamp.) 118 S. m. 32 Bildern und 
28 Tafeln. Woflenbüttel-Hannover 1947. Wolfen- 
bütteler Verlagsanstalt GmbH. (Notdruck). Preis 
brosch. 10,— DM. 

Der Verfasser weist ein staunenswert großes Ge- 
schick auf, Beispiele und beispielhafte Anordnungen 
zu ersinnen und sie in schematischen Skizzen klarzu- 
machen. Das Büchlein stellt eine Fundgrube soleher 
anschaulicher Verdeutlichungen von Regelanordnun- 
gen und ihren Bestandteilen dar. 

Nicht auf der gleichen Höhe an Faßlichkeit und 

Klarheit hält sich der Text, Der Verfasser will ‚für 
den Lernenden“ schreiben. Mir scheint jedoch, daß 
er sich an vielen ‚Stellen nicht darüber im klaren ist, 
wie hoch er in Wirklichkeit die Ansprüche an Vor- 
kenntnisse stellt,: Ich glaube z. B. nicht, daß die Ein- 
führung des Frequenzgangs auf S. 31 auf ein unvor- 
bereitetes Gemüt anders als bedrückend wirkt. Der 
Stil der Darstellung ist an vielen Stellen weniger der 
einer Einführung für Anfänger als vielmehr der 
eines Referates in einem Handbuch. Billigt man die- 
sen, so fallen die meisten der Bedenken fort, die man 
sonst anmelden müßte. Aber auch noch nicht alle. 
Insbesondere die Fassung der Aussagen im Abschnitt 
über die Stabilitätsbedingungen gibt noch zu einigen 
Einwänden Anlaß. Sie im einzelnen anzumerken, ist 
wohl nicht nötig. 
Der Leser jedoch, der schon einige Kritik mitbrin- 
gen kann, findet in dem Büchlein eine überaus anre- 
gende und fördernde und hinsichtlich des Bildmate- 
rials auch erfreuliche Darstellung. 


Karlsruhe, K, Klotter. 


Dr. W. Weizel (o. Prof. a. d. Universität Bonn) 
Einführungin die Physik. 3. Band. Op- 


tik— Atomphysik — Wärme (Meyers kleine 


Handbücher Bd.44)- 2. Auflage, 156 S. mit 98 Abb. 
Leipzig 1947, Verlag: Bibliographisches Institut. Preis 
brosch. 2,90 DM. 

Unter den Lehrbüchern der Physik, die im Hoch- 
schulunterricht benutzt werden, gibt es viele umfang- 
reiche Werke, es fehlt aber an kurzen Zusammenfas- 
sungen. Mediziner, Chemiker, Land- und Forstwirte 


Lehrbuch das Wichtigste auszusue 
Be a | = “ Ir Beben ee Re 
en dure er 
Sushak, Dip Gesetze sollen glaubhaft gemacht werden, 
indem sie folgerichtig auseinander entwickelt werden. 
Zweck des Buches ist es, „Kenntnisse und Verständnis 
zu vermitteln, nicht aber Beweise zu führen“. 
Dresden. Becknagel. 


Dr. Fr. A. Willers_(o. Prof. a.d. Techn. 
Dresden; Elementarmathematik, ein 
Vorkurs zur höheren Mathematik. 
260 8., 172 Abb. Dresden und Leipzig 1948, Verlag 
Theodor Steinkopf. Preis geb. 16,— DM, brosch, 
14,— DM. 


Das Buch entspringt Kursen, die der Verf. vor An- 


fängern zur Vorbereitung auf das Hochschulstudium 
gehalten hat. Es gliedert sich in die Kapitel: Arith- 
metik und Algebra, Goniometrie und Trigonometrie, 
Analytische Geometrie der Ebene, Vektoralgebra, 
die klar und anschaulich geschrieben und mit vielen 
durchgerechneten Beispielen versehen sind. Seine Lek- 
türe setzt, die einfachsten Sätze der Schulmathematik 
voraus, bringt: aber mehr, als üblicherweise an den 
höheren Schulen geboten wird. Das Buch kann allen 
Ingenieurstudenten zum Durcharbeiten neben der 
Kursvorlesung in Mathematik empfohlen werden, zu- 
mal es eine Fülle von schönen Aufgaben mit Lösungen 
bringt, deren Durchreehrung eine gründliche Aneignung 
und Vertiefung der in der Vorlesung behandelten Be- 
griffe verspricht. 


Bonn. Wendt. 


Prof. Dr. Wilhelm Blaschke, Analytische 
Geometrie (Bücher der Mathematik und Natur- 
wissenschaften, herausgegeben von Dr. Henry Poltz). 
152 S. m. 67 Abb. Wolfenbüttel und Hannover 1948. 
Wolfenbütteler Verlagsanstalt GmbH, 1948. Preis 
brosch. 10,50 DM. ; 

Der leitende Gesichtspunkt ist — im Unterschied 
gegen die meisten neueren Darstellungen des Gegen- 
standes — geometrisch und nicht algebraisch. Es 
wird daher ausschließlich im Raum und in der Ebene 
gearbeitet, primär in jedem Falle im Raum. (Eine 
Verallgemeinerung auf die Räume höherer Dimension 
wird nur von der algebraischen Methode, nicht von 
den geometrischen Gegenständen gefordert und bleibt 
unvollzogen,.) Auf projektive Geometrie werden nur 
hier und da Ausblicke gegeben; sie soll in einem wei- 
teren Bande für sich behandelt werden. 

Nur ein Meister der Darstellung wie Blaschke kann 
es überhaupt wagen, auf einen Raum von 100 Seiten 
eine solche, für den Anfänger vielleicht etwas beäng- 
stigende Fülle von Stoff zusammenzudrängen, eine 
solche Fülle sowohl weittragender Ideen wie schöner 
Einzelheiten, Nach den linearen Gebilden werden 
Kugeln betrachtet, lineare Scharen von ihnen, Kon- 
figurationen und Transformationen, vor allem die 
stereographische, die durch reziproke Radien, die 
Geraden-Kugel-Transformationen, die Laguerresche 
Spiegelung und als Anwendung das Berührungspro- 
blem des Appollonius. Einen breiten Raum nehmen 
die Beziehungen zur Kinematik und Statik ein, mit 
einer ausführlichen Besprechung der, Stabwerke, 
Kräftepläne, Trägheitsmomente, an Hand deren die 
Hauptachsentransformation und die Säkulargleichung 
entwickelt wird. Die Theorie der Quadriken wird bis 
hinauf zu den konfokalen Scharen räumlicher Qua- 
driken, ihrem Reyeschen Achsenkomplex, sowie bis 
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auf den Ellipsoiden gefühtt. 0.0 
Verf. setzt materiell und methodisch nichts voraus, 
_ wohl aber eine sehr gründliche Schulung in moderner 


mathematischer Denkweise. Er verwendet überall 


diejenigen algebraischen Hilfsmittel, die die Beweise | 


' am kürzesten und durchsichtigsten erscheinen lassen, 
Vektoren, Gruppen, Matrizen, duale Zahlen, Quatern- 


ionen. Was er davon braucht, definiert und beweist‘ 


EN 


er. Freilich verbietet der knappe Raum jede Wieder- 
holung, jede Beleuchtung eines und desselben Gegen- 

standes von verschiedenen Seiten, jede Erklärung ad 
hominem, so daß der ungeübte Leser vielleicht einige 
Schwierigkeiten zu überwinden haben wird. Für den 


'. reifen Mathematiker ist das Buch materiell’ eine Fund- 
. grube, methodisch eine Bereicherung, in der Darstel- 


lung ein Genuß. 
_ Dresden, 


' Dr.-Ing. habil. Hermann Athen, Ebene und 
sphärische Trigonometrie (Bücher der 
Mathematik und Naturwissenschaften, herausgegeben 

. von Dr. Henry Poltz). 112 S. mit 51 Abb. Wolfen- 
büttel-Hannover 1948, Wolfenbütteler Verlagsanstalt 

- GmbH. (Notdruck.) Preis brosch. 7,50 DM. 
Das vorliegende Buch bringt in systematischer Dar- 
“ stellung in der Hauptsache das aus dem Gebiete der 

Trigonometrie, was für die Anwendungen wichtig ist. 
. Nachdenmi zunächst die Grundbegriffe und die Formeln 

der Goniometrie und ebenen Trigonometrie behandelt 

sind, werden in einem zweiten Abschnitt die Formeln 
der sphäriscnen Trigonometrie zunächst für Eulersche 
Dreiecke abgeleitet. Weiterhin wird dann der Drei- 
. ecksbegriff verallgemeinert und kurz auf Möbiussche 
Dreiecke und den Studysehen Dreiecksbegriff einge- 
' gangen. Hier'ergeben sich auch Beziehungen‘ zu an- 

deren Gebieten der Mathematik wie .der Gruppen- 

‘theorie, auf die verwiesen wird. Überall werden an 
' Beispielen die Anwendungsmöglichkeiten der For- 

meln gezeigt. Ganz den Anwendungen, nämlich denen 

in mathematischer Geographie und Astronomie, ist 

- der letzte Abschnitt gewidmet. | 
. Bei dem geringen Umfang des Bändchens kann na- 
türlich in Theorie und Anwendungen nur das Wesent- 
liche gebracht werden. Dieses ist im allgemeinen in 
‚ übersichtlicher klarer Darstellung gegeben. “Leider 
werden dem Anfänger einige Versehen und Druck- 
fehler auch in den Formeln den Gebrauch des Buches 
erschweren. Um nur einiges auf den ersten Seiten 
zu erwähnen. Der Abschnitt 1, der dem Buch von 

Komerell: Grenzgebiete der elementaren und höheren 

Mathematik entnommen ist, ist hier kaum verständ-. 

lich, da ein wesentlicher Gedanke fortgelassen ist; 

zweimal wird gesagt, daß der Umfang eines‘ Kreises 
von 1 m Radius nahezu 6243 mm (statt 6283 mm) ist; 
auf Seite 17 ist die"Lage des Vektors & falsch ange- 
geben; auf Seite 19 fehlt in der Formel für cos n« der 

Faktor cos? &— sin2«, usw. Bei einer hoffentlich 

bald. nötig werdenden Neuauflage werden sich diese 
kleinen Versehen leicht beseitigen lassen, 


Dresden. 'Willers. 


Dr.C. Schäfer (0. Prof. a. d. Universität Köln), 
Dr. L. Bergmann (Prof., Wetzlar) und Dr. W. Kli@foth 
(Prof., Heidenheim, Grundaufgaben des 


ümmungslinien. und geodätischen Linien 5 


Ott-HeinrichKeller. 2% 


Physikalischen Praktikums, (4. Aufl.), 
VI + 237-8. mit 180 Abb. und 18 Tabellen. Leipzig 
1948, B. G. Teubner. Verlagsgesellschaft. Preis brosch. 


5,80 DM. 


Der Führer durch das Physikalische Praktikum von 
Schaefer-Bergmann-Kliefoth behandelt die wichtigsten . 
Grundaufgaben des Praktikums aus allen Gebieten der 
Physik, ohne Bevorzugung einer Fachrichtung. Die 
Aufgaben, die besprochen werden — insgesamt 65 — 
habensich vielfach bewährt und finden sich mit kleinen 
Abwandlungen in jedem Praktikum. Das Buch kann 


also überall verwendet werden. Eine kleine Einleitung 


gibt einiges über Fehlerrechnung, ein Anhang enthält 
Tabellen von oft gebrauchten Größen. Die Verfasser 
haben angestrebt, daß .das Praktikumsbuch ohne 


zeitraubendes Nachschlagen in anderen Lehrbüchern 


‚hinreichend über die einzelnen Aufgaben unterrichtet. 


Zu diesem Zweck sind bei fast: allen Aufgaben in einem 
ersten Teil die theoretischen Grundlagen, in einem 
zweiten Teil Anordnung und Verlauf des Versuches 
besprochen. Die Darstellung ist ausdrücklich auf 
Anfänger: zugeschnitten, so daß das Buch auch von 


Chemikern, Ingenieuren, Medizinern mit Vorteil be- 


nutzt werden kann, die Physik nur als Nebenfach be- 
treiben. i } “ 


Dresden. : > A.Recknagel. 


Dr. F. Hund (0. Professor an der Universität Jena). = | 
„EinführungindietheoretischePhy- 
J s Al k“*, : 
Meyers kleine Handbücher, Bd. 45 —46. 303 8. mit | 
Bibliographisches_ Institut. 


BandI: Mechanik (zweite ergänzte Aufl.). 


128 Abb. Leipzig 1948, 
Preis brosch. 5,80 DM. 
Die’zweite Auflage der „„Mechanik“ von F.Hund 


-ist im wesentlichen gegenüber der ersten Auflage un- 


verändert geblieben. Das kleine Werk genügt in be- 
sonderem Ausmaß den Wünschen, die man an eine 


Einführung in die theoretische Mechanik stellt.- Der . 


äußere Umfang ist klein genug, so daß der Student 
nicht von vornherein abgeschreckt wird. ‘Der Preis ist 
erfreulich niedrig, ein Vorteil, den nur wenige unserer 
Lehrbücher haben. Dazu tritt die musterhafte Klar- 
heit der Darstellung und die Exaktheit,. die überall 


erreicht ist. Stößt man in einem Lehrbuch auf eine 


sachliche Ungenauigkeit: im.,,Hund‘ kann man sich 
die richtige Auskunft holen. Diese Zuverlässigkeit 
macht das Buch auch für den im Beruf stehenden 
Physiker zum treuen Ratgeber. 


Dresden. A.Recknagel. 

Dr.-Ing. Heinz Wittke (Regierungs-Vermessungs- 
rat der Katasterverwaltung und Lehrbeauftragter der 
Bergakademie Claustal), Vademeküm für Ver- 
messungstechnik. 334 S. mit"133 Abb. und 
24 Tafeln. Stuttgart 1948. J. B. Metzlersche Verlags- 
buchhandlung. Preis geb. 14,50 DM, 

Das kleine, handliche Taschenbuch bringt elementare 
Formeln, Regeln, Vordrucke, Zahlenbeispiele, Tabellen, 


Gedächtnisbrücken und Tafeln für Vermessungs- 


ingenieure, -Markscheider, Bauingenieure usw. Es 
wird insbesondere im Außendienst gute Dienste leisten. 


Dresden. Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen 


EINGEGANGENE BÜCHER . 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche. Besprechung bleibt vorbehalten) 


Dr. $. Flügge (o. Prof. an der Universität Marburg), 
Theoretische Optik (Bücher der Mathematik und 
Naturwissenschaften. Herausgegeben von Dr. Henry 


Poltz), 124 S. mit 31 Abb. Wolfenbüttel-Hannover 


1948. Wolfenbütteler‘ Verlagsanstalt G.m.b.H. (Not- 
druck). Preis brosch. 8, DM. 


Dr. Ing. Kurt Beyer (o. Prof. an der Technischen 
Hochschule Dresden), Die Statik im Stahl- 
betonbau. (Ein Lehr- und Handbuch der Bau- 
statik.) 2. Auflage. Berichtigter Neudruck. 804 8. 
mit 1372 Abb. Berlin 1948, Springer-Verlag. Preis: 
Halbleinen 66,— DM. 5 
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Dr. phil. L. Föppl (o. Prof. a. d. 
München), Are und er a 
keitslehre für i ‚III: ebene 

nungszustand. 192 8. m. 82 Abb. München 1947, 


Leibniz-Verlag. Preis geb. 15,— DM. 3 
"Dr. Wilhelm Schlink (Prof.an der Techn. Hochschule T: 


Darmstadt) und Dr. Ing. habil. Heinrich Dietz (Dozent 
an der Techn. Hochschule Darmstadt), 
sche Statik (Ein Lehrbuch zur Einführung ins 


NACHRICHTEN = 


Berlin: An der Techn. Universität Berlin wurden 
Dr. phil. Fritz Mühlig zum o, Professor für höhere 
Geodäsie, geodätische Astronomie und Geophysik, und 
Dr.-Ing. Isträn Szabd zum o. Professor für Me- 
chanik ernannt. 


Professor Dr. Helmut Hasse wurde zum o, 
Prof. an die Math.-Nat. Fakultät der Universität 
Berlin berufen, ; 


Dr. habil.HeinrichGrell wurde zum Professor 
mit Lehrauftrag an derselben Fakultät und zum Ober- 
assistenten am ersten Mathematischen Institut der 
Universität ernannt. 


Münster: Der Dozent Dr. phil.habil. Karl 
Stein wurde zum apl. Porfessor für Mathematik an 
der Universität Münster ernannt. - 


Dr. Ing. habil. W.Flügge wurde zum Professor 
an der Standford University in Californien USA. 


_ ernannt. r 


\ 
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Bemerkung zu der Arbeit K. Stange: Über 
die Verteilungsdichte der Meß- oder Beob- 
achtungsfehler eines dreidimensionalen Punkt- 
raumes. Z. angew. Math. Mech. 28, S.235 — 243. 

Die Beweisführung gilt unter der Voraussetzung, 
daß der Streuungstensor 7 nicht singulär, also seine 
Veterminante D, von Null verschieden ist. Diese 
Doraussetzung wird jedoch in der Arbeit ($. 235, Ab- 
schnitt 2)mit R > 3 unzureichend formuliert, worauf 
mich Herr Prof. Richter, Haltingen, hinwies. 
Die Bedingung R > 3 ist zwar notwendig, damit D, 
nicht verschwindet, aber keineswegs hinreichend. 
Man muß also ausdrücklich D,=FÜ voraussetzen. 
Der Fall eines singulären Streuungstensors 7 mit 
D, =0 bleibt demnach offen und wird von Herrn R. 
nächstens erörtert werden. 


Bad Sachsa. K.Stange. 


Zu H. Neuber: Vereinfachtes Verfahren zur 
Spannungsberechnung in dünnwandigen pris- 
matischen Hohlkörpern . unter Innendruck. 
Z. angew. Math. Mech. 28 (1948), S. 187—189. 


Herr H. Neuber teilt mit, daß ihm hinsichtlich 
seiner Mitteilung nachträglich noch eine französische 
Arbeit zur Kenntnis gelangte (Marbee: «Sur I’ Equilibre 
d’un anneau ferm6), L’Assoc. Techn. Maritime 1908), 
in welcher im wesentlichen das gleiche Problem er- 
örtert wird. Eine Darstellung der Ergebnisse jener 


Techni- 


 ZUSCHRIFTEN AN DEN HERAUSGEBER 


Eine höhere Festig- Abt 


Göttingen: Der 
der Uni : 
wurde zum 1.10. u nr 


Der Dozent Dr. Theodor Schneider 

am 15.7.1948 zum apl. Professor ernar a 

Clausthal: Der Leite schung 
Sterkrade, 


° G.H.H. Werk Dr. EberhardMe 
wurde am 1.4. 1948 zum o. Professor der 
und Mechanik an der Bergakademie C 


W. Schmeidler, Über die Wärm 
in einem-Körper. Z. angew. Math. Mec 
3.5459. HEEe 
In Formel (2a) sowie in allen fKormeln, 
in denen ein Vekto auftritt, ist das z durch 
ein x zu ersetzen, also z.B. [u x rotv]Jusw. 


Arbeit, insbesondere für den Spezialfall 
metrischer Querschnitte findet sich bei: 
Math. Phys. 61 (1913), S. 246 ff. ae 


Zu C. Weber: Uber elliptische Integrale. 
Z. angew. Math. Mech. 28 (1948), .26. 
Die angegebenen Beziehungen (4) und (6) sind 
Spezialfälle des Additionstheorems für elliptische Inte- 3 
grale. Dies sagt aus, daß 
F(p,k)+F(y,k)=F(e,k) 
und 
.E@k)+E(y,k)=E(o,k) + k’sinpsiny sin o 
wenn 
sino = 
_ sing cosy J1 -ktsintp+sinypcospy 1 -ktsin®g 
1 — k?sin’o sin? y 
oder, was damit gleichbedeutend ist, 
COS 0 = COS P cos y — sing sing 1—k?sinto. S 
Für o= 5 erhält marrdie oben angegebenen Spezial- 
fälle. (Vgl. O.Schlömileh, Comp, d. höh, An- 
alysis, Bd. II (1865), S. 324 ff.) j 
Berlin-Oberschöneweide. Klaus Kirmes, 


- Prof. Weber hat von dieser Zuschrift dankend 
Kenntnis genommen. - 
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